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ABSTRACT 

We present Spitzer Space Telescope IRAC and MIPS observations of a 0.85 
deg^ field including the Corona Australis (CrA) star-forming region. At a distance 
of 130 pc, CrA is one of the closest regions known to be actively forming stars, 
particularly within its embedded association, the Coronet. Using the Spitzer 
data, we identify 51 young stellar objects (YSOs) in CrA which include sources in 
the well-studied Coronet cluster as well as distributed throughout the molecular 
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cloud. Twelve of the YSOs discussed are new candidates, one of which is located 
in the Coronet. Known YSOs retrieved from the literature are also added to the 
list, and a total of 116 candidate YSOs in CrA are compiled. Based on these 
YSO candidates, the star formation rate is computed to be 12 MgMyr"^, similar 
to that of the Lupus clouds. A clustering analysis was also performed, finding 
that the main cluster core, consisting of 68 members, is elongated (having an 
aspect ratio of 2.36), with a circular radius of 0.59 pc and mean surface density 
of 150 pc-2. 

In addition, we analyze outflows and jets in CrA by means of new CO and 
H2 data. We present 1.3 mm interferometric continuum observations made with 
the Submillimeter Array (SMA) covering RCrA, IRS 5, IRS 7, and IRAS 18595- 
3712 (IRAS 32). We also present multi-epoch H2 maps and detect jets and out- 
flows, study their proper motions, and identify exciting sources. The Spitzer 
and ISAAC /VLT observations of IRAS 32 show a bipolar precessing jet, which 
drives a CO (2-1) outflow detected in the SMA observations. There is also clear 
evidence for a parsec-scale precessing outflow, E-W oriented, and originating in 
the SMA 2 region, likely driven by SMA 2 or IRS 7A. 

Subject headings: infrared: general — open clusters and associations: individual 
(Corona Australis, CrA) — stars: formation 



1. Introduction 



The Gould Belt Spitzer Legac y program is a GO-4 pr ogram designed to extend the 
earlier Spitzer Cores to Disks (c2d; lEvans et al.l l2003l 120091 ) program, thereby completing 
a census of star-forming regions within 500 pc. The Gould Belt is a band of stars ari d 
molecular clouds located within ~20° of the Galactic Plane ( iHerschell 118471 : iGouldl Il879l ). 



Although, at a declination of approximately —40°, Corona Australis (CrA) is not technically 
within the Gould Belt, it is discussed in this paper as part of the survey. We present an 
extensive study of the entire Corona Australis star-forming region, investigating the overall 
young population, its spatial distribution along the molecular cloud, and the young stellar 
object (YSO) outflows. 



Rossand ( 119781 ) used star counts to create an extinction map, and identified five clouds in 



CrA, named clouds A-E, noting that CrA is highly elongated, oriented nearly east to west in 
the sky. Cloud A is in the west, and it is th i s clou d which corresponds to the RCr A/Coronet 
region. Large-scale CO mapping by iLoreru (jl979[ ) also shows this elongation, but thei r high 



spatial resolution observations of the velocity field indicated no velocity gradient. iLoren 
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( Il979l ) interpreted this to mean that t he el ongation is not d ue to c ontra ction along the 
rotational axis. However, iHarju et al.l ( 119931 ) argued that the iLorenl (119791 ) data were not 
sufficient, and suggest from their own C^^O observations that the RCrA core is, in fact, a 
fragmented disk, explaining the observed elongation along the major axis. It is somewhat 
surprisi ng that CrA is not associated with the Gould Belt, considering its close distance of 
130 pc0; iMamajek fc Feigelson ( 2001 ) argue that instead it for med as part o f exp anding Sco- 
Cen superbubbles, specifically Loop I, citi ng evidence from thelHariu et al.l (119931) millimeter 
observations, and radio observations from lCappa de Nicolau fc Poppell ( 1l99ll ). 



Studies from the literature have mainly focused on RCrA, the brightest star in the clus- 
ter, the Coronet region, and its population. The variability of the nebula surrounding the 
Herb ig Ae star, R CrA, has been known since the early 1900s (IKnox Shawl Il916l : iRevnolds 
Many years later, the two variable stars R CrA and T CrA were identified by iHerbig 



(Il960l ) to be young, and he then concluded that the associated stars should also be young. 
This prompted an interest in studying the RCrA region, and in 1973, the ffist major op- 
tical and infrared study of the main stars near RCrA was conducted, finding a total of 11 
stars in the young stellar group: TYCr A, S CrA, TCrA, R CrA, DC CrA, VVCrA, KS-15, 



HR7169, HR 7170, Anon l, and Anon 2 (IKnacke et al.l 119731 ). The following year, IRS 1 was 



suggested by IStrom et al.l (119741 ) to be the driving source for the Herbig- Haro (HH) object. 
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Wilking et al. 1992. who ffist mentioned 



IRAS 32). 



Neuhauser fc ForbrichI (120081 ) have recently reviewed the literature on the entire Corona 
Australis star-forming region, although many of the studies cover only a subset of the region 
which we include here as part of our Spitzer IR AC and MIPS study . The most releva nt re- 
cent studies include: deep infrared obser vations (IWi 



limeter and submillimeter observations (IChini et al 



dng et al 



2003 



1997 



Haas et a. 



20051 ) . including Submillimeter Array (SMA) observations (iGroppi et al 



iGroppi et al. 



observations with the Infrared Space Observatory (ISO) (lOlofsson et al 



200. 



ries of papers focusing on Chandra X-ray studies of the region (IForbrich et al. 



2008), mil- 



Nutter et al. 



2007D. m id- infrared 



1999 ) , and a se 



2006, 2007 



^The distance to CrA is based on measurements to several st ars, ranging; from 85 to 190 p c. A compre- 
hensive summary of the various measurements can be found in iNeuhauser &: ForbrichI (|2008[ ): throughout 
this paper we use their suggested distance of 130 pc. 



-A- 



Forbrich &: PreibischI 120071 ). F inally, some spectroscopic work has been done to deter- 



mine association memb erships ( jPattenl Il998l : iNisini et al.l l2005l : ISicilia-Aguilar et al.l 12008 
Mever fc WilkineJboOoh . 



The most massive st ars in CrA are the Herbig Ae/Be stars, RCrA and TYCrA. RCrA 
has a spectral type of A5 (IKnacke et al.lll973l : iMarraco fc Rydgrerull98ll ) and is located at the 
tip of a cometary shaped reflection nebula, N GC6729. TYCrA , located t o the northwest 



of RCrA, is at least a quadruple system (ICasey et al. 



the primary has a spectral typ e of B8-B9 (IHerbig fc Kameswara Rao 



1973 



1995 



IChauvin et al. 12003) where 



1972 



Knacke et al. 



Marraco fc Rydgrenlll98ll ). and is associated with the reflection nebula NGC 6726/7. 
Located 1' south of TYCrA, is HP 176386, which is also surround ed by reflection nebulosity 
from NGC 6727 JXnacke et al.l Il973l: iMarraco fc RvdgrenI Il98ll ). HD 176386 is a binary 
system, first recognized as such by IWilking et al.l (119971 ): HD 17 6386A/B is a visual pai r 
separated by 3. "7 with spectral types of AOV and K7, respectively (IMeyer fc Wilkingll2009l ). 
Additionally, a little further from the Coronet (~12' to the southwest of RCrA) lie two 
B8V stars, HR7169 and HR7170, which are discussed in more detail in App endix I A. 1.111 
Evidence for heating of the molecular cloud by these two B8 stars has been seen (jLorerull979l ). 
making it likely that there is a physical association and that they are there fore located at 
the same distance as RCrA (see further discussion in lNeuhauser et al.ll2000l ). 



CrA is also well-known to harbor many active YSOs with outflows. To date, twenty HH 



Eight obiects (HH 82, 96-101, 104) were detected in the 1970s and 1980s ( 


Strom et al. 


1974; 


Schwartz et al. 


1984; 


Hartiean & Graham 


1987 


; ] 


-leipurth & Graham 


19881), while the 


remaining twelve objects were observed by 


Wang et al. 


(2004 


) in a large optical survey. The 



majority of the detected objects are located close to the Coronet and seem to be driven by 
YSOs inside it or in its outskirts. A few more HH objects are positioned close to HHIOO-IR 
(IRSl), SCrA, VVCrA and IRAS 18595-3712 (IRAS 32), which seem to drive outflows as 
well. On the other hand, there ar e very few and sparse studies oii H2 jets and outflows in 
CrA (see e. g. Iwilking et aP Il99d ICredell Il994l : IPavis et al.l Il999l : ICaratti o Garatti et al. 
20061 ). In these papers, the H2 counterparts of HH99, 100, 101 and 104 were identified and 



studied, and a few new jets in the Coronet were detected (ICaratti o Garatti et al.ll2006l ). 
Davis et al.l feoiol ) cataloged five molecular hydrogen objects (MHOs, MHO 2000-2004), but 
so far, an extensive H2 map of the region has not been made. Thus, a complete census of 
outflows and their driving sources is lacking. 

We present S'jozfeer observations of a 0.85 deg^ region in the Corona Australis molecular 
cloud, identifying the Spitzer-selected YSOs distributed throughout the cloud, as well as 
the outflows and their driving sources. This study includes infrared imaging of a much 
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larger portion of the molecular cloud than many previous studies have included. In § [2], we 
discuss the Spitzer observations and data reduction, including basic statistics for the sources 
detected, as well as the ancillary SMA observations, and H2 observations of outflows. In 
§ El we discuss methods used to select YSO candidates from color-magnitude diagrams 
constructed from the Spitzer data, as well as other methods of selection. We also discuss 
the addition of known YSOs and YSO candidates from the literature because they were not 
classified from Spitzer data due to saturation or other observational issues. The distribution 
of YSOs is discussed in § SJ There we present an extinction map created from near-infrared 
and Spitzer data, an analysis of the spatial distribution of the YSOs in the cloud, and the 
clustering analysis performed using the YSO candidates in CrA. In §[5|, we present an analysis 
of the SMA observations. In § [6l we present an analysis of the jets and outflows, the proper 
motions of the H2 knots, and the driving sources for the outflows and jets seen in CrA. 
Finally, in § [TJ we discuss the overall cloud properties and how they compare with the other, 
c2d and Gould Belt regions, and summarize the results in § [HI 



Observations and Data Reduction 



2.1. Spitzer IRAC and MIPS 



Corona Australis was observed with the Spitzer Inhaied Array Camera (IRAC: lFazio et al. 



2004) at 3.6, 4.5, 5.8, and 8.0 ^m and the Multiband Imaging Photometer for Spitzer (MIPS; 
Rieke et al .1120041 ) at 24, 70, and 160 /im as part of two guaranteed time observation (GTO) 
programs (PID 6, 30784; PI Fazio) as well as the Gould Belt Legacy Survey (PID 30574; PI 
Allen). Tabled] summarizes the program identification numbers, AOR identification numbers 
and observation dates for all of the CrA data included in this paper. The IRAC mapping 
includes one epoch, two dithers, with 12 second integration times per exposure, in high- 
dynamic range mode. The MIPS mapping was done with medium scan rate, one epoch, 
with 148" return and forward leg cross scan steps. The size of the overlap area from IRAC 
band 1 - MIPS band 1 is 0.85 deg^; Figure [1] shows the IRAC and MIPS coverage in CrA. 

The basic calibrated data (BCD) are downloaded from the Spitzer archive and pro- 
cessed by the team using the same, custom pipeline p rocessing prograna s and techniques as 
used by c2d (see the c2d Final Delivery Document, lEvans et al.l 120071 ). Briefly, the data 
were inspected, custom masks were created to identify bad pixels, and corrected for in- 
strur nental effects. Mosaic s were created from the improved data using the MOPEX pack- 
age (IMakovoz et al.l l2006l ) . and source extraction was performed using the c2dphot tool 
( lEvans et al.ll2007| ). which is a derivative of DoPHOT (ISchechter et al.lll993l ). Source lists 
for detections at each wavelength were band-merged together with the 2MASS catalog 
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( jSkrutskie et al.l 120061 ) and cross i dentifications are a ccurate within 2" (see Section 2.4 of 
the c2d Final Dehvery Document, lEvans et aLl 120071 ). Also note that the final catalog is 
"band-filled" (see e.g., the c2d Final Delivery Document) to produce fiux estimates of ob- 
jects that were not found in the original source extraction processing, but were detected in 
the other bands (see further discussion in §[3]). This procedure is described in de tail in the 
delivery documentation for the c2d Final Delivery Document (lEvans et al.l 120071 ). Briefiy, 
c2dphot fits the PSF profile at the position of the known source. Note, however, that all these 
fiux densities, fiagged as band-filled in our tables, should be considered as "bad photometry," 
in the same way as the Sa upper limits in this catalog. 

Statistics for the sources detected in CrA with S/N of at least 3 can be found in Table 
[2J This corresponds to selecting all sources with detection quality "A," "B," or "C" (S/N 
> 7, 5, and 3, respectively) in any of the IRAC bands from the final delivered catalogs (cf. 
the c2d Final Delivery Document). 

Figures |2] and [3] show color mosaics of the entire CrA region (mapped with IRAC and 
MIPS), using IRAC bands 2 (4.5 fim) and 4 (8.0 /um), and MIPS 24 /im (blue, green, and red, 
respectively). Figure [2] includes the Coronet region, the SE and SW filaments, and Figure [3] 
shows the less known region called the "streamer," positioned further to the southwest. 



2.2. Ancillary Submillimeter Array Observations 



The regio ns around the previou s^ 



2006), IRS 7 



Wilking et al 



Wilking et al. 



1997 



19921: IConnellev et al 



V identified YSQs I RS 5 flWilking et al.lll997l : iForbrich et al 
Groppi et al.l l2007h. and IRAS 32 (IRAS 18595-3712; 



jropp^ 
2007t 



Ivan Kempen et al.l 120091 1 were observed with the 
smaI^ (IHo et al.ll2004l ) m the dust continuum near 225 GHz (~ 1.3 mm). Observations of the 



IRS 7 region were made on 2006 August 20 (program 2006-03-S046), while IRAS 32 and IRS 5 
were observed on 2008 June 10 and 14 (program 2008A-S074). In this paper, we present 
only the continuum data for IRS 7, which were obtained as part of a line survey program; 
the results from that program will be reported elsewhere (Lindberg et al., in preparation). 
The data were edited and calibrated in the IDL-based software package MIR|f| 



For the observations of IRS 5 and IRAS 32, the SMA was in its compact-north config- 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 



■^http: / / cfa-www.harvard.edu/~cqi/mircook.html 
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uration, with eight antennas, covering basehnes of 16-139 m. The receiver's two sidebands 
were centered near 221 and 231 GHz. Chunks containing obvious spectral hues were removed 
before forming the continuum, resulting in a total continuum coverage of ~ 3.6 GHz centered 
at 226 GHz. Observations of the quasar 3C 279 were used for bandpass calibration, and the 
quasars J1924-292 and J1937-399 were used for complex gain calibration. 3C 279 was used 
for absolute flux calibration, assuming a flux of 10.3 Jy based on observations at dates near 
in time to these that were flux calibrated using planets and moons0 

For the observations of the IRS 7 region, the SMA was in its compact configuration 
with six antennas, covering baselines of 16-69 m. A two-point mosaic pattern was used to 
cover the region. The receiver's two sidebands were centered near 219 and 229 GHz, with 
each sideband consisting of 24 overlapping chunks of 109 MHz each, or about 2 GHz for 
each sideband. Spectral lines were removed from the data, resulting in a total continuum 
coverage of ~ 3.3 GHz centered at 224 GHz. Observations of the quasar 3C 454.3 were used 
for bandpass calibration, and the quasars J1924-292 and J1957-387 were used for complex 
gain calibration. Uranus was used for absolute flux calibration. 

The MIRIAD software package was used for Fourier inversion of the visibilities, CLEAN 
deconvolution, and restoration with a synthesized beam. A Briggs robust weighting of was 
used when mapping the continuum emission. The data were scaled by the inverse response 
of the primary beam ("primary-beam corrected"), to account for the loss of sensitivity away 
from the phase center. Synthesised beam sizes and 1 a rms sensitivities are: 5'.'5 x 2'.'3 and 
7.4 mJy for IRS 7, 4'.'6 x 2'.'6 and 3.7 mJy for IRS 5, and 3'.'7 x 2'.'2 and 4.4 mJy for IRAS 32, 
respectively. 



2.3. Ancillary Ha Data 



As a complement to the Spitzer an.d SMA data, we make use of multi-epoch H2 (2.12 yum) 
maps to detect jets and outflows in the CrA star-forming region to study their proper motions 
(P.M.s) and identify the exciting sources. Most of the data collected also have or narrow- 
band filter images to remove the continuum from the H2 maps. The relevant information on 
these ancillary data is reported in Table [3l 

Our narrow-band H9 im age archive is composed of dat a collected at ESO-NTT wit h 
SofI ( iMoorwood et al.lll998al ) in 1999 (already published in ICaratti o Garatti et al.ll2006l ). 



http://smal.sma. hawaii.edu/callist/callist.html?data=1256-057 
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and in 2007. Additional images were retrieved from the ESO science archive facihtjif] taken 



at the ESO-VLT with ISAAC flMoorwood et aLlll998bf ) in 2000, 2003, and 2005 



SofI 2007 images were observed in service mode between June and August 2007, and 
cover four regions of CrA (~10'xlO' each), mapping in total an area of ~20'xl7' (about one 
tenth of the area mapped by Spitzer), for a total integration time of 5760 s. Figure H] shows 
the area covered by the H2 mosaic. 

Earlier epoch maps enclose much smaller regions, centered on RCrA {ISAAC 2005 
and SofI 1999, ~5'x7' and 5.2'x5.2' FoV, respectively), on HH99 {ISAAC 2000, ~6.8'x6') 
or HHlOl {ISAAC 2003, 6.3'x6.3'), the details of which can be found in Table H These 
areas are indicated in Figure HI where different polygons are superimposed on the SofI 2007 
mosaic, showing regions mapped in different epochs (yellow: 1999, cyan: 2000, green: 2003, 
magenta: 2005). As a consequence, the proper motion (P.M.) analysis was mostly performed 
on knots inside and close to the Coronet region. 

Finally, additional ISAAC images (H2 and K, filters) around IRAS 18595-3712 (IRAS 32; 
located outside the SofI 2007 map), taken in 2005, were acquired to detect a possible H2 jet 
from that source. 

All the raw data were reduced using IRAF packages applying standard procedures for 
sky subtraction, dome flat-fielding, bad pixel and cosmic ray removal, and making image 
mosaics. Continuum-subtracted images, B.2-Ks and H2-NB(2.09/um), were obtained by sub- 
tracting Kg and NB (2.09 /im) images, appropriately scaled and registered. Each pair of 
images was matched by means of tens to hundreds of field stars (depending on the mosaic 
size and crowding), selected excluding the YSO candidates. The scaling has been done by 
performing relative photometry on the selected stars. The images were not flux calibrated. 



3. YSO Selection and Classification 



There are 45 sources in the observed Corona Australis region which were initially clas- 
sified as Y SO candidate s frorn their Spitzer plus near- infrared colors using the technique 
outlined in iHarvey et al.l (2007, 2008, and references therein). This technique has been suc- 
cessfully adopted for the c2d and Gould Belt Spitzer surveys. Indeed, 33 out of the 45 sources 
in CrA (~73%) have already been confirmed as YSOs in other studies (see § 13. H and Appendix 
|A|) . Briefly, the selection method uses criteria based on a combination of infrared excess with 
a brightness limit, in order to minimize extragalactic contamination. The selection method 



'''http:/ /archive. eso.org/ 
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requires a S/N of 3 or higher detection in all four IRAC bands as well as in the MIPS 24 fim 
band. Figure [5] shows the color-magnitude space used to separate the YSO population in 
CrA from contaminants, which are mainly background giant stars and galaxies. To illustrate 
where the contaminating galaxy population is located in the same color-magnitude space, 
Figure [5] also shows color- magnitude diagrams for t he full and re-samp led Spitzer Wide-area 
InfraRed Extragalactic Survey (SWIRE) catalogs (ISurace et al.ll2004l ). The SWIRE obser- 
vations were processed in exactly the same way as our data were processed in order to be 
used for a direct comparison (except that band-filling was not performed). Note that the 
SWIRE data are considerably deeper than ours, so the SWIRE catalog was trimmed down 
assuming the SWIRE observations had been obtained with sensitivities similar to those of 
our observations (this results in the "re-sampled SWIRE" catalog). A f ull discussion of th e 
specifics of this process can be found in the c2d Final Delivery Document (lEvans et al.ll2007l ). 



Table m lists all YSO candidates sel ected through this method in CrA along with their 
Spitzer fluxes, known names, and classes (lLadalll987l : lGreene et al.lll994j ). Class is determined 
from the spectral slope, a, which is calculated over the widest range possible where data are 
available between 2.2 and 24 /im, and is defined as: 



a 



d log(AF,) 
d\og{X) ' 



The value of a is calculated from a linear fit to the logarithms, taking into account the uncer- 
tainties in the flux measurements. We classify sources with a > 0.3 as Class I, —0.3 < a < 0.3 
as Flat spectrum, —1.6 < a < —0.3 as Class II, and a < —1.6 as Class III. Note that deeply 
embedded protostars. Class sources, which require submillimeter data for identification, 
cannot be distinguished from Class I objects in this analysis. 

It is important to note that the a quoted in the tables and in our catalogs is computed 
by selecting all valid wavelengths (having a quality of detection not "N" or "U" , meaning a 
non-detection or an upper limit, respectively) in the combined epochs between and 24 fim 
(up to 6 possible bands). Any fluxes or their uncertainties that are NaN, zero, or negative 
are excluded. In some cases, for the calculation of a, a band-filled flux has been used (in 
the c2d catalogs this is indicated as having an image type of "-2"). Band-filling typically 
happens for sources that have clear shorter wavelength IR AC fluxes, but are undetected at 
longer wavelengths (see the c2d Final Delivery Document, lEvans et al.ll2007l ). If a source in 
Table H] has a band-filled flux, it will be indicated as such with a footnote. 

The use of the terms "YSO" and "YSO candidate" will be somewhat interchanged 
throughout this paper. To be clear, not all the sources discussed as YSOs are spectroscop- 
ically identified YSOs or confirmed as members of the CrA cloud, although some of them 
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have, through other studies. This is why the term "YSO candidate" is used most frequently 
when describing specific sources. Ample references to the literature are provided in the text, 
appendices, and the tables in order to facilitate further studies of these YSOs and YSO 
candidates. 



3.1. Spitzer Classified YSO Candidates 

There are 45 YSO candidates classified using the standard technique, as discussed above, 
and are listed in Table ID In two cases (CrA-2 and CrA-39), the candidates turned out to 
be known galaxies; they are discussed in Appendix |X] and appear in the table, although are 
not counted in our final list of YSO candidates. There are also three YSO candidates which, 
upon visual inspection, were determined not to be YSOs: CrA-17, CrA-25 (which is a likely 
HH object), and CrA-32 (see descriptions for all in Appendix |A]) . This gives a final total 
of 40 YSO candidates selected by our method with Spitzer. Many of the YSO candidates 
are already known YSOs from the literature, but there are 7 new YSO candidates (CrA-1, 
CrA-7, CrA-9, CrA-22, CrA-24, CrA-36, CrA-37) that have not been selected by any other 
survey, in most cases because the Spitzer observations extend beyond the region included in 
other surveys. CrA-24 is a new YSO candidate which is located within the Coronet, and 
further discussion of it can be found in lA. 1.241 For a more detailed description of each of 
the 45 sources selected using this method, see Appendix \^ 

The spectral energy distributions (SEDs) for the 45 sources are shown in Figure [6] (Class I 
and Flat spectrum candidates). Figure [7] (Class II candidates), and Figure [8] (Class III candi- 
dates). In these figures, the open circles represent the Spitzer and ancil lary data, and the filled 



circle s represent the dereddened fiuxes (using the extinction law of IWeingartner fc Draine 



200ll . with Ry =5.5). In Figures [7] and [H a grey line represents the photosphere of a K7 



main seque nce star, and the das hed black line represents the average SED for T Tauri stars 



in Taurus (ID'Alessio et al.lll999l ). 



3.2. YSO Candidates Classified from 2MASS and MIPS 



In some cases, there are regions covered by MIPS that are not covered in IRAC (or else 
are only covered in IRAC bands 1/3 and not in 2/4, and vice versa). To have as complete a 
sample as possible, we use 2MASS + the MIPS data in order to select additional YSO candi- 
dates. A vs . K., - [24] co l or-ma gnitude diagram for CrA is shown in Figure M Following 
the analysis of iRebuU et al.l (120071 ) for Perseus and lPadgett et al.l (120081 ) for Ophiuchus, we 
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find several new YSO candidates based on this diagram. The specific criteria used to select 
YSO candidates in this way are as follows: < 14 magnitudes, 24 /xm < 9 magnitudes, 
and Kg - [24] > 2. The region of the color-magnitude diagram which meets these criteria 
is outlined in Figure M There are 12 YSO candidates which meet these criteria, and they 
are listed in Table [51 However, as was for the case of CrA-2, source CrA-50 is part of the 
extended galaxy Leda 90315, and so it has been rejected from the hst of candidates (see also 
§ IA.1.21 for a discussion of CrA-2). Thus we yield a total of 11 new YSO candidates from 
this technique, 5 of which have never been identified as candidate YSOs before. 



For the sources discussed in this section and in Section 13.11 the final list of Spitzer- 
selected sources includes 51 YSO candidates, 12 of which are new. Using the values of a 
defined at the beginning of Section [31 we find 10 Class I or Flat spectrum, 35 Class II, 
and 6 Class III candidate YSOs. This high fraction of Class II sources compared with 
other classes, is co nsistent with that seen in the previously studied c2d star- forming regions 



flEvans et al.ll2009[ ). 



3.3. Class III Population in CrA 

While infrared surveys are excellent for identifying embedded YSOs still surrounded 
by circumstellar dust, they are less successful for classifying the Class III YSO population. 
Class III objects do not have the significant infrared excess associated with optically thick 
disks, and thus cannot be clearly distinguished from main sequence stars. That is apparent 
from the analysis discussed in § 13.11 (see Table [4]), where only 6 Class III YSO candidates 
are identified. On the other hand, diskless Class III objects can be readily identified using 
X-ray observations in addition to the infrared (e.g. Grosso et al.lboool ). Therefore in order to 



obtain a census of the young population in CrA, we combine our Spitzer catalog with X-ray 
ROSAT data (Patten et al., in preparation) to select additional Class III YSO candidates. 
The completeness of our selected sample of Class III YSOs will be discussed in § 13. 5[ 

There are 34 ROSAT detections within our Spitzer CrA field. Of those, 3 do not have 
a Spitzer counterpart; two sources are off the main cloud and have Spitzer detections on the 
border of the positional uncertainty circle, so we do not make the association. The third 
(which is located within the high nebulosity region of TY CrA) is likely to be detected, but 
there are two Spitzer detections within the positional uncertainty, making the association 
ambiguous. Three ROSAT detections are the well-known YSOs in CrA known as IRS 2, 
S CrA, and TYCrA. Six are YSO candidates which we classified as Class II or III using 
Spitzer and 2MASS: CrA-6, 11, 14, 16, 30, and 40 (see Table ^. Two Class II sources, 
CrA-52 and CrA-54 (listed in Table [S]), were selected using and 24 /im, as described in 
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§ 13.21 The remaining 20 sources are listed in Table [6] as YSO candidates and are included in 
the clustering analysis later in this paper (see § 14. 2p . All 20 of these sources are classified as 
"star" in our Spitzer catalogs (meaning they have SEDs that are well-fitted by a reddened 
stellar photosphere) and have a values within the range of a Class III YSO, as calculated 
from their Spitzer + 2MASS SEDs. 



3.4. Known YSOs Not Classified with Spitzer 



The Coronet is a well-known young stellar cluster, with several YSOs previously studied 
by several authors at visible and near-infrared wavelengths. Unfortunately this region is very 
bright in the longer wavelength Spitzer bands and many sources are saturated in our images, 
which may result in a band-filled fiux. Therefore we do not have enough SED coverage to 
classify them as YSO candidates (remember from § |3]that we require a S/N of 3 or higher 
detection in all four IRAC bands as well as in the MIPS 24 /im band before even considering 
a source as a YSO candidate). Nevertheless, since we want to analyze the entire CrA 
population, these sources should be included in a comprehensive list of YSO candidates. We 
reviewed the literature to recover those YSOs (and YSO candidates) which are not already 
included in Tables HI [5l and [61 Our findings are listed in Table [3 

In Table [TJ we summarize the known YSOs along with the coordinates and fiuxes of 
their Spitzer counterparts as well as information on how they are classified in the delivered 
catalog. Many of the most massive stars in CrA were added in this way, including RCrA. 
In the case of R CrA, because many of the Spitzer bands are saturated, a is computed from 
the K- and 5.8 /xm bands alone, resulting in a classification of Class III. However, it should 
be noted that RCr A has long been kn own to be a pre-main sequence ( PMS) star with an 



accre tion disk (e.g. iKnacke et al.lll973l ). and an SED slope of a Class II (IHillenbrand et al. 



19921 ). In addition, there are three known YSOs in CrA which do not have a counterpart 
in our Spitzer catalog: IRS 9, SMA2, and Anon 2 (also known as VSSt 14 and HBC675). 
In the case of IRS 9 and SMA2, the nebulosity from RCrA is too bright to distinguish the 
sources, and in the case of Anon 2, it is a diffraction spike fr om S CrA that causes a pro blem. 
IRS 9 is classified as a Class I (an d detected with C/t andra: iForbrich fc PreibischI 120071 ). and 



SMA2 is class ified as a Clas s (IGroppi et al 



debated: both iPattenI (119981 ) and 



2007 



Neuhauser et al 



Whether or not Anon 2 is a YSO is 



(12000) clai m it is not a young member 



(lithium is not detected in its spectrum. iNeuhauser et al.l 120001 ). Unfortunately without the 
Spitzer fiuxes available, we cannot determine a classification for Anon 2. 

Not all sources that have ever been identified as a candidate YSO are included in Table [3 
For example, IRS 3 and IRS 4 are not included because they are spectroscopically recognized 
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as background giants (IMeyer fc WilkingI |2009| ). Additionally, any candidate YSOs in the 
literature which our Spitzer observations classify as a galaxy (noted as "Gale" in the catalog) 
are not included. There is an exception to this, and that is when additional observations 
show that a source classified as "Gale" in our catalog should be a member, as in the case of 
LS-RCrA 1 (discussed at the end of this section). 



Moreover, we compare our candidates to the iHaas et al.l ( l2008l ) list (their Table 3), 
which was selected using a millimeter excess technique. Of those 38 candidates, 8 are al- 
ready included in our tables. There are 5 sources in their list which are classified in our 
Spitzer GB catalogs with SED slopes described as "star" or "cup-down," and their fiuxes 
are reasonable enough (i.e. are not band-filled and have good quality detections) to include 
them as candidate YSOs in Table [71 Note that a "cup-down" classification means that the 
SED has a shape that is c onvex and can not be classified by any other category (see the c2d 
Final Delivery Document, lEvans et al.l 120071 ). However, many of the remaining sources are 
listed in our catalog as "one" or "two," meaning they are only detected in one or two of the 
IRAC bands (generally in the 3.6 /im and/or 4.5 /im bands). Because there is such little 
information in the Spitzer bands for a proper classification, none of these sources has been 
included in our Table [71 

Finally, Table [71 also reports 2 out of 35 sources (namely G 09 CrA-9, 11; from their 
Table 4) classified as candidate YSOs i n iGutermuth et al.l (120091 ) that are not included in 
our previous tables. IGutermuth et al.l ( 120091 ) used Spitzer data of CrA available from the 
GTO PID 6 program (see Table Hj), which is included in the area studied in this paper, but 
covers a much smaller field, and adopted a slightly different selection technique. They found 
35 YSO candidates in the region covered by PID 6, and of those, 25 are also selected by our 
method to be YSO candidates. Two of their sources have been added to our Table [71 G09 
CrA-11, which we classify as "star" and G09 CrA-9, which we classify as "cup-down." Of 
the remaining 8 candidates, 4 are known YSOs which already appear in Table [71 (IRS 7A, 
IRS7B, SCrA, and TCrA), and 4 are classified as "Gale" in our catalog. However, one o f 
the "Gale" sources is a star of late spectral type, LS-RCrA 1 (iFernandez fc Comer6nll200ll ). 
and thus has also been included in Table [71 

Table [71 along with Tables HI [3 and El are intended to make up the most comprehensive 
list of YSOs and YSO candidates in CrA to date. The sources included in these tables 
(minus those sources identified as likely galaxies) will be used to analyze clustering in CrA 
(§ 14. 2p and to determine the driving sources for the outfiows/jets (§ 16. 3p . 

The total number of candidate YSOs in CrA is 116, where 14 are classified as Class I, 5 
are Flat spectrum, 43 are Class II, and 54 are Class HI. All 116 sources are shown overlaid 
on a 24 /im image in Figure [TOl 
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3.5. Completeness 



Evans et al.l (120071 l2009l ) estimated that the c2d (and similarly, GB) selected sample of 
YSO candidates is 90% complete down to a luminosity of O.O5L0, integrated from 1 to 30 /xm 
at a distance of the Serpens star- forming region (260 pc). However, the observing strategy for 
CrA is slightly different from the rest of the GB regions because it was originally observed 
as part of a GTO program. There is only one epoch, as opposed to the usual two, and 
therefore the photometric depth is ~30% less. Accounting for the difference in photometric 
depth and scaling for the distance of CrA (130 pc), we estimate that our 2MASS/Sj>zfeer 
sample, selected using the c2d/GB method (discussed in § 13.11 and reported in Table S]), is 
~90% complete down to ^O.OISL©. 

To illustrate the completeness of the CrA sample, we constructed bolometric lumi- 



nosity functions f or the Class 11 and 111 candidates (Figure [TTl see also iMerm et al.l 12008 



Spezzi et al.l 1201 ll ). First, we estimate the bolometric luminosity of the YSO candidates by 
integrating over all the observed SED data points, assuming a distance of 130 pc to CrA. 
The solid line in Figure [11] indicates the bolometric luminosity distribution of Class 11 and 
III YSO candidates discus sed in ^ 13.11 (and reported in Table Hj) rebinned to a 0.03 size 



bin. Then, following the iHarvey et al.l ( 1200 7l ) method, we apply the c2d completeness cor- 
rection to the samp le of sources selected using the c2d/GB method discussed in § 13. 1[ The 
Harvey et al.l (120071 ) completeness correction compares, for each luminosity bin, the number 
of counts from a trimmed version of the deep SWIRE catalog of extragalactic sources with 
the number of counts for the c2d catalogs in Serpens. The completeness correction estimate 
is then applied to the five molecular clouds observed by c2d. 

For all the regions in the GB survey, this completeness correction translates, since they 
are homogeneous in terms of photometric depth. However, as stated previously, this is not 
precisely the case for CrA. If we assume that the c2d/GB survey samples are homogeneous 
with CrA down to ^O.OISLq, the completeness correction, indicated by dashed lines in 
Figure [TTl provides a completeness estimate of ~90%. We would like to emphasize here 
that CrA is peculiar when compared with the c2d and GB star-forming regions previously 
studied, in terms of the initial observing strategy, and mostly because of various observational 
limitations (discussed in detail in the following paragraphs), which led to the inclusion of 
many additional YSO candidates, selected using different techniques (see § 13. 2^ § 13. 3^ § 13. 4p . 
To address the inclusion of these sources quantitatively, the bolometric luminosity function 
for all sources is overplotted in Figure ITT] (dotted lines). However, because of the addition of 
these sources, the final sample of YSO candidates was not uniformly selected, and therefore 
we cannot speak to the completeness of the sample as a whole. It is clear from Figure [11] 
that many YSO candidates in CrA are not selected using solely the c2d/GB method, and 
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we enumerate the reasons why in the following paragraphs. 



The greatest sources of incompleteness in the standard c2d/GB method (§ 13.11) manifest 
in four ways. First, very low luminosity objects may be lost in the group of extragalactic 
sources. An additional bias originates from the Coronet region, with its bright, diffuse 
nebula. Second, in § 13. 21 we discussed the addition of YSO candidates based on their 2MASS 
+ MIPS 24 fim photometry because they were not covered in IRAC (or else are only covered 
in IRAC bands 1/3 and not in 2/4, and vice versa). These sources (see Tabled]) were not 
selected by the standard technique because the method requires detections in all four IRAC 
bands. Third, as discussed briefly in § 13.41 the brightest sources (e.g. RCrA), which saturate 
the camera arrays, were not selected as YSO candidates by the c2d/GB method due to poor 
photometry and subsequent flux measurements. To remedy such incompleteness, we added 
known YSOs and YSO candidates from the literature (see § 13. 41 and Table[7]). Finally, in order 
to select a sample of YSO candidates with a minimum of contamination from extragalactic 
sources and background stars, the c2d/GB method required at least a small infrared excess 
component. Therefore, as mentioned in § I3.3[ our method fails to detect PMS objects that do 
not have infrared excess, but still have some accretion activity, indicated by X-ray emission. 
Ha emission, etc. Consequently, the sample of Class III candidates in Table S] is heavily 
biased, missing many Class III candidates. To correct this bias, we used the ROSAT data 
to include additional Class III candidates (see Table |6]). 

As discussed in § 13. 3^ there are 34 ROSAT detections within our Spitzer field for CrA, 
of which 31 have Spitzer counterparts. To assess the completeness of the ROSAT All-Sky 
Survey at the distance of CrA, we determine the limiting unabsorbed X-ray luminosities for 
different levels of foreground extinction. Sources listed in the Faint Source Catalog have 
a minimum of six counts, corresponding to a limiting count rate of about 0.03 s~^ for a 
typical exposure time of 218 s. The typical exposure time is the median exposure time of 
Faint Source Catalog sources within 3° of RCrA. To find the limiting unabsorbed X-ray 
fiuxes and luminosities, we simulate X-ray sou rces with ass umed spectra using the Portable, 
Interactive Multi-Mission Simulator (PIMMS; lMukailll993l ). The spectra are assumed to be 
absorbe d APEC thermal plasma spectra. The visual extinction of newly identified T Tauri 
stars in iNeuhauser et al.l (120001 ) is Ay < 1 mag, which is the value that we are assuming for 



this estimate, after conversion into absorbing hy drogen column densities using t he empirical 
relation NH(cm"^) ^2 x 10^^ x Ay (mag) (e.g.. lRvteiill996l : IVuong et al.l 120031). Cons erva- 
tively assuming a plasma temperature of 1 keV (~ 10 MK; e.g., iFeigelson et al.ll2005l ). the 
limiting luminosity is log (Lx) = 30.2 erg s^^ for Ay = 1 mag, assuming a distance of 130 pc. 

Subsequently, we use the cumulative X-ray luminosity di stribution functions derived for 
subsets of the PMS population of the Orion Nebula Cluster ( jPreibisch fc Feigelsonl 120051 ) to 
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estimate the completeness of the ROSAT data for G, K, and M PMS stars. The hmiting 
unabsorbed X-ray luminosity determined above corresponds to a cor npleteness of 60-70% 
for G and K stars and 20% for M stars. For the deep Chandra data, iForbrich &: Preibisch 
(120071 ) quote a considerably better limiting luminosity of 4.3x10^^ erg s~^ for no intervening 
extinction which would be complete for all G, K, and M PMS stars. Assuming a foreground 
extinction of Ay= 10 mag yields a limiting luminosity of 6.1x10^^ erg s~^ for 5-count de- 
tections, still complete for all G, K, and M PMS stars. For our purposes, the Chandra data 
are complete for the PMS population of the inner Coronet, while the much less sensitive 
ROSAT data still yield a good estimate of the outlying YSO population. 



4. Distribution of YSOs 



4.1. Extinction Map and YSO Distribution 



With a census of the YSOs in CrA, we can examine the spatial distribution of the 
sources with respect to each other and the dense gas. Usually such a distribution reflects 
the underlying structure of the cloud, providing us with clues about the physical processes 
responsible for its creation. In particular, it can be of help in discriminating between two 
well-known basic structures: hierar chical-type clusters and centrally-condensed embedded 
clusters (see e.g. jLada fc Ladall2003l . and references therein). In the first case, clouds exhibit 
surface density distributions with multiple peaks and more complex shapes, whereas the 
centrally-condensed clusters have a highly concentrated surface density distribution, with a 
relatively smooth radial profile [p^, ~ r~'*). Moreover, the latter usually correlates with mass 



segregation, with massive stars found near the cluster center (see e.g., iLada fc Ladall2003l . 
and references therein). 

In Figure [TU| all 116 YSO candidates are overlaid on the 24 /zm mosaic, along with con- 
tours showing the extinction map we have created as part of the data pipeline using 2MASS 
and Spitzer data. To make this map, the line-of-sight visual extinction (Ay) is estimated for 
each source classified in the catalog as a "star" (e.g., sources classified as YSO candidates 



are ex cluded) based on its SED from 1.25 to 24 fim, adopting the IWeingartner fc Draine 
(I2OOII ) Ry = 5.5 extinction law, which has been shown to be consistent with data from the 



c2d/GB studies (IChapman et al.ll2009l : IChapman fc Mundyl 120091 ). These line-of-sight ex- 
tinctions are then convolved with uniform ly spaced Gaussian beams (with FWHM = 180") 
to construct the extinction map (refer to lEvans et al.l 120071 . for details). The mean Ay in 
our CrA map is ~5 mag, and the maximum Ay is ~30 mag, which is found in the Coronet 
region. The Ay values in our map have not been corrected for any extinction offsets because 
a nearby off-cloud field was not imaged for comparison, and may be systematically high by 



-17- 



as much as 1-2 mag (see Table 27 in lEvans et al.l 120071 ) . For comparison, the extinction 
map presented by iGutermuth et ahl (120091 . and shown in our Figure [1]) , which is based only 
on 2MASS observations, exhibits a mean Ay of ~0.3 mag and maximum Ay of ~20 mag 
over the same region covered by our Spitzer observations. The lower values can perhap s be 
explained by the use of a different extinction law, Ry = 3.1 ( jRieke fc LebofskyI 119851 ). for 
the near-infrared observations. Correcting for the different extinction laws, the near-infrared 
extinction map would have Ay values that are ~ 25% higher. 

As can be seen in Figure [TIJ| the morphology of the entire region is mainly a clustered one, 
extending toward the southeast, although the stars are mainly clustered near the Coronet, 
which indeed harbors a Herbig Ae star (RCrA) in its very center. It is striking how the 
Class I sources (red -|-) are found in the highest extinction regions, highly clustered in the 
Coronet, the Class II sources (green diamonds) are spread a little further out, and the 
Class III sources (blue squares), spread even further out into the molecular cloud, forming 
a sort of halo around the Class I and Class II sources. A distribution where the Class I 
sources are more densely clustered than the Class II sources is consistent with observations 



of other nearby star-forming regions such as Lupus, Serpens, and NGC 1333 (e.g. iMerfn et al. 



2008; Gutermuth et al. 2008, 2009; Winston et al. 2009 



Bressert et al.ll2010l ). A quantitative 



analysis of the clustering in CrA will be discussed in the next section, § 14.21 



4.2. Clustering Analysis 



Using all 116 YSO candidates selected and/or compiled in this paper (see Figure [TOj) 
and listed in Tables jH El El and [71 w e perfo rm the clustering analysis identical in method 
to that pres ented in IGutermuth et al.l ( l2009l ) for comparison. The sample sizes are different 



because IGutermuth et al. 



fl2009h used only the Spitzer data of CrA available from the GTO 
PID 6 program (see Table [1]), which covers a much smaller field of v iew than the region 
we include in our study. For this smaller field, IGutermuth et al.l ( 120091 ) classified a total of 
35 YSO candidates: 25 Class II and 10 Class I, for a Class Il/Class I ratio of 2.5. They 
found one core of more than 10 stars in CrA, consisting of 24 members. For the core of 24 
members, there is a Class Il/Class I ratio of 2.0. 



In their analysis, which we also employ, IGutermuth et al.l ( l2009l ) isolated the dense 
"cores" of clustered star-forming regions via analysis of the Minimal Spanning Trees (MST) 
of the S'j>zteer-identified YSO positions in each region. The MST is defined as the network 
of lines, or branches, that connect a set of points together such that the total length of the 
branches is minimized, containing no closed loops. Once the MST is constructed for all YSO 
positions, a critical branch length, Lcru, is adopted whereby all branches in the network that 
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are longer than Lcrit are excluded. In the case of the iGutermuth et al.l (120091 ) analysis, Lent 
is computed based on a two-line fit to the cumulative distribution function of the lengths of 
the MST branches. After the branches longer than Lcru are excluded, smaller sub-networks 
generally remain connected. Any of these sub- n etworks that are c omposed of 10 or more 
members (a common, but arbitrary choice, e.g. lAUen et al.l 120071 : IGutermuth et al.l |2009|) 
are defined as "cluster cores." The size of a core is quantified in two ways, by its circular 
radial size, Rcirc, and by its effective radial size, Rhuu- The circular radial size, Rcirc, is 
defined as the maximal radial distance to a group member from the cluster center point, 
where the cluster center point is the geometric center of the convex hull polygon enclosing 
each gro uping (the geometric cen ter is the position of the center of mass of the convex hull 
polygon; IGutermuth et al.ll2009[ ). The effective radial size, Rhuii, is defined as the square 
root of the area (divided by vr) of the convex hull polygon that contains each grouping. The 
ratio Rcirc/^huu is the aspect ratio of elongated distributions, i.e. it is used quantify the 
degree to which a core is elongated. 

When we analyze our list of all 116 YSO candidates (which includes Class I, Flat 
spectrum. Class II, and Class III candidates), we find one significantly sized core consisting 
of 68 members, which is located in the Coronet region. The ratio of Class II to Class I 
candidates of the entire sample is: 2.26 (43 Class II, 19 Class I and Flat spectrum sources); 
for the core of 68 members. Class Il/Class I = 1.77. If we include cores that consist of 
fewer than 10 members, the next two largest groups consist of 6 members (the region near 
IRAS 32) and 5 members (a region to the northwest of TYCrA). The CrA cluster is clearly 
dominated by the cluster core around the Coronet. 



Figure IT2] shows the residual MST sub- networks for the YSO positions after those longer 
than Lcrit have been excluded. Black points represent each source, and the branches (grey 
lines) connect those sources that are more closely spaced than the critical length measured 
in the cumulative distribution of MST branch l e ngths , which is shown in Figure [T31 The 
critical branch length found in IGutermuth et al.l ( 120091 ) for CrA, corrected for a distance of 
130 pc (a distance of 170 pc was used in their paper) was 0.116 pc (see their Figure 1 for 
CrA). Using all 116 YSO candidates, we obtained a Lcrit = 0.164 pc (see the top panel of 
Figure US]) • 

We performed the same analysis for the 62 YSO candidates that had either Class I, 
Flat spectrum, or Class II SEDs (i.e. the Class III sources were excluded). Again, only 
one significant core was found, consisting of 34 members. The bottom panel of Figure IT2] 
shows the MST of the YSO positions for these sources, and in this case, we obtained a Lcrit 
= 0.152 pc (see the bottom panel of Fig ure [T3|). Table |8| sumrn arizes the results of this 



comparison, as well as the results of the IGutermuth et al 



(|2009[ ) study, including values 
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of Rcirc, Rhuii, the aspect ratio Rcirc/R-Lii' the density for each core. For CrA, when 
looking at the core of 68 members, the value for Rcirc = 0.59 pc, with an aspect ratio of 
2.36, indicating that the cluster is elongated. To illustrate this elongation, the convex hull 
for each grouping is shown (solid black line) in Figure [T2i 

The addition of new mer nbers from this paper cha nges the membership of the main CrA 
core from 24 members in the iGutermuth et al.l (120091 ) study, to 68 members in our study of 
all classes of YSOs. However, the mean densities are quite similar since our analysis covers 
a larger area and includes known YSOs from the literature (see Table [H]). In addition, the 
elongation of the core changes only slightly, from an aspect ratio of 2.51 to 2.36. A careful 
look at the differences in the top and bottom panels of Figure [12] shows that the Class III 
population extends further to the west, indicating an evolutionary stage gradient to the west 
in CrA. This confirms our qualitative analysis of § 14. 1^ obtained from looking at the spatial 
distribution of YSOs (see also Figur e [TOj). and the noted elongation of the CrA cloud fr om 
previous studies, as discussed in § [1] (INeuhauser et al.ll2000t iMamajek fc Feigelsonll200ll ). 



5. Submillimeter Analysis 



Six YSOs were detected in the dust continuum at 225 GHz by the SMA, and are listed 
with their flux densities in Table [HI and shown in Figures and [TH The YSOs IRS 5, IRS 9, 
CX0 34, and CrA-24 were not detected, above a 3a rms level of 11 mJy (IRS 5) or 22 mJy 
(IRS 9, CX0 34, and CrA-24). The detected emission in all cases is compact and centrally 
concentrated, and is thus assumed to be primarily due to compact disk emission. Point 
source and 2-D Gaussian models were fitted to the detected YSOs, and for all sources but 
IRAS 32, the results were consistent with a point source, indicating unresolved emission on 
size scales of about 400 AU, consistent with disk sizes measured for o ther embedded YSOs 
( J0rgensen et al.l 120071 12009| ). and for T Tauri stars in nearby regions (lAndrews fc Williams 
20071 ). For IRAS 32, an extended component is also present (more flux is recovered from 
fitting a Gaussian than from a point source alone), and is assumed to b e due to inner envelop e 
emission; this extended component was also seen in infrared imaging by lWilking et al.l ( 1l992l ). 
To remove this component, IRAS 32 was re-imaged using only baselines > 20 kA, with the 
result being consistent with a point source. More detailed model fitting of IRAS 32, including 
2D radiative transfer, will be presented in a later paper. The flux densities reported in 
Ta bleware cons i stent with the values for SMA 2 (= RS9) and SMA 1 (= IRS 7B) reported 
by lGroppi et al.l (120071 ). who observed with the SMA at a slightly higher frequency of ~271 
GHz. 



If the disks are optically thin, the continuum fluxes provide a direct measure of the 
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masses. For marginally opticall y thick emission, possible for such disks (|j0rgensen et al. 



20071 : [Andrews fc Williamsl 120071 ). the mass becomes a lower limit to th e true mass. The cor- 



rections are usually small, however. Previous studies of T Tauri stars ( lAndrews fc Williams 



20071 : iLommen et al.l 120071 ) indicate that the ratio of optically thick-to-thin emission is typi- 



cally of order 0.2-0.3. This translates into an underestimation of the mass by up to a factor 
of two. Given the other assumptions (constant disk temperature, no radiative transfer mod- 
eling), this is a small c orrection. To compute disk masses, we use equations (3) and (4) from 



J0rgensen et al.l (120071 ) with their assumptions for the dust temperature (30 K) and 1.3 mm 



opacity (assuming a gas-to-dust ratio of 100). For the assumed distance of 130 pc for CrA, 
the disk mass at ~ 225 GHz (1.3 mm) is 



M = 0.236 M 



© 



F 



(2) 



where F is the flux density measured in Jy. The results for the detected sources are listed 
in Table [91 The uncertainties listed in that table are those returned by the fitting procedure 
in MIRIAD, and are similar to the 3a rms values reported earlier. So although the flux 
density of RCrA reported in this table is only 1.5 times the reported uncertainty, it is a well 
detected source in the maps (about 5cr). 

For the non-detected sources, the upper limits to the disk mass are 0.003 Mq for IRS 5 
and 0.005 Mq for IRS 9, CrA-24, and CX0 34, that is, only a few Jupiter masses (Mj), 
and signific antly less than the rn inimum mass solar nebula (MMSN), of 10 - 70 Mj (~0.01 
- 0.07 Mq; IWeidenschillingI 119771 ). In these cases, the masses are either very low, the dust 
has been processed into bodi es significantly l arger than a millimeter, and are thus mostly 
invisible to our observation s (ILommen et al.l 120071) . or the dust is absent, perhaps due to 
clearing in a close binary ( 1 Jensen et al.l 119961 ) . It is known that IRS 5 is a close binary 
system (~100 AU ), and so this last sceriario is a real possibil i ty (see further d iscussion of 
IRS 5 in § IA.1.191 Ichen fc Grahamlll993l : iForbrich et allbood : Ichoi et al.l[2008[ ). The disk 
masses reported for the detected sources, aside from RCrA, are consistent with the MMSN. 

The masses reported here are similar to masses of Class and Class I protostars 
derived in an identic al manner (0.017 - 0.089 Mq, and 0.009 - 0.076 Mq, respectively; 
J0rgensen et al.l 120091 ) . For Class II YSOs (T Tauri stars) also observed with the SMA, 
but anal yzed with more detailed mo dels, the masses are similar or slightly higher (0.003 - 
0.18 Mq; [Andrews fc WiUiamsl l2007l). The fluxes reported for the Class II YSOs are similar 
to those reported here (lAndrews &: Williamsll2007l ) . and lie at a similar distance (Taurus and 
Ophiuchus), so calculating masses using equation ([2]) would produce masses similar to our 
results. 
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Observations at centimeter wavelengths show that the emission we have observed is 
most hkely due to thermal dust emis sion, with essent i ally no contaminatioii from free-free 
emission, except perhaps for RCrA (jChoi et al.ll2008t iMiettinen et al.l 120081 ) . All sources, 
other than RCrA, are either weakly detected at cm wavelengths (< 1 mJy), or not at all, 
with spectral indices between 4.9 and 8.5 GHz that are fiat or negative. Although RCrA is 
also only weakly detected at cm wavelengths, it has a positive spectral index in the range 
0.4-0.9, which implies that up to 5 mJy of the emission observed at 1.3 mm could be due 
to free-free emission, and so the disk mass should be sli ghtly reduced to 0.0 08 Mq, which is 
consistent with the upper limit of 0.012 derived by iGroppi et al.l ( 120071 ). 



6. Outflows and Jets 



6.1. Outflows and Jets Detected in H2 and Spitzer maps 



Our H2 survey maps an area of about 20'xl7', which spans from ~5' east of the Coronet 
region to ~4' west of S CrA, and from ~3' north of TYCrA to ~3' south of HH 101 (see 
Figure H]). The mapped region is about one tenth of the Spitzer area, and is shown in 
Figure [TBI as a green polygon on an IKAC/ Spitzer three-color map of CrA. 

Our H2 continuum-subtracted mosaic reveals nearly 100 knots, some of them coincident 
with, or close to previously known HH objects. These knots are shocked regions along the 
outflows, usually traced by H^ and/or Fe II emission lines in the near- infrared (see e.g., 
Caratti o Garatti et al.ll2006[ ). and ionic species in the optical (as O I, S II, and Ha), also 
known as HH objects, which t race dissociative sh ocked regions with lower extinction and 
higher temperatures (see e.g., lEisloffel et al.lll994l ). Once the continuum is removed from 
the H2 narrow-band images, what is left is the emission line alone, and thus knots and jets 
can be unambiguously identified. We detect 8 H2 counterparts out of 14 HH objects located 
in the mosaic FoV, while no evidence of molecular emission is observed from the remainin g 
6, namely HH96, HH729, HH730, HH734, HH736, HH860. Following [Pavis et al.l feoioh . 
we name the detected H2 flows MHO 2000-2014, followed by letters indicating knots and 
sub-structures inside each flow. It is worth noting that such a grouping is based on the 
flow morphology a lone, and does not necessarily imply a real association to a single outflow. 
Davis et al.l (120 lOl ) designated five MHOs (2000-2004) in the Coronet region, on the basis of 
morphology alone, because no PM analysis was available. Usually such a grouping is effective 
for isolated outflows and sources, where the outflows do not overlap and it is straightforward 
to recognize the outflow structure (chains of knots more or less aligned, bow-shocks indicating 
the outflow direction, etc.), which also points to the exciting source. On the other hand, 
it is often particularly difficult to disentangle different flows in crowded regions such as the 
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Coronet, where several flows overlap, and the exciting sources are enclosed in a relatively 
small area. Here, we also use our P.M. analysis (see § 16.21 and Appendix |B]) to recognize 
and group each flow. However, our analysis and thus our groupings are not resolute for 
several cases (see e.g., § 16.31 and Appendix [B|) due to the crowding i n the Coronet and 
the uncertainties in our PM measurements. As a result, we keep the iDavis et al.l (120101 ) 
nomenclature, adding ten more newly detected MHOs, i. e. MHO2005-2014, and include 
letters to pinpoint recognized knots and sub-structures within each MHO. We stress that 
such a grouping does not necessarily imply a real association to a single outflow. 

The majority of the flows originate from the Coronet region and are radially distributed 
around it, with a larger concentration in the east and north-east regions (namely MHO 2000- 
2002, 2009, 2011-2014), and west and south-west regions (namely MHO 2003-2008), with 
MHO 2005 positioned around S CrA. Finally, MHO 2010 is positioned around IRAS 18595- 
3712 (IRAS 32). We identify about 14 Hs flows in CrA. 

Most of the known HH objects, and all the detected H2 objects, are visible in our 
IRAC maps. Indeed, the IRAC maps con tain molecular hydrogen lines that may be shock- 
excited in protostellar outflows (see e.g. iDe Buizer fc Vaccal I2OIOI . and references within). 
Generally, the H2 lines in band 2 (4.5 fim) and, partially, in band 1 (3.6 /im) have excita- 
tion con^itioiis_guite_s^^ tracing shocks at T~2000K (see 
e. g. ICaratti o Garatti et al.ll2006l : iTakami et al.ll2010[ ). In particular, IRA C band 2 contains 



bright molecular hydrogen lines and can be used to detect knots (see e. g. ICveanowski et al. 



200al2009l : lYbarra fc Ladall2009r ). often called "extended green objects" (EGOs, ICvganowski et al. 
20081 ) . for the common coding of the 4.5 /im band as green in the IRAC three-color images. 
To better identify these features in the IRAC/Spitzer mosaics, we thus constructed IRAC 
three-color images, using 3.6 /im, 4.5 /im, and 8 /im (blue, green, and red, respectively, see 
Figure [TBI) . and identifled the new knots/EGOs by means of their colors and morphologies. 
Usually the outflow morphologies in band 2 and, partially, in band 1 resemble those observed 
in our H2 images, mostly delineating relatively high-excitation regions as bow shock heads 
or clumpy structures/knots along the flow axes. On the other hand, band 4 and 3, which 
contain relatively low-excitation H2 lines, often delineate wing s and wakes of bow shocks, 
and in general, milder shock regions (for a detailed analysis see lTakami et al.ll2010l ). 



In addition to the H2 detected features, our Spitzer maps reveal three more flows posi- 
tioned outside the H2 mosaic, and named Spitzer outflow 1, 2, and IRAS 32 outflow, which 
include 30 additional new knots. Because we lack spectral conflrmation, these flows could 
not be named MHO. 



H2 and Spitzer newly detected knots and all previously known HH objects (including 
those not detected in H2) are labeled in our H2 map (see Figure H]) and in our three-color 
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IRAC mosaic (Figure [16]) as blue, red, and green crosses, respectively. Additionally, a list of 
known HH objects of the CrA star-forming region is reported in Table [10], and a complete 
list of knots detected in our Spitzer and H2 images is given in Table [TT] 

In Appendix |B] Sections [B.ltlB.2t and lB.3t a description of individual MHOs and Spitzer 
flows is given. 



6.2. Proper motions of H2 knots 

To derive proper motions from our multi-epoch set of images, we followed the technique 



described in detail by lCaratti o Garatti et al.l (120091 ). As shown in Figure H] the latest epoch 



map {SofI 2007) encloses all the remaining four, thus we identify the SofI 2007 map as the 
common reference frame for our PM analysis. All earlier epoch maps were scaled and aligned 
to the reference frame with sub-pixel accuracy. The match was done using several field stars 
observed in each map (see § 12. 3p . and optical distortions were corrected by means of the 
geomap (with a polynomial fit of 3rd order) and geotran routines in IRAF. The resulting 
errors are represented by the residuals of the fits, and range between ~0.2 and ~0.4 pixels 
(i. e. ~0'.'06-0'.'12), depending on the seeing and the number of stars in each map. 

Knot shifts were determined between image pairs (i. e. last epoch and previous epoch 
maps) using a cross-correlation method. Knots were identified in the continuum-subtracted 
maps, and a rectangle was defined around each knot, enclosing its 3 a" contour. The first 
epoch image was then shifted with respect to the last epoch image, with a sub-pixel accuracy 
of 0.1 pix and then multiplied by it. For each shift (x,y), the total integrated flux (f) in the 
rectangle around each knot was measured in the product image. For the final shift for each 
knot, we used the position of the maximum of f(x,y), determined via a Gaussian fit. 

Systematic errors for shift measurements were obtained by modifying the size and shape 
of the rectangle that defines each object. The resulting range of values gives a systematic 
error, which depends on the S/N, shape, and fiux variability of the knot. This error is usually 
comparable to the alignment error. Thus, each single shift error was derived combining the 
alignment error of the two epoch images and the uncertainties in the centering routine. 
Usually, for each pair, the final error was around 20% of the P.M. value. It is worth noting 
that the number of P.M. measurements differs from knot to knot, depending on the number 
of matching maps. Thus knots located in the Coronet, with three or four P.M. measurements, 
usually have smaller errors. 

In Figure [TT] we show the fiow charts of the CrA star-forming region. To display MHO 
motions for the entire region, proper motions, indicated by arrows, are reported in units of 
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500 yrs. The results of our P.M. measurements are given in Table [T2| which lists the knot 
ID, the measured P.M. in arcsecyr"^, the derived tangential velocity (ftan) in kms~^ at a 
distance of 130 pc, and the position angle (P.A.) in degrees of the P.M. vector. P.M.s were 
measured for knots observed in at least two different epochs. Moreover Table [T^ reports only 
P.M.s with S/N ratios larger than 2a. As a result, 44 measurements are reported, namely 
those from MHO 2000, 2002, 2003, 2007, 2013, the majority of knots in MHO 2001, 2009, 
2011, 2012, 2014, and some knots in MHO 2004, 2006, and 2008. 

The derived P.M. values range between 0.024 and 0.256 arcsecyr^^, corresponding to a 
ftan between 15 and 158kms~^ (at a distance of 130 pc). The uncertainties vary from 0.006 
to 0.045 arcsecyr"^ (i.e. ~4 to ~28kms~^), depending on the knot S/N ratio, the number 
of epochs, and the accuracy in the alignment. Such uncertainties, in turn, reflect on P.A. 
errors, which span from 2° to 27°. 

A quick inspection of Figure [T71 confirms that all the studied knots, except MHO 2008 D, 
are launched by YSOs inside or close to the Coronet. Several flows originate in the eastern 



Nutter et al. 


2005; 


Groppi et al. 


2007) 



A detailed analysis of P.M.s in single outflows is reported in Appendix |Bl There we 
show in more detail the flow charts of each region, where P.M.s (in 100 yr) and their error 
bars are indicated by arrows and ellipses, respectively (Figures | 



Finally, detailed discussion on outflow occurrence in CrA, and possible driving sources 
can be found in § 16.31 



6.3. Outflows and possible driving sources 



Our H2 survey along with the Spitzer IRAC maps and proper motion analysis shows at 
least 14 different H2 flows in CrA. This number increases to 17 if we also include HH objects 
that do not have H2 counterparts but likely originate from one of the detected YSOs. How- 
ever, the number of flows could be larger than this because we are not considering some of the 
HH objects which do not have a viable candidate driving source. Outflows and their possible 
driving sources were selected on the basis of the outflow morphology and proper motion anal- 
ysis (see § 16. II and Appendix [B|) . Candidate driving sources, their possible flows and position 
angles, are reported in Table [131 Also blue- and red-shifted lobes (hereafter blue and red 
lobes, respectively) have been tentatively indicated considering ra dial velocity measurements 



from this work (see ^ [5] and Figure [T5|) and frora the literature (jHartigari fc Graham 



Levreault 1988 : Anderson et al. 1997[ " Davis et al. 1999; Groppi et al 
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van Kempen et al.ll2009l ). Several question marks indicate the uncertainty in matching flows 
and sources. A clear match was obtained for those flows close to the driving source (IRS 1, 
IRS 2, IRS 6, SMM2, SMA2, IRS7B, CrA-41), or driven by isolated YSOs (IRAS 32). 

The uncertainties in Table [T3] are mainly due to the particular morphology of the Coronet 
region, as well as the uncertainties in the P.M. measurements. The majority of the flows 
originate from sources inside, or close to, the Coronet, in particular from the regions around 
SMA2 and IRS 7A, where four known YSOs are located (namely SMA2, IRS 7A, B and 
CX0 34). However, we note that at least six different overlapping flows are launched from 
this region of the Coronet, likely ind icating the presence of other embedded YSOs (see e. g., 
Groppi et al.l 120071 : 1 Choi et al.l l2008l . see also Figures [18] and |22]) . The overlap of the flows 
and the errors on the P.M.s do not allow us to pinpoint the driving sources, with the possible 
exceptions of SMA 2, IRS 7A, and IRS 7B, where the detected jets appear to be launched from 
or close to the positions of the sources (see Figures [221 and [23l and Appendix lB.l|) . The SMA 2 
flow should include MHO2001J, A, and MHO 2009 to the east, and MHO 2001 Q to the 
west. The IRS 7B flow should include MHO 2002 C (see Figure [22]), MHO 2012, and possibly 
HH736 in one lobe (see Figure [T5]) . and, possibly, HH731B in the other (see Figure [57] and 
Table [12]). 

There is also strong evidence of a very extended (~lpc) east-west wide precessing 
flow (~15°), with blue and red lobes positioned east and west of the Coronet, r espec- 
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whereas MHO 2004 D, Spitzer outflows 1 and 2, and MHO 2005 (see Figures [391 [401 and[4T]) 
match the more extended CO red lobe. The driving source is likely to be one of the youngest 
sources in the region around SMA 2 and IRS 7A. SMA 2 would be the obvious candidate, 
but it already drives another jet, which is also east-west oriented, but with smaller size and 
precession angle (~4°) (MHO 2002 and MHO 2009). Ho wever we n o te th at SMA 2 could 



harbor a binary or multiple system, as e.g. reported by IChoi et al.l ( |2008[ ). In this work, 
they detect three different sources, namely B9a/b in the centimeter continuum and CT 2 in 
the 7 mm continuum, positioned a few arcseconds north of IRS 7A, and at least two outflows 
(FPM 10 and FPM 13). We observe at least three different precessing flows escaping from 
the SMA 2 - IRS 7A region (namely knots MHO 2001 M, N, and O). Another possibility is 
that this large outflow is driven by IRS 7A. 



Another large CO flow extends roug hly SSW to NNE (iHartigan fc GrahamlEgSTl iLevreaultl 



19881 : [Anderson et al.lll997l : iKned 120051 ). with HH 101 (MHO 2007) coincident with the blue 
lobe. HH99 (MHO 2000) was previously believed to be the count erflow, but both geometric 
and kinematic considerations ruled out this hypothesis (see e. g., [Hartigan &: Graham[[l987[ : 
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Davis et al.lll999l ). More recent CO observations ( lKneell2005[ ) place the peak of the red lobe 
far to the north-west, i.e. close to the position of the newly discovered MHO 2014, which is 
likely the counterfiow, as both morphology and P.M. analysis also suggest. IRS 5 and 5N are 
positioned between the two lobes, thus one of the two is likely the driving source. Based on 
our SMA results, IRS 5N seems to be the most likely candidate, since no continuum emission 
from IRS 5 is detected (see also SI51 IA.1.191 and IA.1.20p . 



There is no clear evidence of the HH 99 B (MHO 2000 A) counter-flow in our analysis. 
The RM. analysis indicates that MHO 2000 A and B (HH99B and A) are likely part of two 
distinct flows. MHO 2000 A could be driven by CrA-24, one of the newly discovered Class I 
YSOs, or by IRS 9 (see also Appendix [B] and Figure [T8l) . MHO 2000 B seems to be part o f 
a larger flow driven by IRS 6 (see also. Appendix Rl [B| Figure \TE\ and IWang et al.l |2004| ) . 
Another Class I YSO, namely CrA-43 (or SMM2) in the Coronet outskirts (see e.g. Figures 
H] and [22]), is also driving an H2 outflow, SW-NE oriented (P.A. ~125°), which includes 



MHO 2002 B, G, and possibly H, and HH 734 (see also. Appendix |Al El and IWang et al. 
20o3). 



On the basis of our morphological and kinematical analysis (see Appendix[B]), in Table [T^ 
we tentatively associate the remaining observed H2 and HH objects with nearby detected 
YSOs. These possible matches are labeled in the table with a question mark. It is worth 
noting that we clearly detect only one lobe in several flows. I ndeed this has alre ady been 
observed in relatively crowded and confused regions (see e.g., iDavis et al.l l2009l ). but our 
sample has an extremely high rate of single- lobe flows (9 out of 17), which is particularly 
disturbing. Although this can be ascribed both to the morphology of the region, the high 
extinction (see Section HTTj) . and the limited sensitivity of our survey, it could also indicate 
mismatches between flows and sources. 



6.4. Incidence of outflows vs number & class of YSOs 



A more detailed analy sis of Table 



early stage YSOs (see also, iDavis et al. 



T5] co nfirms that H2 outflows are mostly driven by 



20091 ). Indeed, Class source SMA 2 drives an H2 



flow, as well as 10 out of the 14 YSO candidates classified as Class I. All of the Class I 
YSOs located near or in the Coronet seem to drive outflows. Three out of 43 Class II YSOs 
have H2 flows, which becomes 6 out of 43 Class II YSOs if we also consider flows with HH 
objects. None of the Class HI sources shows any sign of outflow activity, indicating that 
accretion activity has reached completion in this stage, as one would expect. On the other 
hand, the pe rcentage of H2 and atomic jets detected in Class II sources by this, and previous 
surveys (e. g. IWang et al.ll2004j ) seem to be too small (~14%) for this class of objects, because 
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accretion should still be ongoing in Class II sources, and thus outflows should be expected. 
This would then suggest that these surveys are not sensitive enough to detect such outflows, 
or that there are mechanisms, for example the cloud magnetic field, that inhibit the outflow 
activity in more evolved YSOs. Nevertheless our data confirm that the accretion/ejection 
activity decreases with protostellar age. The H2 outflow sources are protostars (Class O/I 
sources) rather than disk-excess sources (mostly Class II T Tauri stars), though a few of the 
latter source sample do seem to power H2 flows, and HH objects. However, if we exclude 
the IRS 6 flow, the H2 emission in these Class II flows is extremely reduced, whereas the HH 
and ionic emission are still well detected. 



7. Overall Cloud Properties 



7.1. Disk Fraction 



We have compiled a fairly complete estimate of the Class HI population in CrA (see 
§ 13. 5p . and as a result, can examine the disk fraction of YSO candidates. There are 90 
YSO candidates of all classes within the Ay = 2 mag contours; of those, 42 are Class HI 
YSO candidates. Theref ore the disk fractio n in C rA is 53%. Recent disk fraction estimates 
for CrA include 53% by iLopez Marti et al. who also included a population of disk- 



less sources, and 67% by iMeyer &: Wilkingl (j2009i) . who obtained the disk fraction from an 
extinction-limited sample of sources and used spectroscopic information to estimate an age 
range of 0.3—3 Myr from their H-R diagrams. D isk fraction values simi lar to ours have 
been i nferred for young clusters NGC 7129 (54%; ICutermuth et al.l booj ) an d NGC 2264 
(52%:|Haisch et al.ll200ll ). Based on comparison with theoretical PMS tracks, iHaisch et al. 
( I2OOII ) estimate an age of ~2.5 Myr for NGC 2264, whic h is consistent with the age range 
of 0.3-3 Myr estimated for CrA jMever fc Wilkingjl2009h . 



7.2. Star Formation Rate and Efficiency 



The mass of the CrA cloud can be estimated, using the Spitzer + 2MASS extinction 
map discussed in § 14.11 To comp are our results to those for the c2d and GB clouds, we 
employ the same method used by iHeiderman et al.l (120101 ). They use the lowest resolution 
extinction map (FWHM = 270") to calculate Mgas,cioud by summing up all the pixels (SAy) 



above an A y = 2 mag th resho 
cm~^mag~^ (jDrainell2003[ ) f or a 



d, convert to column densit y using N/^/Ay = 1.37x10^^ 



Weingartner Sz Draind ( 120011 ) Ry = 5.5 extinction law, and 



obtain Mgas,cioud = 279 Mq ( IHeiderman et al.ll2010l ). By following this procedure we obtain 
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the same value for the cloud mass and the area of the cloud within the Ay = 2 mag contour 
(i.e. 0.6 deg^), which translates to Adoud = 3.03 pc^ at a distance of 130 pc. By dividing 
Mgas.cioud by the area of the cloud, we obtain the average gas surface density of the cloud, 

The range of gas surface densities for the c2d and GB clouds 

45-140 M0PC-2, and CrA falls in the 



Sgas,cloud = 92.1 M(T,PC 



discussed in iHeiderman et al.l ( I2OIOI ) is Sgas,cioud 
middle of that range, similarly to Chamaeleon I, Auriga, and Perseus. 



Because we include additional YSO candidates not included in the Heiderman et al. 



(I2OIOI ) study (i.e. Nysctot = 48, versus 41), the resulting star formation rate (SFR) for CrA 
is a factor of ~1.2 higher than theirs. Note that only Class I, Flat spectrum, and Class II 
sources located in the area of the extinction map above t he Ay = 2 mag contours are con- 
sidered. Following the analysis of IHeiderman et al.l ( 120101 ). we assume a mean YSO mass of 
0.5±0.1 Mq over a peri od of star formatio n of 2±1 Myr, based on e stimates of the time spa n 
of the Class II phase (lEvans et al.l 120091 ) . Using equation (8) in IHeiderman et al.l ( I2OIOI ). 
we compute the SFR surface density to be Ssfr = 3.96 M0yr~^kpc~^, which translates to 
a SFR of 12.0 M0Myr~^. CrA has a relatively low SF R, compared with Perseus, Serpens 
and O phiuchus (96, 56, and 73 MQMyr"^, respectively; lEvans et al.ll2009l : IHeiderman et al. 



2OIOI ). and is similar to the Lupus clouds, in particular Lupus I II, Lupus V, and Lupus 



VI (17, 11, and 11 MQMyr ^, respectively; IHeiderman et al. 



2010f ). The SFR surface den- 



sity, however, is the highest of all the clouds, wit h only Serpens showing similar values 
(E oFR = 3.29 Mf7^yr~ ^kpc~^; Heiderman et al.l 2010 ). Similar comparisons were also made 



by iLada et al.l (I2OIOI ) , where a useful graphical comparison of CrA with the c2d and GB 
clouds can be found (ILada et al.ll2010l . particularly, their Figures 2 and 4). 



For comparison with the other c2d/GB regions, we also compute the star formation 
efficiency (SFE) for CrA using the following relation: 



SEE 



M, + M{cloud) 



(3) 



The total mass of the stars, is estimated by multiplying the total number of stars (Nyso.tot 
= 48) by an assumed average stellar mass of O.5M0. Using the entire cloud mass computed 
above (Mgas.cw = 279 Mq), we obtain a SFE for CrA of ~0.08. fevans et al.l J2OO9I ) compiled 
SFE for several of the c2d clouds (see their Table 4), which range from 0.03 (Chamaeleon 
I) to 0.063 (Ophiu chus). The overall value of SFE for the GB cloud IC 5146 is ~0.05 
( iHarvey et al.ll2008l ). In general, CrA has the highest SFE of all the other clouds observed 
in c2d/GB surveys. This result is still valid even if we use the smaller, NYso,tot = 41 value 
used in the IHeiderman et al.l (12010! ) study, so the effect is not merely a result of the addition 
of YSO candidates (although this may represent a useful measure of the uncertainties that 
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the st udies for the other c2d/GB regions may face due to incomplete YSO hsts). iKirk et aL 
( 120091 ) found an SFE of 8% for one of the GB regions, the Cepheus Flare region, when 
looking only at cores where groups of YSOs were located (and an SFE of 1% for the quiescent 
protostellar cores). In the case of CrA, perhaps we are seeing a small cloud of gas efficiently 
producing a cluster of stars, and no additional regions of diffuse gas diluting the SFE. 



The high SFE in CrA has been noted in the past (e.g. 45%, 40% and 14%: IWilking et al. 



19861 : iHariu et al. 



1993 



Tachihara et al.l |2002|. respec tively) . The very high estimate by 



both the IWilking et al.l ( 1l986l ) and iHarju et al.l (119931 ) studies is due to the fact that they 



used only the main Coronet core sur rounding RCrA to esti r nate SFE (spec i fically , the area 
corresponding to the 2 ^um survey by lTaylor fc Storeylll984j ). IWilking et al.l ( 119921 ) compare 
the SFE of CrA with that of the p Ophiuchi cloud, which is similar in terms of distance and 
angular extent; also, neither cloud has a star more massive than BOV. They note that CrA 
has fewer YSOs than p Oph, and attribute t hat to the overall lower mass of both the low- 
and high-density molecular gas in the cloud (jPame et al.lll987l : lLorenlll979l ). 



Tachihara et al.l (120021 ) go a step further and compare the SFEs for a large sample of 



179 cores in Taurus, the p Oph cloud, the Ophiuchus north region, the Lupus clouds, L1333, 
Southern Coalsack, the Pipe nebula, and CrA star-forming regions. They find an average 
SFE of roughly ~ 10% for all the cores, with an SFE o f 14% for all of CrA, an d 40% for the 
cluster-forming core of CrA. Using C^^O observations, iTachihara et al.l (120021 ) note that the 
regions with cluster-forming cores (in addition to CrA, this includes p Oph, Lupus III and 
Chamaeleon I) have similar cloud structures, i.e., a head, where the active star formation 
occurs, and a tail that extends from it, suggesting that star formation has been triggered 
by an external shock. In order to explore this observation, they measure the average line 
widths, AV, for the starless cores which have not been disturbed by outfiows from YSOs, 
and therefore represent the initial turbulent velocity in the cloud. The SFE, estimated using 
the total cloud mass from the ^'^CO clouds, is then plotted versus AV (see Figure 11 in 



Tachihara et al.ll2002l ). For CrA, they estimate a cloud mass by applying the ratio of the 
total C^®0 core mass to the ^^CO cloud mass to be 30%, and estimate SFE to be ~4%. 



In general, ITachihara et al.l ( 120021 ) find that the regions actively forming stars, and 
whose ^^CO-based SFE is greater than 2% (which includes CrA), all have AV < 0.7 km s~^ 
(conversely, the less active regions have AV > 0.7 km s~^). Based on these results, they 
suggest that a star-forming region with an initial turbulent velocity that is low will form 
stars spontaneously; if hit with a shock wave, the shock wave will then trigger the formation 
of a cluster in the head, forming a head-tail structure. Our observations and clustering 
analysis of CrA (see § 14. 2p are consistent with this picture: we see an elongated cluster core 
which is dominated by star formation in the "head," the Coronet, with a more extended 
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population of PM S stars, and a quiescent clump in the southeastern "tail" of the cloud 
( iHarju et al.l Il993l ) . As mentioned in § [H it has been shown that Cr A is moving away from 
the Upper Centaurus Lupus (UCL) association with an expansion velocity similar to that 
of the Sco-Cen superb ubbles (IMamajek fc Feigelsonll200ll ). specifically the H I shell of Loop 
I, suggesting that the expanding H I shell collided with the CrA cloud a few million years 
ago, triggering star formation in the Coronet, which is the part of the cloud facing the UCL 
( jHarju et al.l 119931 ) . The "tail" of CrA is located downstream from the UCL association, 
which is also consistent with the triggering scenario. 



8. Summary 

In this paper, we present the Spitzer Space Telescope, Submillimeter Array, and H2 
observations, combined with ROSAT and 2MASS, of the Corona Australis YSO population, 
its distribution, and outflows. 

• A total of 116 YSO candidates are identified, where 14 are classified as Class I, 5 are 
Flat spectrum, 43 are Class II, and 54 are Class III. These candidates were selected 
based on Spitzer, 2MASS, ROSAT, and Chandra observations, as well as an extensive 
search of the literature. Of the 116 YSO candidates, there are 12 which have not been 
selected as candidates by any previous study, where 3 are classified as Class I/Flat, 8 
are Class II, and 1 is a Class III. 



Six of these YSOs: RCrA, IRS5N, IRS7B, SMA2 (RS9), CrA-43 (SMM2), and 
IRAS 32, were detected in the dust continuum at 225 GHz by the SMA. All sources 
except IRAS 32 have emission that is compact and centrally concentrated, assumed 
to be primarily due to compact disk emission. Disk masses for these sources were 
computed and range from 0.009-0.076 Mq. SMA data of IRAS 32 show an extended 
component, assumed to be due to inner envelope emission. The computed disk mass 
for IRAS 32 is 0.024 M©. 

An extinction map created using 2MASS and Spitzer data is presented and compared 
with the distribution of YSOs. The measured Ay over the entire field has a mean of 
~5 mag, with a maximum of ~30 mag, peaking on the Coronet. Using this extinction 
map, we calculated the cloud mass in the area of the e xtinction map above the Ay = 2 



contours to be 279 Mq, confirming the result found in lHeiderman et al.l ( I2OIOI ). There 
is a clear radial spread of sources from the Coronet, with Class I sources clustered most 
tightly in the center, surrounded by Class II sources. Class III sources can be found 
spread throughout the entire field. 
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We perform a clustering analysis mirroring that of iGutermuth et al.l (120091). us i ng the 



116 YSO candidates identified in this paper, and confirm the IGutermuth et al.l (|2009[ ) 
result that the CrA cluster is elongated (having an aspect ratio of 2.36) with a circular 
radius of 0.59 pc and mean surface density of 150 pc~^. By looking at the MST branches 
in the case where all YSOs are included compared with the case where the Class III 
sources have been excluded, it is clear that there is an evolutionary stage gradient in 
CrA, where the older Class III population extends to the west of the Coronet. 

The star formation rate is calculated to be 12 M0Myr~^, similar to that of the Lupus 
clouds. The SFR surface density, however, is quite high: Esfr = 3.96 M0yr~^kpc^^, 
similar to that of Serpens. A disk fraction is also calculated, and found to be 53% for 
CrA. The star formation efficiency for CrA is ~8%, which is the highest of any of the 
c2d/GB regions surveyed, and is perhaps due to the triggering of star formation by an 
expanding H I shell associated with the Sco-Cen superbubbles. 

Our II2 multi-epoch survey maps an area around CrA of about 20'xl7' (i.e. ~10% of 
the Spitzermap), and the continuum-subtracted mosaic reveals nearly 100 H2 knots, 
and we study P.M.s in 44 of them. The derived P.M. values range between 0.024 
and 0.256 "yr^\ corresponding to a Vtan between 15 and 158 kms"^ (at a distance of 
130 pc). 

We identify at least 17 outflows with their candidate driving sources. All the Class O/I 
sources inside the Coronet, and the majority of Class I sources in the CrA region 
drive H2 flows. Most of the known HH objects, and all the detected H2 objects, are 
visible in our IRAC maps. Outflow morphologies in the 3.6 and 4.5 /im bands resemble 
those observed in our H2 images, mostly delineating relatively high-excitation regions, 
whereas the 5.8 and 8.0 /im bands mostly delineate wings and wakes of the bow shocks 
and, in general, milder shock regions. 

There is clear evidence for a parsec-scale precessing outflow, E-W oriented, and orig- 
inating in the SMA2 region: IRS 7A or SMA2 are likely the driving sources. Fur- 
thermore we likely identify the HHlOl lobe counterpart (MHO 2014), which is likely 
driven by IRS 5 or IRS 5N. By means of Spitzer mapping, we identify new flows in 
those regions not covered by our H2 maps: Spitzer outfiows 1, and 2, which are likely 
part of the parsec-scale precessing outflow, E-W oriented, and the IRAS 32 outflow, 
which extends ~0.8 pc and appears to precess as well. 
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A. Appendix A 



A.l. Notes on individual Spitzer-identihed YSO candidates 



In many cases, e.g. CrA-7, the ISOCAM source listed in lOlofsson et al.l (119991 hereafter 
099) corresponds to one of our Spitzer YSO candidates, but is not exactly coincident with 
the position of the Spitzer detection. However, each case was checked and it was found 
that the ISOCAM sources listed as a match in this paper are within the pointing errors 
for both telescopes (and in no case is there another infrared source close enough to deem it 
questionable). 



A. 1.1. CrA-1 



CrA-1 is too far afield to have been observed before; it is not included in the Forbrich fc Preibisch 



( 120071 ) Chandra study. CrA-1 is located in the western part of the CrA molecular cloud called 
"the streamer." It is not detected at 70 fim, but is bright at 24 fim; we classify it as a Class 
II, and it is a new YSO candidate. 
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A. 1.2. CrA-2 (Leda Galaxy 90315) 

CrA-2 was originally classified as a Class I YSO candidate from the 5j>zteer plus 2MASS 
photometry. However, by means of visual inspection, it proved to be a galaxy, and a quick 
search in SIMBAD found it to be Leda galaxy 90315. Interestingly, as was described in Sec- 
tioning! another candidate YSO, CrA-50, is also included within the area which constitutes 
this extended source. Therefore, both CrA-2 and CrA-50 are clearly misclassified as YSO 
candidates and have been removed from our sample. 



A.1.3. CrA-3 (ISO-CrA55) 

CrA-3 was detected with ISOCAM by 099, and is referred to as ISO-CrA55 in their 
survey. From the ISOCAM photometry, it was found to have a mid-infrared excess, and 
therefore is listed as one of their YSO candidates. 

Like CrA-1, it is located in "the streamer" so it has not been observed with Chandra. 
It is not detected at 70 /im, but is quite bright at 24 /xm, and, according to Spitzer colors, 
has a Flat spectrum SED. 



A. 1.4. CrA-4 (DENIS-P J185950. 9-370632, ISO-CrA63) 



CrA-4 is an association member, a young binary brown dwarf in CrA (IBouv et al ."2004h. 



Known as DENIS-P J185950.9-370632, it was discovered to be a binary bv lBouv et aL (20041 
using high resolution Hubble Space Telescope (HST) observations. The source was elongated 
in the WFPC2, ACS and STIS images and the companion comes up clearly on the WFPC2 
and ACS images after the primary PSF was subtracted. From the images, they found it 
to be a common proper motion pair with a separation of 0".060. The magnitude difference 
between the two components indicates a mass ratio of ~75%. 

High and low resolution optical spectra of CrA-4 were also ob tained, and the sp ectral 
features show that it is young, and likely an association member ( iBouy et al.ll2004j ). The 
high resolution spectrum shows Ha in emission, and with a strength that makes it a likely 
accretor; there is also weak Li I absorption, which indicates a young age. 

CrA-4 was detected with ISOCAM by 099, and is referred to as ISO-CrA63 in their 
survey. It is not li sted as a candidate YSO by 099 (i.e. it does not exhibit a mid-infrared color 
excess), however iBouy et al.l ( 120041 ) analyzed the archival images and found there to be a 
mid-infrared excess in the 5.0—8.5 ^m. band. This is consistent with our S^^zfeer observations. 



which detect this source in all IRAC bands and at 24 /im, classifying it as a Class II YSO 
candidate based on the shape of its SED. 



A. 1.5. CrA-5 (ISO-CrA76) 

This source was detected with ISOCAM by 099, and is referred to as ISO-CrA76 in 
their survey. From the ISOCAM photometry, it was found to have a mid-infrared excess, 
and therefore is listed as one of their YSO candidates. 

CrA-5 is classified as a Class I YSO candidate by our method and is outside the Chandra 
field of view. Like CrA-1 and CrA-3, it is not detected at 70 ^m, but is quite bright at 24 /zm. 
Along with CrA-6 (Class II) and CrA-7 (Class III), these three YSO candidates surround 
"Spitzer Outflow 1" which is discussed in more detail in § IB. 31 (see Figure [39!) . 



A. 1.6. CrA-6 (CrAPMSS, CP g2, ISO-CrA 



First observed in the infrared by iGlass fc PenstonI (jl975l ). CrA-6 is liste d as source 
[GP] g2 in their Table 1. It was also detected in X-rays by EINSTEIN (CrAPMS 8: IWalter etal 
19971 ) who also obtained spectra and classified it as M3V. iPattenI (119981 ) classifies it as a 
likely association member: he finds a spectral type of M5, and detects Ha in its spectrum, 
a ROSAT X-ray counterpart, and a VRI color excess. 

Additionally, CrA-6 was detected with ISOCAM by 099, and is referred to as ISO- 
CrA 88 in their survey. From the ISOCAM photometry, it was found to have a mid-infrared 
excess, and therefore is listed as one of their YSO candidates. Most recently, it was one of 
the sources included in a kinematic study of the extended R CrA association (referred to as 
(GP75) RCrAg2 in their Table 13; iFernandez et al.l 120081 ) . and found to be a member. All 
of this is consistent with our classification of CrA-6 as a Class II YSO candidate from its 
Spitzer photometry; we detect this source in all IRAC and MIPS bands out to and in cluding 
70 ixm. Finally, CrA-6 was recently found to have a companion (IKohler et al.l 120081 ). 



A. 1.7. CrA-7 (ISO-CrA 93) 

CrA-7 was detected with ISOCAM by 099, and is referred to as ISO-CrA 93 in their 
survey. From the ISOCAM photometry, it was not found to have a mid-infrared excess. 
However, we classify CrA-7 as a Class HI YSO candidate based on its Spitzer colors, and so 



it is a new YSO candidate. Bright 24 /im emission is detected, but nothing is detected at 70 



A. 1.8. CrA-8 (CrA-444) 



YSO candidate CrA -8 was previously identified as a candidate low-mass member of CrA 
(ILopez Marti et al.ll2005l . their source CrA-444) as pa rt of their optical st u dy of the very low 
mass population of CrA. Using a method outlined in iLopez Marti et al.l (|2004[ ). CrA-8 was 
classified from its optical spectrum as a brown dwarf candidate with spectral type of MS. 5. 
We classify CrA-8 as a Class II YSO candidate from the shape of its SED, detecting it out 
to, and including, 24 fim. 



Lopez Marti et al.l (120051 ) postulate that CrA-444 could have a very close companion. 



which they call CrA-444b. This source, approximately 2". 3 away, is only visible in their deep 
/-band exposure. We do not see any conclusive evidence for CrA-444b in any of our Spitzer 
images, due to the lower spatial resolution. 



A. 1.9. CrA-9 

CrA-9, located to the northwest of the Coronet region and TY CrA, has not been previ- 
ously observed, thus it is a new YSO candidate, classified as a Class II. The source is located 
very close to the edge of the 70 yum map, so there is no catalog entry in this band for it. 
However, the MIPS 70 /im image clearly shows emission, and we estimate the flux using 
aperture photometry, finding a 70 /xm fiux of 154.2 ± 19.8 mJy (using a 2 pixel aperture 
radius) . 



A. 1.10. CrA-10 (CrA-432) 



Lopez Marti et al.l (120051 ) observed this source (referred to as CrA-432 in their Table 3) 



and found a photometric spectral type classification of M7 (1 subclass error on the spectral 
type). Located near CrA-9, we classify CrA-10 as a Class II YSO candidate from the shape 
of its SED; it is detected in all bands out to, and including, 70 fim. 



CrA-10 was also observed by ISicilia-Aguilar et al.l ( l2008l . also called CrA-432) with the 
Spitzer Iniiaied Spectrograph (IRS); the silicate feature at 8—13 nm is not detected in this 
source and there is only a marginal detection of the gas lines in the 13—14 /im IRS band. 
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Another identifier for tliis source is DENIS-P J190059. 7-364711. 



A. 1.11. CrA-11 (HD 176269, HR7169) 



CrA-11, more commonly known as HD 176269, is an optically d etected pre-raain se- 
quence star with a spectral type of B8V that was first observed by iKnacke et al.l (119731 ) 
as HR7169 to be a possible member of the young stellar association near RCrA (~12' to 
the SW). Another B8V star ~12" away, HR7170 (HD 176270), was also hsted as a possible 
member and we include it in T able [71 These two s ources were detected, but not spatially 
resolved, by the IR AS (IRAS 18: lwilking e7aDll992h and EINSTEIN satellites (CrAPMS 10; 
Walter et al]ll997[ l. 



HD 176269 was observed in the near-infrared by lGlass fc Penstoru (119751 ) and is referred 
to in their Ta ble 1 as source /; HD 176270 was also observed, and listed as source k. Later, 
PattenI Jl998h does not classify the HD 176269/HD 176270 pair (which he calls source R05) 
as a likely association member: he finds a spectral type of B9V for the source and detects 
a ROSAT X-ray counterpart, however the source does not fulfill enough criteria for him to 
consider it an association memb er. Note that these stars are far enough afield so that they 
were not included as part of the iForbrich fc PreibischI ( 120071 ) Chandra study. 



HD 176269 was also detected with ISOCAM by 099, and is referred to as ISO-CrAllO 
in their survey. ISO-CrA 110 was found to have a mid-infrared excess based on those ob- 
servations so it is considered one of their YSO candidates. HD 176270 is referred to as 
ISO-CrA 111 in the same study, and it was not found to have a mid-infrared excess. 

These studies are consistent with the Spitzer observations; we classify HD 176269 as a 
Class HI based on the shape of its SED. Although we do not classify HD 176270 as a YSO 
candidate, it very likely is also a Class HI based on its colors and SED shape (see Table [7]). 
The SEDs of the two stars look similar out to 8 /im. However, at 24 //m, HD 176269 is much 
brighter than HD 176270. Neither source is detected at 70 fim. 



A. 1.12. CrA-12 (V667CrA, CrA-4110) 



CrA-12 (or V667CrA; iKukarkin et al.l Il972l) is a kn own variable star, detected with 
Chandra fJ190116.26-365628.4: lForbrich fc Preibischll2007l ) and with ISOCAM (ISO-CrA 123; 
099). CrA-12 was not found to have a mid-infrared excess based on the ISOCAM observa- 
tions. However, based on our Spitzer observations, it is classified as a Class II YSO candidate. 
It does not have a 70 fim detection, but is quite bright at 24 fim. 



-37- 



Lopez Marti et al.l (|2005[ ) observed this source (referred to as CrA-4110 in their Table 3) 



and found a photometric spectral type classification of M5 (erro r of 2 subclasses on the spec - 
tral type). The same spectral classification was also obtained by lSicilia-Aguilar et al.l ( 120081 ) . 
who classify it as a weak-line T Tauri star (wTTs), due to its small Ha equivalent width 
(EW) and the fact that it has no other accretion indicators (see their Table 5). A SpitzerlKS 
spectrum was obtained for CrA-12 but only a short-low (5—14 n m) spectrum was extracted ; 
the silicate feature from 8—13 fj,m was not detected in this source (ISicilia-Aguilar et al.ll2008l ). 



Another identifier for this source is DENIS-P J190116.3-365628. 



A. 1.13. CrA-13 (CrA-466) 



CrA-13 is likely ISO-CrAl27 (099), which is slightly offset (~11") from the CrA-13 
Sj>zteer coordinates. ISO-CrA 127 does fall within the extent of the very bright 24 /im source, 
though, so it is highly likely this is the same source. We detect CrA-13 in all IRAC bands 
and in both the 24 and 70 fim MIPS bands, and classify it as a Class II from its SEP. T his 
source was also detected with Chandra (J190118. 90-365828. 4: iForbrich fc PreibischI 120071 ) . 



Lopez Marti et al.l (120051 ) observed this source (referred to as CrA-466 in their Table 3) 



and found a photometri c spectral type classificatio n of M4.5 (error of 2 subclasses on the 
spectral type). However, ISicilia-Aguilar et al.l (120081 ) . who also obtained an optical spectrum, 
classified it as an M2 based on the strength of its spectral features. Moreover they identify 
it as a classical T Tauri star (cTTs) because it has strong Ha emission (see their Table 5; 
it is also referred to as CrA-466 in their paper). They also classify CrA-13 as a transition 
object (TO) due to i ts lack of infrared excess at wavelengths shorter than 6 /xm. However, 



Ercolano et al 



(I20091) run several models to fit the SED of the source using the spectral type 



of M2 given by lSicilia-Aguilar et al.l (120081 ) . and find that its SED is well fit by an untruncated 
disk model, and is likely not a TO. A Spitzer IRS spectrum was obta ined for CrA-13 and 



shows silicate emission in both the 8—13 fim and 20—30 /im regions (ISicilia-Aguilar et al. 
2008h . 



Another identifier for this source is DENIS-P J190118.9-365828. 



A. 1.14. CrA-14 (HH 101 IRS 1, G-94) 



CrA-14 is an infrared source first observed by iReipurth fc Wamstekerl (|l983l ) to be 
ass ociated with HH 1 01, an d they call it IRS 1. HH 101 was observed again, in the optical, 
bv iHartigan fc Ladal Jl985[ l who called it HH 101 7. 
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This source was also detected with ISOCAM by 099, and is referred to as ISO-CrA 137 
in their survey. ISO-CrA 137 was not found to have a mid-infrared excess based on those 
observations. This is consistent with our Spitzer observations which classify this source as a 
Class III YSO candidate based on its SED; it is faint at 24 and n ot detected at 70 ^^m. 
This source was also detected with Chandra (J190129. 01-370148. 8; iForbrich fc Preibisch 



20071 ). making it a likely YSO candidate. It was observed in multiple epochs with Chandra 



and was fou nd to have some var iable X-ray emission as well as enhanced activity in one of 
the epochs (IForbrich et al.l 120071 ). 



CrA-14 was also observed by ISicilia-Aguilar et al.l (|2008[ ). which they call G-94, who 
obtained an optical spectrum and classified it as M3.5 based on its spectral features. They 
classify it as a wTTs due to its small Hct EW and the fact that it has no other accretion 
indicators (see thei r Table 5). A Spitzer IRS s pectrum was obtained for G-94, but no silicate 
emission was seen (ISicilia-Aguilar et al.l 120081 ). 



Another identifier for this source is DENIS-P J190129.0-370148. 



A. 1.15. CrA-15 (IRS14) 



CrA-15 was first observed by lTaylor fc Storeyl (Il984l . hereafter TS84), and called IRS 14 
(TS 2.9), with the 3.9-m Anglo- Australian Telescope, using its infrared photometer-spectrometer 
sys tem. This sarne sour ce was obs erved with incr easingly better arrays in the ne ar-infrared 
bv lWilking et al.l Jl986l ). and later Iwilking et al.l (source 185809.8-370224; ll997|). However, 
each time, IRS 14 was classified as a field star due to its lack of a near-infrared excess. Based 
on Spitzer photometry, we classify IRS 14 as a Class II YSO candidate; it is not detected at 
70 /im but its detection at 24 /im gives it a mid-infrared excess consistent with a Class II 
SED. 099 saw this mid-infrared excess as well (their source ISO-CrA 139, see their Ta- 
ble 1) and therefore al so considered it a YSO caii didate. IRS 14 was detected with Chandra 
fJ190132.34-365803.1: [Forbrich fc Preibischll2007h . 



IRS 14 was also observed by ISicilia-Aguilar et al.l ( l2008l ). which they call G-87, who 
obtained an optical spectrum and classified it as M3— M4 based on the strength of its spectral 
features compared with other, known M stars. They classify it as a wTTs due to its small 
Ha EW and the fact that it has no other accretion indicators (see their Table 5). They also 
classify IRS 14 as a TO due to its lack of infrared excess at wavelengths shorter than 6 ^m. 
However, as was discussed in Section IA.1.131 for CrA-13, lErcolano et al.l (120091 ) found that 
the SED for IRS 14 was also well fit by an untruncated disk model, and is likely not a TO. 
A Spitzer mS spectrum was obtained for IRS 14 and shows silicate emission in the 8—13 //m 
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regio n as well as marginal detections of gas lines in the 13—14 /im region (jSicilia-Aguilar et al. 
2008h . 



Another identifier for this source is DENIS-P J190132.2-365803. 



A. 1.16. CrA-16 (IRS13) 
CrA-16 was also observed by TS84, and is known as IRS 13 (TS2.8). Li ke IRS 14, this 



source was obser ved with increa singly better arrays in t he nea r-infrared by IWilking et al. 



(119861), and later IWilking et al.l (source 185811.4-370206; ll997|). Unlike IRS 14, though, a 
near-infrared excess was seen for IRS 13 and so it was considered a cTTs and association 



member. IRS 13 was also detected with Chandra (J190133. 84-365745. 0: lForbrich fc Preibisch 
2007h . 



IRS 13 and IRS 14 are about 30" apart, and CrA-17 (discussed in the next section) 
is just a little bit further away. Based on Spitzer photometry, we classify IRS 13 as a 
Class II YSO candidate; it is not detected at 70 /im but its detection at 24 /im gives it 
a mid-infrared exces s consi stent with a Class II SED. An optical spectrum was obtained by 
Sicilia-Aguilar et al.l (|2008[ ). who refer to IRS 13 as G-85, and classify it as M2— M3 based 
on the strength of its spectral features compared with other, known M stars; it also has 
strong Ha emission, confirming it as a cTTs (see their Table 5). A Spitzer IKS spectrum was 
obtaine d, and as for CrA-13, show s silicate emission in both the 8—13 /xm and 20—30 fim 
regions (ISicilia-Aguilar et al.ll2008[ ). 



A. 1.17. CrA-17 



CrA-17 was first observed in the near-infrared by lWilking et al.l (185813.8-370225: Il997l ). 
Based on the Spitzer catalog, we classify CrA-17 as a Class II YSO candidate. The catalog 
lists a 24 fim fiux for CrA-17, but that fiux was obtained by band-filling, and after examining 
the images, it appears that the detec tion may not be real . In a ddition, although CrA-17 
falls within the region covered by the iForbrich fc PreibischI ( 120071 ) Chandra study, it is not 
detected in X-rays. Due to these reasons, we have removed CrA-17 from our list of candidate 
YSOs. 
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A. 1.18. CrA-18 (G-65) 



CrA-18 is l ocated just nor th of TYCrA (see Figu re [2T1) . and was first observed in the 
near-infrared by lWilking et al.l (source 185818.2-365603: 119971 ) . Based on S'pztzer photometry 
we classify CrA-18 as a Class II YSO candidate; it is not detected at 70 /xm but its detection 
at 24 /im gives it a mid-infrared excess consist ent with a Class II SEP. T his source was also 
detected with Chandra f J190140.40-365142.4: [Porbrich fc PreibischI l2007l ) . 



An optical spectrum was obtained by ISicilia-Aguilar et al.l (120081 ) . who refer to CrA-18 
as G-65, and classify the source as Ml— M2 based on the strength of its spectral features 
compared with other, known M stars; it also has strong Ha emission, confirming it as a 
cTTs (see t heir Table 5). Also, a s was discussed in Sections lA. 1.131 and lA.l.TSl for CrA-13 
and IRS 14, lErcolano et al.l (120091 ) found that the SED for this source was also well fit by an 
untruncated disk model, and is likely not a TO. A Spitzer IRS spec trum was obtained for 
CrA-18, but no silicate emission was seen ( ISicilia-Aguilar et al.ll2008l ). 



A. 1.19. CrA-19 (IRS 5, MMS 12, SMM 4 ) 



Best kn own as IRS 5, this was one of the first infrared sources to b e associated with 
CrA ( TS 2.4; iTaylor fc Storeylll984[). IRS 5 was first seen to be a binary by I Chen &: Graham 
(119931 ). and later confirmed by iNisini et al.l ( l2005l ). who found a separatio n of ~78 AU 



between the two components (at a distance of 130 pc). iNisini et al.l (|2005[ ) also found a 



spectral type for IRS5A of K5— K7V (IRS5B was too faint to classify). Additionally, IRS 5 
was detected with the VLA at 6 cm fcal led VLA 7:lBrownl 19871). and shows s i gnific ant short 



and long-term variability in the radio (ISuters et al.lll996l ). iFeigelson et al.l ( 1l998l ) detected 



circularly polarized continuum emission at centimeter wavelengths in IRS 5, demonstrating 
that radio emission from protostars can sometirnes ar ise from non-thermal processes. IRS 5 
is also a millimeter sou rce, MMS 12 ( IChini et al.ll2003l ). and was observed at 450 and 850 /im 
by iNutter et al.l (120051 referred to in their survey as SMM 4), seen as an extende d sour ce in 
their maps. Using photometry available in the literature at the time, IChen et al.l ( 119971 ) use 
a distance of 130 pc to estimate Thoi = 403 K and Lhoi = 0.9 Lq for IRS 5 (which they refer 
to as TS2.4 in their Table 2). 

The Chandra X-ray emission for the binary is only margin ally resolved, so IRS 5 A/B i s 
listed as a single source in the catalog of iForbrich fc PreibischI (J190148. 02-365722. 4; l2007l ). 
However, a later study by iHamaguchi et al.l ( 120081 ) used a subpixel repositioning technique 
in order to resolve the X-ray emission from multiple Chandra observation s. They found tha t 
IRS 5A is fiaring, while IRS 5B is quiescent. IRS 5A/B is also resolved by lChoi et al.l ( 120081 ) . 



-41 - 



called CHLT3 (IRS 5 A) and CHLT4 (IRS5B), in their centimeter imaging survey of CrA. 

With Spitzer, we classify IRS 5 as a Class I YSO candidate; it is detected in all bands 
out to, and including, 24 fim. We do not resolve the binary. In addition, IRS 5 is one of the 
sources which we observed with the SMA (see Section [5]), and surprisingly, there is no clear 
detection, although IRS 5N does show strong detection (see Figure [H] and the following, 
^ IA.1.20p . 



A. 1.20. CrA-20 (IRS5N) 



IRS 5N was first seen in the radio ( iForbrich et al.ll2006l . listed as sou rce 5), and then later 
in X- r ays with Chandra, an d referred to as IRS 5N (J190148.46-365714.5: lForbrich fc Preibisch 



20071 : iForbrich et al.l 120071 ) . Based on Spitzer photometry, we classify IRS5N as a Class I 



YSO candidate; it is detected in all bands out to, and including, 24 /im. IRS 5N was one of 
the sources we observed with the SMA, and it has a strong detection (see §[S]and Figured!]). 
We compute a disk mass for IRS5N of O.O23M0 (see Table [9]). 



A. 1.21. CrA-21 (IRS 8) 



CrA-21 was also observed by TS84, and is known as IRS 8 (TS2.2). Like IRS 13 and 
IRS 14, this source was obser ved with incre asingly better arrays in the n ear- infrared by 
Wilking et allJl986l ). and later Iwilking et al.l (source 185828.8-365834;ll997|). Unlike IRS 14, 
though, a near-infrared excess was seen for IRS 8 and so it was considered an association 
member. We classify CrA-21 as a Class II YSO candidate; it is not detected at 70 /xm but 
its detection at 24 /xm gives it a mid-infrared excess consistent with a Class II SED. 



Meyer fc WilkingI ( 120091 ) also classify this source as a likely association member from H- 



and K-band spectra, finding it to be consistent with a G5 star. A dditionally, it was detected 
with Chandra rjl90151. 11-365412.5: librbrich fc Preibischll2007h . 



A.l. 



CrA-i 



CrA-22, which is located directly south of the Coronet, i s too far afield to have been ob- 
served before, nor is it included in the IForbrich fc Preibisch! (120071 ) Chandra study. Another 
source falls very close to CrA-22 in the IRAC bands, however the catalog flux detections of 
CrA-22 are good, and do not appear to contaminate the photometry. CrA-22 has a MIPS 



detection at 24 fim, but is not detected at 70 /xm. We classify CrA-22 as a Class II, and it 
is a new YSO candidate. 



A. 1.23. Cr A- 23 (Star A) 



First detected by iGrahanj (119931 ). this source, referred to as "Star A" (which was not 
detected by IRAS), was found to be very bright in Ha, weaker in i?— band, and very weak, 
only just visible in the original S II image. They recognized it as a likely a strong emission- line 
star, but it was too faint to obtain a spectrum. 



Later, Star A was also detected in the infrared J, H, fC'-bands by IWilking et al.l (119971. 



source 185831.1-370456), who classified it as a brown dwarf candidate. iLopez Marti et al. 
(120051 ) found a photometric spectral type classifi cation of M8.5 (1 subclass error on the 



spectral type) for Star A (they call it CrA-465). ISicilia-Aguilar et al.l (120081 ) obtained an 
optical spectrum (also referred to in their paper as CrA-465), classifying it as M7.5 based on 
the strength of its spectral features; it also has strong Ha emission which is double peaked, 
confirming it as a cTTs (see their Table 5). With Spitzer, we detect Star A in all IRAC 
bands and at 24 /um, but not at 70 //m. We classify it as a Class II YSO candidate based 
on its SED. 



Another identifier for Star A is DENIS-P J190153.7-370033. 



A.l.^ 



CrA-24 



We classify CrA-24 as a Class I YSO candidate based on its Spitzer SED. It is definitely 
detected up through 8 /im with a rising SED, however it is difficult to tell whether or not the 
24 /im MIPS detection is an actual detection due to contamination from the nearby R CrA 
nebula. It is not dete cted in our SMA data (see Figure [T^. but it is observed by Chandra 
fJ190155.61-365651.1: [Forbrich fc PreibischI l2007h . Due to its position, CrA-24 could be a 
good candidate driving source for HH99 (MHO 2000 A; see Figure [T8|) . Thus we include it 
as a new YSO candidate. 



A. 1.25. CrA-25 



CrA-25 is classified as a Class I YSO candidate based on its Spitzer SED, but that 
SED is based on almost all band-filled fiuxes. Looking at CrA-25, it is clear that it is an 
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extended source, a knot that is associated with an HH object, located within 3" of HH732A 
( jWang et aLll2004j ). Therefore CrA-25 has been removed from our hst of candidate YSOs. 



A. 1.26. CrA-26 (WMB 185844.2-370304) 



CrA-26 was detected by IWilking et al.l ( 119971 ) in the infrared J, H, X'-bands (named 
185844.2-370304). It is classified with Spitzer as a Class II YSO candidate and is well 
detected through 8 fim. There is a flux listed at 24 /im, but it is suspect, due to the fact 
that it falls within a diffraction spike; the 24 /im flux has been band-filled. After inspecting 
the SED (see Figure [7]), we determine that it looks like a YSO candidate out to, and including 
8 /im. Therefore we retain CrA-26 as a candidate Class II YSO. 



A. 1.27. CrA-27 (WMB 185848.1-365808) 



CrA-27 was detected by IWilking et al.l ( 119971 ) in the infrared J, H, X'-bands (named 
185848.1-365808). Also detected by 099 (ISO-CrA 142), it did n ot exhibit a mid-i n frared 
color excess and was not listed as a candidate YSO. However, iMeyer fc WilkingI (120091 ) 
obtained an infrared spectrum and determined it to be a likely association member, "late- 
type" > K4. With Spitzer, we classify CrA-27 as a Class III YSO candidate, which is 
consistent with the other studies. 



A. 1.28. CrA-28 (ISO-CrA 143, G-14) 



CrA-28 was detected by IWilking et al.l ( 1l997l ) in the infrared J, H, K'-bands (named 
185849.3-370733). It was also detected with ISOCAM by 099 (ISO-CrA 143), was found to 
have a mid-infrared excess, and therefore is listed as one of their YSO candidates. CrA-28 
is similar to CrA-30 (see § lA. 1.3011 because they both have small H — K colors and similar 
SEDs (099). On the basis of our Spitzer data, we classify this source as a Class II YSO 
candidate; it is detect ed out to, and including, 24 /xm. CrA-28 also has a Chandra detection 
(J190211.99-370309.4; IPorbrich fc PreibischI EoOtI ) . 



This source was also observed by lSicilia-Aguilar et al.l ( 120081 ) . called G-14, who obtained 
an optical spectrum and classified it as an M4.5 based on the strength of its spectral features. 
They classify it as a wTTs because it has weak/no Ha emission (see their Table 5). They 
also class ify it as a TO due to its lack of infrared excess at wavelengths shorter than 6 /im. 
However, lErcolano et al.l ( 120091 ) run several models to fit the SED of G-14 using the spectral 



-44- 



type of M4.5 given by ISicilia-Aguilar et al.l (|2008[ ). and find that its SED is well fit by an 
untruncated disk model, and is likely not a TO. A Spitzer IRS spectrum w as obtained and 
it shows silicate emission in the 8—13 fim region (ISicilia-Aguilar et al.ll2008l ). 



A. 1.29. CrA-29 (ISO-CrA 145) 



CrA-29 was detected by IWilking et al.l ( 119971 ) in the infrared J, H, iT'-bands (named 
185852.0-370456). It was also detected with ISOCAM by 099 (ISO-CrA 145), was found to 
have a mid-infrared excess, and therefore is listed as one of their YSO candidates. With 
Spitzer, we classify this source as a Class II YSO candidate, and it is detected out to, and 
including, 24 jj,m. 



A. 1.30. CrA-30 (Ha U; HBC680) 



CrA-30 was detected as an emission-line star bv lMarraco &: RvdgrenI (Il98ll Ha 14) and 
is also listed in the Herbig & Bell Catalog as HBC 680 Jflerbig fc Belllll988l). It has also been 



noted as an IRAS detection, source IRAS 18591-3702. fe.g. IWeintraublll99ol: IWilking et al.l 



1992). Using photometry available in the literature at the time, Chen et al. ( 19971 ) 
distance of 130 pc to estimate Ti,oi = 2604 K and L^oj = 1.2 for Ha 14. 



use a 



This source was detected by IWilking et al.l ( 1l997l ) in the infrared J, H, X'-bands (named 
185904.3-370238), and by 099 (ISO-CrA 155). From the ISOCAM photometry, it was found 
to have a mid-infrared excess, and therefore is listed as one of their YSO candidates. Ha 14 
has a Chandra detection (J190227 .05-365813.2; Forbrich fc Preibischll2007l ). and in addition, 
it was found to have a companion (IKohler et al.ll2008l ) . A Spitzer IRS sp ectrum was obtained, 



and the source, referred to as G-1, shows silicate emission at 8—13 /xm (ISicilia-Aguilar et al. 
2008h . 



Most recently. Ha 14 was one of the sources included in a kin ematic study of the ext ended 
RCrA association (referred to as HBC 680 in th eir Table 13: iFernandez et al.ll2008l ). and 
found to be a member; it was also observed by iMeyer fc Wilking] ( l2009l ) and determined 
from its infrared spectrum to be a MO member of the association. All of these observations 
are consistent with our Spitzer observations which show it to be well-detected out to, and 
including, 24 /im, classifying it as a Class II YSO candidate from its SED. 



-45 - 



A. 1.31. CrA-31 (Ha 3; ISO- CrA 1 59 ) 



CrA-31 is an emission-line star of M type (IMarraco fc Rydgrenlll98ll . Ha 3). IChen et al. 

(119971 ) use a distance of 130 pc to estimate Tboi = 1875 K and Lboi = 0.6 L© for Ha 3. This 



source was also detected with ISOCAM by 099 (ISO-CrA 159), was found to have a mid- 
infrared excess, and the refore is listed as one of the ir YSO candidates. In addition, it was 
detected with Chandra (IForbrich fc PreibischI 120071 . catalog number J190233. 07-365821.1). 
Our Spitzer observations classify this source as a Class II YSO candidate and it is detected 
all the way out to, and including, 70 /xm. 



A. 1.32. CrA-32 

We originally classified CrA-32 as a Flat spectrum YSO candidate based on its Spitzer 
SED. However, its MIPS 24 fim flux has been band-filled and is likely not very accurate 
due to the fact that CrA-32 is close to the very bright 24 /im source, CrA-44 (also known 
as IRAS 32, see § lA. 1.441) . In the IRAC bands, there is also another source in very close 
proximity making CrA-32 look quite extended. Due to the complicated nature of this source, 
CrA-32 has been removed from our list of candidate YSOs. 



A. 1.33. CrA-33 (IRAS 32d) 



CrA-33 is classified as a Flat spectrum YSO candidate based o n its 5'pzteerSEp, an d is 
well detected out to, and including, 24 /im. It was first observed by lWilking et al.l (119921 ) as 
part of the IRAS 32 outflow, and is referred by to them as IRAS 32d. It was also detected 
with ISOCAM by 099, and is referred to as ISO-CrA 185 in their survey. From the ISOCAM 
photometry, it was found to have a mid-infrared excess, and therefore is listed as one of their 
YSO candidates. Located to the southeast of the main IRAS 32 outflow (see Figures [15] and 

CrA-33 is too far afield to have been observed by many other surveys, e.g., it is not 
included in the iForbrich fc PreibischI (120071 ) Chandra study. 



A. 1.34. CrA-34 (ISO-CrA 198) 



This source was detected with ISOCAM by 099, and is referred to as ISO-CrA 198 in 
their survey. From the ISOCAM photometry, it was found to have a mid-infrared excess, 
and therefore is listed as one of their YSO candidates. 
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Located to the north of IRAS 32, CrA-34 is to o far afield to have b e en ob served by 
many other surveys, e.g., it is not included in the iForbrich &: PreibischI ( 120071 ) Chandra 
study. Based on the shape of its Spitzer + 2MASS SED, we classify it as a Class III YSO 
candidate and it is detected in all bands up to, and including, 24 fim. 



A. 1.35. CrA-35 (ISO-CrA201) 

This source was detected with ISOCAM by 099, and is referred to as ISO-CrA201 in 
their survey. From the ISOCAM photometry, it was found to have a mid-infrared excess, 
and therefore is listed as one of their YSO candidates. 

Located very close to IRAS 32, just to the southeast, CrA-35 is too far afield to have 



been observed by many other surveys, e.g., it is not included in the iForbrich &: Preibisch 



(120071 ) Chandra study. Based on the shape of its Spitzer + 2MASS SED, we classify it as a 



Class II YSO candidate and it is detected in all bands up to, and including, 24 //m. 



A. 1.36. CrA-36 



Also located to the southeast of IRAS 32, CrA-36 is too far afield to have been observed 
by many other surveys, e.g., it is not included in the IForbrich fc PreibischI (l2007l ) Chandra 
study. Based on the shape of its Spitzer + 2MASS SED, we classify it as a Class II YSO 
candidate and it is detected in all bands up to, and including, 70 /im. CrA-36 has not been 
classified as a YSO candidate before, therefore it is a new YSO candidate. 



A. 1.37. CrA-37 



CrA-37 was possibly detected with ISOCAM by 099, referred to as ISO-CrA232 in their 
survey. ISO-CrA 232 is ~8" away from CrA-37 which is close to the pointing uncertainties 
for both telescopes; there are no other infrared sources nearby so it is likely that this is 
the same source. From the ISOCAM photometry, ISO-CrA 232 was not found to have a 
mid-infrared excess, and therefore is not listed as one of their YSO candidates. 

CrA-37 is too far afield to have been observed by many other surveys, e.g., it is not 
included in the IForbrich fc PreibischI ( l2007l ) Chandra study. It is located to the southeast of 
IRAS 32. Based on the shape of its SED, we classify it as a Class I YSO candidate and it is 
detected in all bands up to, and including, 70 /im. CrA-37 has not been classified as a YSO 
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candidate before, therefore it is a new YSO candidate. 



A. 1.38. CrA-38 (IS0-CrA13) 



This source was detected with ISOCAM by 099, and is referred to as ISO-CrA 13 in 
their survey. From the ISOCAM photometry, it was found to have a mid-infrared excess, 
and therefore is hsted as one of their YSO candidates. 

CrA-38 is also too far afield to hav e been observed by many other surveys, e.g., it is 
not included in the lForbrich fc Preibischl (120071 ) Chandra study. Like CrA-1 and CrA-3, it is 
located in the western part of the CrA molecular cloud, "the streamer." Based on the shape 
of its Spitzer + 2MASS SED, we classify it as a Class III YSO candidate and it is detected 
in all bands up to, and including, 24 /xm. 



A. 1.39. CrA-39 (6dFGS gJ190023. 5-371224 ) 

CrA-39 was originally classified as a Class II YSO candidate from the Spitzer + 2MASS 
photometry. However, a closer look at the images reveals that it is a galaxy (indicated in 
SIMBAD to be an emission-line galaxy, 6dFGS gJ190023.5-371224). Therefore, this source 
is clearly misclassified as a YSO candidate and has been removed from our sample. 



A. 1.40. CrA-40 (VSStlS, ISO-CrA 134) 



CrA-40 was first observed at 2.2 fim by IVrba et al.l (Il976a] . known as VSSt 18) and 
again later, bv lwilking et al.l (119971 1 in the infrared J, K'-bands (named 185802.9-370344). 
CrA-40 was also detected with ISOCAM by 099, and is referred to as ISO-CrA 134 in their 
survey. From the ISOCAM photometry, it was found to have a mid-infrared excess, and 
therefore is listed as o ne of their YSO candidates. Additionally, it was detected with Chandra 
(J190125. 75-365919.3: IForbrich fc Preibischl l2007h . 



Most recently CrA-40 was classified by iMeyer &: Wilkingl (120091 ) from its infrared spec- 
trum as a K3 member of the association. All of this is consistent with our Spitzer observa- 
tions, which detect this source in all IRAC bands and at 24 fim, classifying it as a Class II 
YSO candidate. 
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A.1.4i- CrA-41 (Ha2; HBC677) 



CrA-41 is an emission- line star, first observ ed by iMarraco &: RydgrenI (jl98ll Ha 2), and 



listed in the Herbig & Bell Catalog as HBC 677 ( iHerbig &: Belllll988). This infrared source is 



known to be a part of HH 100, which was obse rved in the opt i cal bylHartigan fc Ladal (119851 ) 
who called it HH 100 1, and was detected by IWilking et al.l ( 119971 ) in the infrared J, H, K'- 
bands (named 185819 1-37 0418). CrA-41 also has a C/ianc?ra detection (J190141.62-365953.1; 
Forbrich k PreibischI 120071 ) . 



Most recently, CrA-41 was one of the sources included in a kinematic study of th e 



extended R CrA association (referred to as HBC 67 7 in their Table 13: 



^ernandez et al.ll2008l ) 



and found to be a member; it was also classified by lMeyer &: Wilkine J2OO9I ) from its infrared 
spectrum as a M2 member of the association. All of these observations are consistent with 
our Sj»zfeer observations, which detect this source in all IRAC bands and at 24 ^m, classifying 
it as a Class II YSO candidate based on the shape of its SED. 



A. 1.42. CrA-42 (IRS6) 



CrA-42 is a binary source, known as IRS 6 A an d 6B, first no t ed as a binary by lNisini et al 



( 120051 ): we do not resolve it in our Spitzer images. iNisini et al.l ( 120051 ) found a separation of 
97 AU (< 1") between the two components (at a distance of 130 pc), and a spectral type for 
IRS 6A of M2V (IRS 6B was too fai nt to classify ) Pos itionallv, our YSO candidate CrA-42 
is closer to the coordinates given by lNisini et al.l ( 120051 ) for IRS 6A, but we will discuss both 
IRS6A and 6B here (referring to it as "IRS 6," since it is unresolved in most studies). 



IRS 6 was first observed by TS84 ( TS 2.3 ), and then by IWilking et al.l (Il986l ). and later 
Wilking et al.l (source 185828.2 -370058: 119971 ). Using photometry available in the literature 
at the time, IChen et al.l ( 119971 ) use a distance of 130 pc to estimate Tboi = 1055 K and Lhoi 
= 0.1 Lq for IRS 6 (which they refer to as TS2.3 in their Table 5). 

Along with several Class I sources in the Coronet cluste r (including IRS l, IR S 2, IRS 5, 
IRS 7, IRS 9, and RCrA), IRS 6 was first seen in X-rays bv lKovama et al.l (119961) with the 
ASCA satellite. Interestingly, it was not detected in X-rays with ROSAT (iNeuhauser fc Preibisch 



19971) . although the ROSAT band has a softer energy range than ASCA. iForbrich et al. 
(j2006 ) find that IRS 6 is detected as a weak X-ra y source with Chandr a and XMM-Newton 
and it is referred to as IRS6A in Table 2 of Forbrich fc Preibisch (source J190150.45- 



365638.1; 12003) 



Wilking et al, 



( I1997I ) noted that IRS 6 is associated with a reflection nebula, and 



Wang et al, 
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( 120041 ) infer that it could be the driving source for a giant outflo w including HH objects 
HH99, 730, 860, and the HH 104 D jet (see Figure [28]). However, ICaratti o Garatti et al. 
( 120061 ) discount IRS 6 as the driving source for HH 99 based on the tangential velocity and 
position angle of HH99 (see their Section 7.11 and their Figure 15). On the basis of our anal- 
ysis, we suggest that IRS 6 is driving a large outflow which includes MHO 2000 B, HH99 C, 
and MHO 2004 B-D (see also Appendix [B] and Table [13]). 



I n the radio, IRS 6 has b een seen as a very weak source at 3.6 cm and 6 cm (IForbrich et al. 

20061: iMiettinen et al.ll2008l ). Unlike IRS5A/B, IRS6A/B was not resolved bv lchoi et al. 
(120081 ) in their centimeter imaging survey of CrA, in which they refer to IRS 6 as CHLT 5. 
In fact, like the previous studies, the single source is only weakly detected in their 3.5 cm 
continuum map, and flux values for densities at both 6.2 and 3.5 cm are not listed in their 
Table 3. 



Most recently, IRS 6 was classified by iMeyer &: Wilking] (120091 ) from its infrared spec- 
trum as a Ml member of the association. All of these observations are consistent with our 
Spitzer observations, which detect this source in all bands up to, and including, 24 fim, and 
classify it as a Class II YSO candidate. 



A. 1.43. CrA-43 (SMM2, CHLT 15, WMB 55) 



CrA-43 was first observed by I Wilking et al.l (119971 ) in the in frared J, H, -K '^ -band s (named 
185836.1-370131). This source co rresponds to SMM 2 from the iNutter eTaP ( l2005h SCUBA 
450 and 850 /xm survey of CrA. iNutter et al.l (120051 ) also refer to the near- infrared coun 
terpar t of SMM 2 as WMB 5 5 because it is source number 55 in Table 1 of IWilking et al. 



(119971 ). iGroppi et al.l (120071 ) mention this source within their study of SMA observations 



of the IRS 7 region (of which SMM 2 is nearby), however their field did not cover the 
SMM 2 region. Although C rA-43 is part of the field covered by Chandra, it is not detected 
jForbrich fc PreibischI [2007f ) . 



In lHaas et al.l (120081 ). CrA-43 is identified as a candidate YSO using a millimeter excess 
technique, called A3 51, with a 1. 2 mm flux o f 101.8±18.4 MJy/sr. It i s also weakly detected 
in 3.5 cm continuum emission by lChoi et al.l (referred to as CHLT 15; 120081 ) . but not in the 
6.2 cm continuum. A close-up of CrA-43 and the surrounding flo ws can be seen in F igures Wl 
and 1231 as well as in Figure 5 ([SII] narrow-band images) of I Wang et al.l (120041 ). clearly 
indicating that it is driving an outflow (see also Appendix [B] and Table [T3|) . 



In our Spitzer images, CrA-43 is well-detected out to, and including 24 //m, and is 
classifled as a Class I YSO candidate. It is also detected with the SMA (see Figured!]), and 
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we compute a disk mass for CrA-43 of O.O32M0 (see Table [9]). 



A. 1.44. CrA-44 (IRAS 32c) 



CrA-44 is classified as a Class I YSO candidate based on its Syitzer SEP, and is well 
detected out to, and including, 160 /im. It was first observed by lWilking et al.l (119921 ) as part 
of the IRAS 32 outflow and is referred to by the m as IRAS 32c (IR AS 18595-3712). Using 
photometry available in the literature at the time. lChen et al.l ( 119971 ) use a distance of 130 pc 
to estimate = 148 K and L^oi = 1.3 Lq for IRAS 32 (see their Table 2). It was also 
detected with ISOCAM by 099, and is referred to as ISO-CrA 182 in their survey. From the 
ISOCAM photometry, it was found to have a mid-infrared excess, and therefore is listed as 
one of their YSO cand idates. IRAS 32 has a near-infrared ref lection nebu l a first reported by 
Connelley et al.l (120071 ) . and named the "Isabelle-Nebula" by iHaas et al.l (120081 ) . which they 
observed in the H- and K^-bands. This reflection nebula can be easily seen in the infrared 
Spitzer bands and ISAAC/ VLT narrow-band 2.12 /xm image, and it likely delineates the 
outflow cavity (see Figure H2|) . 



There have been other Spitzer studies of IRAS 32 as well. ISeale fc Looneyl (120081 ) show 
that the outflow cavity has an opening angle averaging 66° ± 4. Bec ause IRAS 32 is nearly 
edge-on, when correcting for inclination angle, ISeale fc Looneyl ( 120081 ) find a similar average 
opening angle of 6 5° ± 4. In addition , a Spitzer IRS spectrum of IRAS 32 shows a solid CH4 
feature at 7.7 /im (lOberg et al.ll2008l ). 



IRAS 32 is also a millimeter source, MMS23 (IChini et al.l l2003l ). and was observed 
at 450 and 850 fim by iNutter et al.l (120051 . re ferred to in their survey as SMM8). Using 
the Atacama Pathfinder Experiment (APEX). Ivan Kempen et al.l (120091 ) observed IRAS 32 
in the submillimeter lines of CO, HCO^ and their isotopologues, finding it to be a truly 
embedded source with an exceptionally strong outflow. This is confirmed by our SMA, II2, 
and Spitzer observations (see §13 Appendix |B| and Figure USD, which show a parsec-scale 
precessing outflow. Using the SMA data, we compute a disk mass for IRAS 32 of O.O24M0 
(see Table |9]) . 



IRAS 32 is too far afield to have been observed by many other surveys, e.g., 
included in the iForbrich fc PreibischI (12007! ) Chandra study. 



it is not 
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A. 1.45- CrA-45 (VSStlO, IRAS 34) 



CrA-45 was first observed at 2.2 um by IVrba et al.l (1976al . known as VSSt 10), and is 
an IRAS source, named I RAS 34 by IWilking et al.l ( 119921 ) . Using photometry available in 
the literature at the time, IChen et al.l ( 119971 ) use a distance of 130 pc to estimate Tboi = 
712 K and Uoi = 0.5 Lq for VSSt 10 (see their Table 2). 

Located to the southeast of IRAS 32 (CrA-44), CrA-45 is on the outskirts of the main 
CrA cloud core, a nd is the closest YSO candidate to the quiescent "Core C" discussed in 
Harju et al.l (119931 ) at about ~8' to the northwest. It is too far afield to have been observed 



by many other surveys, e.g., it is not included in the iForbrich fc PreibischI (120071 ) Chandra 
study. Our Spitzer observations classify CrA-45 as a Class II YSO candidate, and it is well 
observed out to, and including, 70 /im. 



B. Appendix B 

B.l. Notes on individual MHOs - North and east regions, Coronet. 

B.1.1. Morphology 

Figure [18] shows a SofI II2 image of north-eastern flows detected outside and inside the 



Davis et al. 


1999; 


Caratti Garatti et al. 


2m. 


e. IIH733, see 


Wilkine: et al. 


1997; 


Davis et al. 



I999I : IWang et al.l 120041 ) . and four newly detected outflows, namely MH O 2009, 201 



MHO 2011, and 2014 are partially coincide nt with HH733, and 732 (see 



respectively. Finally, MHO 2001 and 2002 (ICaratti o Garatti et al.ll2006 



Wang et al 



Davis et al. 



-2014. 



2004) 



2OIOI) 



positioned inside the Coronet are labeled with blue crosses in Figure [TSl 

A close up view of MHO 2009, 2011, and 2012 is given in Figure [IS MHO 2009 is 
composed of a chain of precessing knots (A-F), roughly elongated west-east for ~ 1.2'. About 
20" northward, MHO 2011 (knots A-G) shows a curved shape, elongated towards ENE, and 
is more extended (~ 3.6'). Knots D and F are the brightest structures. The first is a large 
bow-shock, coincident with HH735, while the second has a point-like shape. 

Roughly 25" northward, we detect three faint aligned knots (MHO 2011 A-C) oriented 
SW-NE (length ~1.7'). HH736, not detected in our H2 images, is positioned farther north- 
eastward, and although not well aligned with the knots, could be part of the fiow. Further 
north (see Figure ITS!) we detect MHO 2013 and 2000, molecular counterparts to HH733 and 



- 52 - 



HH99. A close up view of the two flows is given in Figure |20l Both flows are roughly 
parallel and elongate d towards the NE. M HO 2013 shows a morphology identical to the [SII] 
emission observed by IWang et al.l ( 120041 ). with three brighter condensations located at the 
curved head of the jet. MHO 2000 is composed of two bright structures, a heading bow- 
shock (A - HH99B) and a knot (B - HH99A), with a rough bow-shape, oriented towards 
MHO 2013. It is not clear whether or not this last is associated with the main flow or with 
the MHO 2013 flow, although our P.M. analysis marginally indicates that the latter is the 
case (see Appendix IB. 1.21 and Figure | 



About 5' northward of the Coronet and ~3.5' to the east of TYCrA, we discover 
MHO 2014, the brightest MHO in CrA (see Figure EB- The MHO extends for roughly 
1', and displays an extremely complex morphology, as shown in Figure [5T], where H2 and H2 
continuum-subtracted images are reported. Three arc-shaped structures are visible in the 
eastern part of the object (knots A-C), while the brightest knots (D-F) are located in the 
western part, partially coincident with HH732A and C. Such a puzzling structure could be 
caused by the interaction of several bow-shocks from a precessing flow, or the overlapping 
of different flows. Notably, HH732B, not detected in H2, is located ~30" to the north, 
suggesting that more than one flow could be present. 

Flow detection and P.M. study inside the Coronet region is particularly difficult, due 
to the presence of RCrA and the extremely bright and diffuse nebulosity around it, which 
cannot be completely removed in our H2 continuum-subtracted images (see e. g. Figure | 



Nutter et al. 


2005; 


rial 112 iviiuuD 

Groppi et al. 


2007) 



(]2006[ ). and named MHO 2001 and 2002 by IPavis et al.l fcOlOf ). Figure |22] (upper panel) 
shows a SofI H2 continuum-subtracted image of the region, where we detec t more knots 



posit ioned a few a rcseconds NE and E of the SMA 2 and IRS 7A region (see e. g. , iGroppi et al. 



20071 ). Following IPavis et al.l (]2010l ). the new knots are named MHO 2001 from K to Q, 
although it is clear that MHO 2001 and 2002 consist of several overlapping flows. Knot Q 
is the only source of H2 emission detected westward of RCrA. MHO 2001 is an elongated 
chain of scattered knots (roughly directed SW-NE), with several knots grouped around the 
SMA 2 and IRS 7A regions, which are the focal point of several flows. A closer look with 
ISAAC/VLT is given in Figure |22] (lower panel). The H2 image is shown in logarithmic 
scale to display simultaneously the bright regions close to SMA and the fainter knots of 
MHO 2002. Unfortunately, an ISAAC H2 continuum-subtracted image is not available for 
this region. The ISAAC/VLT image reveals the presence of at least five different flows in 
this region, possibly more. Knots C, I, K, P, O are possibly part of a flrst flow directed 
NNE (displayed with white arrows). A second flow could consist of MHO 2001 A-E, L, 
N , directed NE (cyan arrows). MHO 2002 A is the terminal bow-shock of a third outflow 
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(roughly west-east oriented, P.A.~90°) delineated by green arrows in the figure. It is not 
clear whether MHO 2002 C could be part of this fiow, or, as a faint H2 trail seems to suggest 
could be driven by IRS 7B (yellow arrows). MHO 2001 J is the jet of MHO 2002 A, and knot 
M could be part of t he flo w as well (green arrows). This flow matches well the centimeter 
maps by IChoi et al.l (120081 ). who detected a west-east jet l aunched by SMA2. This is also 
observed in the CO 3-2 maps by (Ivan Kempen et al.l 120091 ). who detected a west-east flow 
(P.A.~ 90°). The blue lobe is positioned to the east of SMA2 region, while the red lobe 
is to the west. Another flow (red arrows), SE-NW oriented (P.A.~125°), and composed of 
MHO 2002 B, MHO 2001 G, and possibly H, seems to originate from CrA-43 (SMM2). 



B.1.2. P.M.s 

A quick inspection of Figure [T71 conflrms that all the studied knots, except MHO 2008 D, 
are launched by YSOs inside or close to the Coronet. As already noted in § 16. H several flows 
originate in the eastern part of the Coronet, in particular in a region a few arcse conds around 



SMA2, which encloses s everal early stage YSOs (see Appendix |X] and e.g., iNutter et al. 



20051 : iGroppi et al.ll2007f ). This can be seen in more detail in Figure El MHO 2002 A is 
the terminal bow-shock of a west-east flow (MHO2001J, M), which has MHO 2001 Q as 
a flow counterpart. MHO 2002 C is likely part of another flow from IRS7B. MHO 2002 B 
and MHO 2001 G are the south-eastern and north-western lobes of the CrA-43 (SMM2) 
outflow. Then MHO 2001 D-N are part of one or more outflows roughly moving towards the 
north-east and emanate from the SMA 2 region. 

Figure m shows the flow charts of MHO 2009, 2011, 2012 and part of MHO2001 A and 
B, which move from west to east and from WSW to ENE, being indeed the prosecution of 
the flows of Figure El Figure ES shows P.M.s in MHO 2000 and 2013. Our multi-epoch 
images allow us to get a better measurement of the P.M. for MHO 2000 A, which appears to 
originate from a region around IRS 9 and CrA-24. On the other hand, MHO 2000 B is moving 
eastward and might not be related to this flow, but rather to MHO 2013, which appears as 
an expanding shell moving towards ENE, in this case IRS 6 is likely the driving source. 

MHO 2014 (see Figure [26!) is roughly moving towards NNE, and it is likely an expanding 
bow-shock, which originates from the western outskirts of the Coronet. Unfortunately, the 
lack of a second-epoch image for knots D and E does not allow us to completely reconstruct 
its motion. 
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B.2. Notes on individual MHOs - South and west regions 



B.2.1. Morphology 

Figure [771 shows a SofI H2 image of flows located to the south-west of the Coronet. There 
are four MHO flows, roughly NE-SW oriented, escaping from the Coronet (MHO 2003, 2004, 



2006- 2008). MHO 2003, 2004 flWilking et al.lll997l : ICaratti o Garatti et al.ll2006l : lOavis et al. 



2OIOI ) and the newly detected MHO 2006, 2007 and 2 008 are the molecular counterparts of 
HHIOO, HH104, HH97 HHlOl, and HH 731 fsee e. g.. [strom et al.lEoTi : iHartigan fc Graham 
19871 : IWang et al.ll2004h . respectively. 



Figure [28] shows an ISAAC H2 continuum-subtracted image of MHO 2004 (i. e. HH 104 C 
and D). The molecular object, positioned about 20" SW of IRS 6, appears as a faint curved 
jet with three bright condensations (namely B-D), NE-SW oriented, and a brighter bow- 
shock (A or HH 104 C) placed ~8" NW off-axis from the jet and east-west oriented. The 
bow-shock should not be associated with the main jet, which is likely driven by IRS 6. 

MHO 2003 (see Figure [27]) is positioned about 20" SW of IRS 1 (or HH 100 IR), which is 
likely the driving source. A close up view of the flow (see Figure [2^ with ISAAC shows two 
bright asymmetric arcs (A and B) close to the source, and two fainter trails about 10" and 
40" farther to th e SW. The H2 emission o verlaps well with the CO 3-2 blue-shifted outflow 
ejected by IRS 1 dvan Kempen et al.ll2009l ). The flow has a RA.~203°. 



MHO 2008 (A-C), located about 1' to the SW and well aligned with the previous flow 
(see Figure [27]), could be part of MHO 2003. On the other hand, MHO 2008 D is a bow-shock 
east- west oriented and should not be part of the flow (see Figure [30]) , but it is likely driven 
by the Class II source CrA41. 

MHO 2006 (A-C) is a trail of faint precessing knots (P.A.~206°), located about 15" to 
the SW of IRS 2 (see Figure [3l]), which is the driving source. MHO 2006 D (see Figure [27]) 
is a faint patch and could be part of the flow. 

About 6' south-westward of the Coronet, we detect MHO 2007, the molecular counter- 
part of HHlOl (see Figure [27]). The MHO extends for roughly 40", NNE-SSW elongated, 
with two faint arcs in front, likely a bow-shock structure and several bright scattered knots 
behind (A and B are particularly bright). The MHO is shown in Figure [32] where an ISAAC 
H2 continuum-subtracted image is shown. 

The last newly detected flow, MHO 2005 (knots A-I), is located ~4' to the west of the 
Coronet and extends for about 6' around S CrA (see Figure [33]) . roughly showing a broad 
S-shape. Knots C and E, the brightest knots of the flow, along with D, are the molecular 
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counterpart of HH 82 B and A, respectively. Knot F is positioned about 1' south and might 
not be part of the flow. The HH 7 29 knots, which should be driven by S Cr A (see e.g., 



Strom et al.lll974j : IWang et al.l 120041 ). are not detected in our H2 images. 



B.2.2. P.M.s 

Figure [35] shows the MHO 2004 A P.M., confirming that the bow-shock is not related to 
the main jet, but is moving WNW, possibly associated with the same fiow as MHO 2001 Q. 

Finally, our kinematical analysis reveals three or four different flows in the CrA south- 
western region, as shown in Figures [Ml EH and [381 MHO 2003 and knots A-C of MHO 2008 
are likely part of the same flow (Figure [36!) moving towards the SW. On the other hand, 
knot D stands alone moving towards the SSW (Figure [371) . likely driven by the nearby 
YSO candidate CrA-41, and MHO 2006 (same flgure) represents another flow. MHO 2007 
(Figure [55]) is roughly aligned with the MHO 2003-6 flow, with similar tangential velocities. 
However, it is not clear whether or not it could be part of it, because of the different 
orientation of the leading bow-shock, which rather points to the SSW. 



B.3. Notes on individual Spitzer flows 

The flrst outflow detected in our Spitzer maps, named Spitzer outflow 1, is located on 
the western edge of the CrA star forming region, about 17' from the Coronet. Figure [5^ 
shows the morphology of the outflow, which is composed of several knots and bow-shocks, 
roughly east-west oriented. These knots likely originate from a precessing jet, as indicated by 
their scattered placement. No obvious counterpart on the eastern edge of CrA is observed. 

Spitzer outflow 2, located ~ 1' west of IRS 11 (see Figure WJh . is a curved jet which in- 
cludes MHO 2005 A and B, but is not detected in our H2 SofI map, because it is slightly below 
the detection limit. It is worth noting that Spitzer outflow 1 and 2, along with MHO 2005, 
could be part of a larger east-west precessing flow (see Figure [4T|) . possibly driven by a YSO 
inside the Coronet. On the other hand, as seen in Figure WT\ the bulk of Spitzer outflow 1 is 
well aligned with IRS 2, and thus it is not clear which is the driving source. 

The IRAS 32 outflow, the third Spitzer flow discovered (P.A.~39°), originates from 
IRAS 32, and has a clear precessing morphology, as seen in Figure [121 The outflow extends 
for ~20.3', i.e. it has a projected length of about 0.77 pc. Knots from A to E are part of 
the north-eastern lobe, while knots from F to N, along with HH861A, B, and HH862 are 
located in the south-western lobe. The shock-excited nature of these knots is revealed by 
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our H2 image taken with ISAAC /VLT (lower left panel of Figure H2|) . centered around the 
YSO. Due to the small FoV, we could only identify the closest knots (i. e. D-G), which thus 
we also labeled as MHO 2010. Our H2 and images also show the two outflow cavities, 
illuminated by scattered light from the YSO. In between the two cavities, a dark band is 
well visible, possibly indicating the presence of a disk. 
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Fig. 1. — Coverage map for the Spitzer data taken in CrA. The thick black hne shows the 
IRAC coverage in bands 1 and 3 (3.6/5.8 /xm) and the thick black dotted line shows the 
MIPS 24 /im coverage. The underlying greyscale is a near-infrared extinction map created 
using 2MASS point sources. 
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Fig. 2. — Three-color image of the main region in CrA, including the Coronet (bright white 
saturated region toward the center), made from the Spitzer IRAC 4.5 fim (blue), 8.0 /im 
(green), and MIPS 24 (red) bands. 
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Fig. 3. — Three-color image of the streamer to the west of the main CrA region shown in 
Figure El made from the Spitzer IRAC 4.5 /xm (blue), 8.0 /im (green), and MIPS 24 /im 
(red) bands. 
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Fig. 4. — SofI (2007) H2 mosaic. Different polygons show regions mapped in different epochs 
(Yellow: 1999, Cyan: 2000, Green: 2003, Magenta: 2005). Known HH objects are indicated 
with green crosses. Newly detected H2 and Spitzer knots are also displayed as blue and red 
crosses, respectively. 
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Fig. 5. — Color-magnitude diagrams showing the color space where YSOs are selected from 
Spitzer data along with the f ull SWIRE and re-sa mpled SWIRE catalogs. Based on the c2d 
YSO identification scheme by lHarvey et al.l (120071 ). where red and black dot-dashed lines are 
hard and fuzzy identification limits, respectively. 
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Fig. 6. — SEDs are shown for the Class I and Flat spectrum sources in CrA selected from 
the Spitzer data as discussed in Section 13.11 Open circles represent the Spitzer and ancillary 
data. The number that appears in each panel corresponds to the object number in Table 
m Note that sources CrA-2 (a known galaxy, see discussion in §A.1.2p and CrA-32 (see 
discussion in §A.1.32p have been removed from our sample of YSO candidates. 
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Fig. 7. — SEDs are shown for the Class II sources in CrA selected from the Spitzer data as 
discussed in Section 13.11 Open circles represent the Spitzer and ancillary data, filled circles 
represent the dereddened fluxes, the grey line represents the stellar photospheric emission 
expected from the best-fit NextGe n model, and the dash ed black line represents the average 
SED for T Tauri stars in Taurus (iD'Alessio et al.lll999l ). The number that appears in each 
panel corresponds to the object number in Table HI Note that sources CrA-17 (see discussion 
in §A.1.17p . CrA-25 (see discussion in §A.1.25p . and CrA-39 (a known galaxy, see discussion 
in §A.1.39p have all been removed from our sample of YSO candidates. 
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Fig. 8. — SEDs are shown for the Class III sources in CrA selected from the Spitzer data as 
discussed in Section 13.11 Open circles represent the Spitzer and ancillary data, filled circles 
represent the dereddened fluxes, the grey line represents the stellar photospheric emission 
expected from the best-fit NextGe n model, and the dash ed black line represents the average 
SED for T Tauri stars in Taurus (ID'Alessio et al.lll999l ). The number that appears in each 
panel corresponds to the object number in Table |H The fit to CrA-11 is not very good in the 
visible wavelength range because it has a known spectral type of B8V (see §A.l.lip . whereas 
the model assumes a spectral type of K7 (which is a reasonable approximation for most of 
the sources). 
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Fig. 9. — A Ks vs. - [24 /iin] color- magnitude diagram for CrA. Symbols represent the 
following classifications: "Gale", or galaxies (green squares), YSO candidates that appear 
in Table H] (filled red circles), and YSO candidates that turned out to be known galaxies 
(filled red squares). The phase space used to select star+dust sources as candidate YSOs, as 
discussed in §3.2^ is shown by the dashed lines. Those star+dust sources that were rejected 
as candidate YSOs are represented by filled blue squares. The star+dust sources that met 
the criteria to be considered YSO candidates, and which appear in Table [5l are annotated 
and represented by filled blue circles. 
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Fig. 10. — All 116 YSO candidates are overlaid on a 24 /im image; the image inset is the 
Chandra X-ray image of the Coronet (the blue square near the center is RCrA). The Class I 
(+) and Flat spectrum (x) sources are shown in red, Class II (diamond) sources are shown in 
green, and Class III (square) sources are shown in blue. The contour shown is the extinction 
map made from 2MASS + Spitzer data, with levels Ay=[5, 10, 15, 20]. This extinction map 
has a resolution of 180". Most Class I/Flat sources are located in the high density cluster, 
while Class II and III sources are found spread further out. 



- 73 - 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



1 4 - CrA 



12 



10 



8 - 



2 - 



- 1 I - 







I 1 1 1 I I I I I I I I I I I I I 1 1 1 I I 1 1 I I 



-4 



-3 -2 -1 

log Luminosity (L0) 



Fig. 11. — Bolometric luminosity function for the YSO candidates in CrA selected using the 
c2d/GB method (discussed in § 13. H solid line), with the completeness correction included 
for comparison (dashed line), as discussed in § 13.51 In addition, the bolometric luminosity 
function for all YSO candidates is shown (dotted line). 
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Fig. 12. — The residual MST sub-networks for the YSO positions after those longer than 
Lcrit have been excluded for the case of all 116 YSO candidates (upper panel) and for the 
62 Class I, Flat spectrum, and Class II sources (lower panel). The branch lines connecting 
sources that are closer than the critical branch length (see Figure fT3ll are shown as grey lines. 
The convex hull (solid black line), which is discussed in § 14. 2[ is overlaid on both figures. 
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Fig. 13. — (Upper panel) MST critical branch length for all 116 YSO candidates. (Lower 
panel) MST critical branch length for the 62 Class I, Flat spectrum, and Class II sources 
(Class III sources are excluded). 
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R.A, [J2000] 



Fig. 14. — SMA observations of the RCrA region. The image is the Spitzer 4.5 fim image. 
Two SMA mosaics are shown. On the right, indicated by the large dashed black circle, 
are the SMA observations of IRS 5 at 226 GHz. On the left, indicated by the large dashed 
oval-like shape, are SMA observations of the IRS 7 region at 224 GHz. The circle and oval, 
respectively, indicate the 50% sensitivity levels of the observations. In the IRS 5 region, 
the red contours indicate emission at the 6a, 9a, 18a levels, where a = 3.7 mJy/beam 
(4'.'6 X 2'.'6). In the IRS 7 region, the green contours indicate emission at the 3a, 5a, and 7a 
levels, , while the yellow contours indicate emission at the 15a, 27a, and 39a levels, where 
a = 7.4 mJy/beam (5'.'5 x 2'.'3). Sources with black labels were detected by the SMA, 
while those with grey labels (in parentheses) we re not. R Cr A is detected at the 5a level. 
The sou rces IRS 7B and SM A 2 (Radio Source 9: iBrownl 119871 ) are equivalent to SMAl and 
SMA 2 (iGroppi et al.l 120071 . 271 GHz). The synthesised beams are shown as grey ovals at 
lower right (IRS 5) and lower left (IRS 7). Crosses indicate artifacts in the Spitzer data.. 
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Fig. 15. — SMA observations of the CrA IRAS 32 region. The location of YSO candidate 
CrA-33 is also shown, although it was not detected with the SMA. The image is the Spitzer 
4.5 /xm image. The blue and red contours represent blue- and red-shifted integrated ^^CO 2-1 
emission, respectively. The white contours represent 226 GHz continuum emission. Contour 
levels for the blue- and red-shifted emission are 3a, 6a, etc, where la = 1.0 Jy km/s and 
0.7 Jy km/s, respectively. The CO is integrated over LSR velocities of [0,4.2] km/s and 
[6.7,11.6] km/s for blue- and red-shifted emission, respectively. The contours of continuum 
emission are 3a, 6a, etc, where la is 4.4 mJy/beam. The large dashed circle shows the 
primary beam size of the SMA at 226 GHz, while the small grey oval at lower left shows the 
synthesised beam of 3'.'2 x 2'.'2. 
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Fig. 16. — Spitzer three-color image of the CrA star- forming region, made from the IRAC 
3.6 /um (blue), 4.5 //m (green), and 8 /im (red) images. The green polygon shows the area 
covered by our SofI H2 map. Green, blue, and red crosses indicate the position of HH 
objects, H2, and Spitzer knots, respectively. 
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Fig. 17. — H2 flow chart of the CrA star-forming region. Proper motions in 500 yrs are 
indicated by arrows. 
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Fig. 18. — SofI B.2 image of the north-eastern flows detected outside the Coronet. 
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Fig. 19.— Close up view of MHOs 2009, 2011, and 2012 at 2.12 ^m, from a Sofl/NTT 
image (see Figure [T8|) . Labels and crosses indicate names and positions of detected knots 
and sub-structures along the flows. 




Fig. 20.— Close up view of MHOs 2000 and 2013 at 2.12 /zm, from a ISAAC/VLT image 
(see Figure [T8|) . Labels and crosses indicate names and positions of detected knots and 
sub-structures along the flows. 
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Fig. 21. — Close up view of MHOs 2014 at 2.12 /iin, from a Sofl/NTT image (see Figure [TSl) . 
Labels and crosses indicate names and positions of detected knots and sub-structures along 
the flows. 
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Fig. 22. — (Upper panel) SofI H2 continuum-subtracted image of MHO 2001 and 2002 in 
the Coronet region. Labels and crosses indicate names and positions of detected knots and 
sub-structures along the flows. (Lower panel) Close up view of the region at 2.12 fim, from 
a ISAAC/VLT image in logarithmic scale. 



- 84 - 




Fig. 23. — H2 flow chart of MHO 2001 and 2002. Proper motions in 100 yrs and their error 
bars are indicated by arrows and elhpses, respectively. 
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Fig. 24.— H2 flow chart of MHO 2009, 2011, 2012 and part of MHO 2001. Proper motions 
in 100 yrs and their error bars are indicated by arrows and elhpses, respectively. 
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Fig. 25. — H2 flow chart of MHO 2000 and 2013. Proper motions in 100 yrs and their error 
bars are indicated by arrows and elhpses, respectively. 




Fig. 26. — H2 flow chart of MHO 2014. Proper motions in 100 yrs and their error bars are 
indicated by arrows and elhpses, respectively. 
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R.A. CJ2000) 

Fig. 27. — SofI H2 image of the south-western flows detected outside the Coronet. MHO 
2001 knots inside the Coronet are labeled with blue crosses. 




Fig. 28. — ISAAC H2 continuum-subtracted image of MHO 2004 (see Figure 1271) . Labels and 
crosses indicate names and positions of detected knots and sub-structures along the flows. 
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Fig. 29. — ISAAC H2 continuum-subtracted image of MHO 2003 (see Figure [27]). Labels and 
crosses indicate names and positions of detected knots and sub-structures along the flows. 
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Fig. 30. — ISAAC B.2 continuum-subtracted image of MHO 2008 (see Figure [271). Labels and 
crosses indicate names and positions of detected knots and sub-structures along the flows. 
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Fig. 31. — ISAAC B.2 continuum-subtracted image of MHO 2006 (see Figure [271). Labels and 
crosses indicate names and positions of detected knots and sub-structures along the flows. 
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Fig. 32. — ISAAC H2 continuum-subtracted image of MHO 2007, the molecular counterpart 
of HH 101 (see Figure [27]). Labels and crosses indicate names and positions of detected knots 
and sub-structures along the flows. 
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Fig. 33. — SofI H2 image of the western flows detected outside the Coronet. 
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Fig. 34. — SofI H2 continuum-subtracted image of MHO 2005 (see Figure [33 




Fig. 35. — H2 flow chart of MHO 2004. Proper motions in 100 yrs and their error bars are 
indicated by arrows and elhpses, respectively. 
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Fig. 36.— H2 flow chart of MHO 2003 and part of MHO 2008. Proper motions in 100 yrs 
and tlieir error bars are indicated by arrows and ellipses, respectively. 
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Fig. 37.— H2 flow chart of MHO 2006 and part of MHO 2008. Proper motions in 100 yrs 
and their error bars are indicated by arrows and elhpses, respectively. 
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Fig. 38. — H2 flow chart of MHO 2007. Proper motions in 100 yrs and their error bars are 
indicated by arrows and eUipses, respectively. 
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Fig. 39. — Spitzer three-color image of Spitzer outflow 1, made from the IRAC 3.6 /im (blue), 
4.5 //m (green), and 8 /im (red) images. Labels indicate the position of the newly detected 
knots (EGOs). 
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Fig. 40. — Spitzer three-color image of Spitzer outflow 2, made from the IRAC 3.6 //m (blue), 
4.5 //m (green), and 8 /xm (red) images. Labels indicate the position of the newly detected 
knots (EGOs) and MHO 2005 A and B. 
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Fig. 41. — Spitzer three-color image of the region around Spitzer outflow 1, 2 and MHO 2005. 
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Fig. 42. — Spitzer three-color image of the IRAS 32 (IRAS 18595-3712) outflow, made from 
the IRAC 3.6 jjm (blue), 4.5 /im (green), and 8 /xm (red) images. Labels indicate the position 
of the newly detected knots (EGOs), MHO 2010 and previously known HH objects. Lower 
left panel displays a close up view of the inner region of the outflow at 2.12 fim, showing the 
inner knots, the outflow cavities, and a possible disk. 
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Table 1. Spitzer Observations of CrA 



Instrument 


AOR ID 


PID 


Observation Date 


IRAC 


0003650816 


6 


2004-04-20 




0017672960 


30784 


2006-09-25 




0017673472 


30784 


2006-09-25 




0027041280 


30574 


2008-05-10 


MIPS 


0003664640 


6 


2004-04-11 




0017673216 


30784 


2007-05-30 




0017673728 


30784 


2007-05-29 




0027042816 


30574 


2008-10-23 



Table 2. Detection Statistics (S/N>3) in CrA 



Detected with/in... Number of Detections 

■IRAC Band 1 73254 

.IRAC Band 2 50096 

.IRAC Band 3 11023 

.IRAC Band 4 6381 

.MIPS Band 1 1919 

.MIPS Band 2 65 

.MIPS Band 1 and Band 2 54 

.2MASS (total) 8974 

.all 4 IRAC bands, but not 2MASS 350 

.2MASS alone 4 

.IRAC, but not 2MASS 51745 

.MIPS Band 1 and 2MASS Kg 557 



Table 3. Journal of Observations - Ancillary Near-Infrared Data 



Date of Obs. 


Telescope/ 


Filter 


Resolution 


Seeing 


Exp. Time 


FoV 


Centered on/ 


(d,m,y) 


Instrument 


Band 


(" pixcl-i) 


(") 




('X') 


Target 


2007 08 23 


ESO-NTT/SofI 


Ha, K, 


0.29 


1.3 


1440, 290 


10x10 


HH 734 


2007 08 23 


ESO-NTT/SofI 


Ha, K, 


0.29 


1.3 


1440, 290 


10x10 


HH 730 


2007 08 23 


ESO-NTT/SofI 


Ha, K, 


0.29 


1.3 


1440, 290 


10x10 


HH 101 


2007 06 26 


ESO-NTT/SofI 


Ha, Ks 


0.29 


1.1 


1440, 290 


10x10 


HH 99 


2005 06 17 


ESO-VLT/ISAAC 


Ha, K, 


0.15 


0.8 


500, 45 


5x2.5 


IRAS 32 


2005 06 16 


ESO-VLT/ISAAC 


Ha 


0.15 


0.8 


4800 


5x7 


R CrA 


2003 07 12 


ESO-VLT/ISAAC 


Ha, NB 2.09 


0.15 


0.9 


1440, 290 


6.3x6.3 


HH 101 


2000 07 19 


ESO-VLT/ISAAC 


Ha 


0.15 


0.9 


1200 


6.8x6 


HH 99 


1999 06 06 


ESO-NTT/SofI 


Ha, K, 


0.29 


0.8 


600, 100 


5.2x5.2 


R CrA^' 



^Data already published in lCaratti o Garatti et aL I 1I2OO6I ') 



CO 



Table 4. Sources Classified as YSO Candidates in CrA based on Spitzer IRAC and MIPS 



ID 




RA 




Dec 






3 . G fiTci 


4 5 fjm 


5. 8 fj,m 


8.0 ;im 


24 fim 


70 fJ.ni 




Other Names'* 




(J2000) 


(J2000) 




[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


(Class) 




CrA-1 


18 


56 


39 


76 


-37 


07 


20 


8 


47.4±2.4 


38.8±2.0 


37.7±1.8 


42.6±2.1 


35.6±3.3 




-1.17 (II) 




CrA-2 


18 


57 


52 


63 


-37 


14 


40.1 


2.04±0.12 


8.90±0.44 


35.8±1.8 


93.6±4.6 


171. ±16. 


1690±196 


1.73 (G) 


Leda Galaxy 90315 


CrA-3 


18 


59 


43 


92 


-37 


04 


0.11 


1.69±0.09 


2.22±0.11 


2.93±0.20 


3.58±0.19 


6.97±0.66 




-0.18 (F) 


ISO-CrA 55^ 


CrA-4 


18 


59 


50 


95 


-37 


06 


31 


6 


4.69±0.25 


3.93±0.20 


3.16±0.17 


2.79±0.19 


1.93±0.22 




-1.53 (II) 


DENIS-J185950. 9-370632 


CrA-5 


19 


00 


15 


55 


-36 


57 


57 


7 


1.07±0.06 


2.84±0.15 


5.79±0.3a 


8.04±0.40 


13.3±1.2 




0.51 (I) 


ISO-CrA 76' 


CrA-6 


19 


00 


29 


07 


-36 


56 


03 


8 


51.6±2.7 


36.0±1.8 


26.1±1.3 


18.5±0.9 


36.0±3.3 


52.0±6.4 


-1.57 (II) 


CrAPMS 8^; GP g2^ 


CrA-7 


19 


00 


38 


94 


-36 


58 


14 


7 


40.8±2.2 


33.8±1.7 


25.8±1.3 


17.8±0.9 


4.49±0.43 




-2.08 (III) 


ISO-CrA 93^ 


CrA-8 


19 


00 


45 


31 


-37 


11 


48 


2 


7.73±0.40 


6.66±0.33 


5.73±0.29 


6.85±0.34 


14.6±1.4 




-0.90 (II) 


CrA-444* 


CrA-9 


19 


00 


58 


05 


-36 


45 


05 





56.3±3.3 


43.8±2.3 


31.7±1.5 


25.1±1.2 


179. ±17. 


154. ±20.'' 


-1.04 (II) 




CrA- 10 


19 


00 


59 


75 


-36 


47 


11 


2 


4.05±0.22 


3.95±0.20 


3.46±0.18 


3.87±0.20 


4.36±0.41 


35.3±6.0 


-1.03 (II) 


CrA-432* 


CrA- 11 


19 


01 


03 


26 


-37 


03 


39 


4 


473. ±29. 


335. ±18. 


243. ±12. 


174. ±9. 


225. ±21. 




-1.75 (III) 


HD 176269; HR 7169^ 


CrA-12 


19 


01 


16 


29 


-36 


56 


28 


3 


8.69±0.45 


7.81±0.38 


6.77±0.33 


6.08±0.29 


7.46±0.70 




-1.23 (II) 


V 667 CrA; CrA-4110'* 


CrA-13 


19 


01 


18 


95 


-36 


58 


28 


4 


36.5±1.9 


34.8±1.8 


36.2±1.7 


35.0±1.7 


52.8±4.9 


54.6±7.2 


-0.92 (II) 


CrA-466*; ISO-CrA 127' 


CrA- 14 


19 


01 


29 


01 


-37 


01 


48 


4 


19.5±1.0 


13.8±0.7 


9.62±0.48 


6.61±0.32 


2.93±0.29 




-2.12 (III) 


HH 101 IRS l""; G-94^ 


CrA-15 


19 


01 


32 


31 


-36 


58 


03 


Q 


21.2±1.1 


18.9±0.9 


16.4±0.8 


13.8±0.7 


16.2±1.5 




-1.19 (II) 


IRS 14^; TS 2.9*; G-87'^ 


CrA- 16 


19 


01 


33 


85 


-36 


57 


44 


g 


68.3±5.1 


70.9±3.7 


66.6±3.2 


71.9±3.6 


131. ±12. 


173. ±30. 


-0.65 (II) 


IRS 13*; TS 2.8*; G-85'' 


CrA-17 


19 


01 


36 


26 


-36 


58 


03 





0.816±0.047 


0.678±0.041 


0.513±0.128 


0.340±0.a60 


5.73±0.92'^ 




_ 


removed from list 


CrA-18 


19 


01 


40 


41 


-36 


51 


42 


3 


49.0±2.7 


42.0±2.1 


41.3±2.1 


46.7±2.8 


118. ±14. 




-0.70 (II) 


G-65^ 


CrA-19 


19 


01 


48 


03 


-36 


57 


22 


2 


389. ±24. 


798. ±50. 


1270. ±76. 


1810. ±107. 


4500. ±570.'' 




0.78 (I) 


IRS 5*; TS 2.4*; MMS 12^; SMM 4'" 


CrA-20 


19 


01 


48 


46 


-36 


57 


14 


7 


8.56±0.51 


26.2±1.3 


52.8±3.5 


80.9±4.2 


1110. ±155. 




1.41 (I) 


IRS 5N" 


CrA-21 


19 


01 


51 


12 


-36 


54 


12 


4 


30.3±1.5 


32.6±1.7 


37.2±1.8 


45.2±2.2 


39.9±3.7 




-0.68 (II) 


IRS 8*; TS 2.2* 


CrA-22 


19 


01 


51 


86 


-37 


10 


44 


7 


3.75±0.19 


3.32±0.17 


2.79±0.15 


2.27±0.13 


1.79±0.21 




-1.47 (II) 




CrA-23 


19 


01 


53 


75 


-37 


00 


33 


9 


3.16±0.17 


2.96±0.15 


2.74±0.17 


3.01±0.16 


3.09±0.34 


113. ±40.'' 


-1.09 (II) 


Star A'2 


CrA-24 


19 


01 


55 


60 


-36 


56 


51 


1 


11.0±1.6 


30.3±3.2 


54.0±5.2 


93.1±13.5'' 


297.±50.'= 




0.60 (I) 




CrA-25 


19 


01 


57 


27 


-36 


52 


05 





0.588±0.179'= 


1.30±0.3S'= 


1.43±0.21 


1.06±0.32'^ 


3.22±0.47'' 




HH 


removed from list 


CrA-26 


19 


02 


06 


80 


-36 


58 


41 





2.29±0.12 


2.22±0.11 


2.35±0.14 


2.48±0.14 


3.11±0.64'' 




-0.84 (II) 




CrA-27 


19 


02 


10 


44 


-36 


53 


44 


5 


6.94±0.36 


4.55±0.23 


3.55±0.20 


2.12±0.11 


0.962±0.177 




-2.24 (III) 




CrA-28 


19 


02 


12 


00 





03 


09 


4 


6.15±0.31 


5.04±0.25 


4.26±0.22 


4.62±0.24 


4.63±0.45 




-1.27 (II) 


ISO-CrA 143'; G-14'' 


CrA-29 


19 


02 


14 


63 


-37 


00 


32 


9 


9.21±0.49 


11.4±0.57 


12.2±0.59 


11.8±0.6 


14.2±1.3 




-0.64 (II) 


ISO-CrA 145' 


CrA-30 


19 


02 


27 


07 


-36 


58 


13 


1 


300. ±17. 


242. ±14. 


256. ±13. 


248. ±13. 


150. ±14. 


• • ■ 


-1.37 (II) 


Ha 14'^; HBC 6S0'* 


CrA-31 


19 


02 


33 


07 


-36 


58 


21 


2 


356. ±19. 


335. ±20. 


347. ±17. 


339. ±17. 


238. ±22. 


167. ±18. 


— 1.09 (II) 


Ha 3'^; ISO-CrA 159' 


CrA-32 


19 


02 


56 


82 


-37 


07 


19 


4 


0.314±0.018 


0.374±0.021 


0.351±0.039 


0.234±0.038 


5.21±0.92'' 






removed from list 


CrA-33 


19 


03 


01 


03 


-37 


07 


53 


4 


7.04±0.36 


10.2±0.5 


15.3±0.8 


16.9±0.8 


32.9±3.1 




0.06 (F) 


IRAS 32d'^ 


CrA-34 


19 


03 


09 


16 


-36 


57 


22 





149. ±8. 


121. ±6. 


92.2±4.4 


68.2±3.4 


20.4±1.9 




-2.02 (III) 


ISO-CrA 198' 


CrA-35 


19 


03 


11 


84 


-37 


09 


02 


1 


19.6±1.0 


14.9±0.7 


11.7±0.6 


10.4±0.5 


11.7±1.1 




-1.46 (II) 


ISO-CrA 201' 


CrA-36 


19 


03 


24 


29 


-37 


15 


07 


7 


44.9±2.5 


44.9±2.3 


52.1±2.7 


61.3±3.0 


132. ±12. 


96.1±9.7 


-0.39 (II) 




CrA-37 


19 


03 


55 


24 


-37 


09 


35 


9 


1.43±0.07 


1.95±0.10 


2.63±0.14 


3.84±0.20 


20.8±1.9 


78.1±8.3'' 


0.36 (I) 




CrA-38 


18 


58 


56 


40 


-37 


07 


37 


9 


519. ±28. 


298. ±17. 


264.±13. 


176. ±9. 


51.3±4.7 




-2.21 (III) 


ISO-CrA 13' 


CrA-39 


19 


00 


23 


47 


-37 


12 


24 


2 


3.81±0.33 


2.21±0.18 


1.65±0.15 


1.96±0.27 


6.92±0.72 


943. ±97.'' 


-0.95 (G) 


6dFGSgJ190023.5-37 


CrA-40 


19 


01 


25 


75 


-36 


59 


19 


1 


259. ±18. 


223. ±12. 


216. ±11. 


170. ±8. 


103. ±10. 




-1.38 (II) 


VSSt 18"^ 


CrA-41 


19 


01 


41 


62 


-36 


59 


52 


7 


298. ±17. 


295. ±16. 


293. ±18. 


324.±16. 


293. ±27. 




-1.04 (II) 


Ha 2'*; HBC 677'*" 


CrA-42 


19 


01 


50 


48 


-36 


56 


38 


4 


106. ±7. 


143. ±9. 


179. ±10. 


196. ±12. 


348. ±35. 




-0.39 (II) 


IRS 6* 


CrA-43 


19 


01 


58 


54 


-36 


57 


08 


5 


22.5±1.3 


72.1±3.8 


136. ±8. 


195. ±11. 


807.±78. 




0.88 (I) 


SMM2"'; CHLT 15''^ 


CrA-44 


19 


02 


58 


67 


-37 


07 


35 


9 


6.30±0.71 


17.2±1.2 


17.3±1.1 


12.9±a.7 


2040. ±193. 


28500. ±2870 


1.66 (I) 


IRAS 32c'5 


CrA-45 


19 


03 


16 


09 


-37 


14 


08 


2 


223. ±12. 


224. ±12. 


256. ±13. 


278. ±14. 


676. ±63. 


680. ±69. 


-0.42 (II) 


VSSt lO'"; IRAS 34'^ 



"References are as follows: lOlofsson et al.l i 19991 1; HWalter et al.l jl997D ; ' iGlass fc PenstoiJ jl975D ; iLopez Marti et al.l i2005l ); iKnacke et al.l i 19731 ); " ^Reipurth fc WamstekeJ il983l) ; 



' iSicilia-Aguilar et all i200gD: IXaylor &c StorevI JlQS^T: ^Chini et all 120031): ^^Nutter et al.| {200^ ; ^ ^Forbrich Preibis"^ j2007l ): ^ ^Grahaml j 19931) : ^ iMarraco &: RvdgreU il98lD : 
^ iHerbig &c Belli jl98a ) : ^ Wilking et al.Nl992| ): ^ Wrba et al.| il976a| ):"^ ' iChoi et al.N2008| ). 

^The flux for CrA-9 was computed using aperture photometry from the 70 fJ.ai image; the catalog entry for this source shows it as a non-dctcction. See Scction |A.1.9] 

"^The flux for this source at this wavelength has been "band-filled" (described in ■ 

'^This 70 i-im flux is listed in the catalog with a "P" flag, meaning that there were multiple catalog detections near this 70 fim source. The catalog source chosen to match the 70 fim detection 
gave the most "consistent" SED, and was also the one closest to the 70 i-im detection (see the c2d Final Delivery Document, [E'vans et al.||2007| ). 
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Table 5. Sources Classified as YSO Candidates in CrA based on 2MASS and MIPS 24 /im 



ID 




RA'» 




Doc 






3.6 fJ,m 


4.5 fj,m 


5.8 


8.0 


24 


70 p,m 






Other Namcs^ 




(J2000) 




(J2000) 




[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


(Class) 




CrA-46 


18 


55 


56 


32 


-37 


00 


07 


1 










15.6±1.5 




-0.11 


(F) 




CrA-47 


18 


56 


40 


28 


-36 


55 


20 


8 










4.46±0.43 




-1.04 


(II) 




CrA-48 


18 


57 


07 


86 


-36 


54 


04 


4 










2.53±0.26 




-1.21 


(11) 




CrA-49 


18 


57 


22 


48 


-37 


34 


42 


5 










3.0D±0.30 




-0.96 


(11) 




CrA-50 


18 


57 


52 


22 


-37 


14 


45 


8 


0.883±0.053 


0.567±0.031 


0.403±0.060 


0.286±0.098'^ 


9.87±1.87 




-0.66 


(G) 


Leda Galaxy 90315 


CrA-51 


18 


59 


34 


26 


-37 


21 


41 


5 


1.26±0.11 




0.848±0.083 




5.69±0.55 




-0.49 


(11) 




CrA-52 


19 


00 


01 


58 


-36 


37 


06 


2 










124.0±11.5 




-0.95 


(11) 


CrAPMS 7I; HBC 674^; Ha 6^ 


CrA-53 


19 


01 


11 


49 


-36 


45 


33 


8 


5.99±0.31 




3.49±0.19 




4.35±0.42 




-1.31 


(11) 


LMlO-1* 


CrA-54 


19 


01 


39 


15 


-36 


53 


29 


4 


93.4±4.9 


59.4±3.0 


67.4±23.9'' 


131.0±48.4'^ 


602. ±101. 




-1.01 


(11) 


HD 176386B 


CrA-55 


19 


01 


55 


23 


-37 


23 


41 













1420. ±134. 


1540. ±161. 


-0.51 


(11) 


DG CrA; HBC 289^; IRAS 24^; 
IRAS 18585-3728 


CrA-56 


19 


02 


16 


66 


-36 


45 


49 


4 


15.2±0.8 




7.3±0.4 




10.4±1.0 




-1.39 


(11) 


ISO-CrA 146''; CrA-4109''; 
DENIS-P J190216. 6-364549 


CrA-57 


19 


03 


25 


48 


-36 


55 


05 


3 


12.4±0.6 




5.74±0.29 




5.43±0.52 




-1.60 


(11) 


ISO-CrA 217" 



^Coordinates from our Spitzer catalog. 

''References are as follo ws: ^Walter etlj] )1997D ; ^Herbig fc Belli 119861) ; ^Marraco Sz RvdgrenI il981l ); lL6pez Marti et al.l I201CI) ; IWilking et al.l i 19921 ); lOlofsson et al.l i 19991 ); 
iLopez Marti et al.l i2005l ). 

''This flux has a quality of detection of "D" meaning it has a 2 < S/N < 3. Fluxes that do not have this note all have a quality of detection of "A" (S/N >7), "B" (S/N >5), or "C" 
(S/N > 3). 



Table 6. Class III YSO Candidates from Spitzer + ROSAT 



Spitzer ID 


ROSAT ID-'' 


Offset'' 




RA"^ 




Dec 






3.6 


4.5 f-Ln\ 


5.8 fim 


8.0 fj.ia 


24 /im 


ct 




Other Names^ 




RX # 






(J2000) 




{J2000) 




[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 








CrA-58 


J1855. 3-3730 


6 


26 


18 


55 


18 


91 


-37 


29 


56 


7 


1350. ±77. 


804. ±50. 


610. ±34. 


338. ±18. 


34.9±3.2 


-2 


86 


HD 175073; LTT 7499 




J1857.5-3732 


2 


54 


18 


57 


34 


19 


-37 


32 


33 


3 










0.826±0.144 


-2 


76 




CrA-60 


J1857.7-3719 


13.1 


18 


57 


45 


15 


— 37 


19 


37 


9 


32. Oil. 6 


22.7±1.1 


15.7±0.8 


9.00±0.43 


0.882±0.150 


-2 


68 


[P98c] ROle^ 


CrA-61 


J185802. 3-365348 


6 


29 


18 


58 


01 


82 


-36 


53 


45 


4 










1.87±0.21 


-2 


82 


CrAPMS 5^ 


CrA-62' 


J1858. 7-3706 


6 


24 


18 


58 


43 


33 


-37 


06 


27 


3 


5980. ±690. « 


2370. ±286. 


3120. ±158. 


1960. ±111. 


228. ±21. 


-2 


67 


e CrA; HR 7152; R04'- 




J lOOOOO. Z-oD41:UOji 


4 


35 


18 


58 


56 




-36 


40 


29 












n Q'/R-t-n on"/ 
u. y i Ditu. 1 






[i yocj txi 1 c 


CrA-64 


J1859. 2-3711 


1 


04 


18 


59 


14 


82 


-37 


11 


32 


6 


21.6±1.3 


14.6±0.9 


9.15±0.52 


5.90±0.32 


0.979±0.138 


-2 


74 


RO2I; CrAPMS 6^; Ha 11^ 


CrA-65 


J190039. 0-364805 


8 


17 


19 


00 


39 


30 


-36 


48 


12 


3 


27.2±1.4 


18.0±0.9 


12.4±0.6 


7.45±0.37 


0.898±0.146 


-2 


64 


RI4I; CrAPMS 9^ 


CrA-66 


J190109. 4-364751 


4 


00 


19 


01 


09 


70 


-36 


47 


52 


9 


118. ±7. 




57.1±2.8 




3.58±0.34 


-2 


79 


VSS VIII-27* 


CrA-67 


J1901. 4-3704 


8 


60 


19 


01 


25 


63 


-37 


04 


53 


7 


11.5±0.6 


8.56±0.42 


5.94±0.29 


3.70±0.18 


0.371±0.163''J 


-2 


40 


ISO-CrA 133^ 


CrA-68 


J190126. 8-365911 


4 


88 


19 


01 


27 


15 


-36 


59 


08 


5 


37.1±2.0 


26.2±1.3 


18.4±0.9 


11.0±0.5 


1.12±0.18' 


-2 


62 


R07al; ISO-CrA 135^ 


CrA-69 


J190128. 6-365933 


1 


82 


19 


01 


28 


71 


-36 


59 


31 


8 


95.1±5.5 


66.1±3.4 


47.7±2.3 


28.9±1.4 


2.86±0.29 


-2 


66 


ISO-CrA 136^; CXO 12'' 


CrA-70 


J190134.6-370058 


3 


27 


19 


01 


34 


84 


-37 


00 


56 


5 


258. ±15. 


152. ±8. 


111. ±6. 


67.5±3.4 


7.87±0.74 


-2 


74 


CrAPMS 1^; VSSt 13^ 


CrA-71 


J190139. 0-370207 


3 


88 


19 


01 


39 


32 


-37 


02 


07 


4 


25.9±1.4 


16.8±0.9 


11.5±0.6 


6.65±0.33 


0.718±0.193' 


-2 


64 


[GP75] R CrA 0** 


CrA-72 


J190140. 3-364429 


4 


20 


19 


01 


40 


55 


-36 


44 


32 





91.9±5.0 




40.9±2.1 




2.92±0.30 


-2 


72 


VSS VIII-26* 


CrA-73 


J190200. 1-370227 


4 


93 


19 


02 


00 


11 


-37 


02 


22 


1 


13.7±0.7 


9.21±0.46 


6.38±0.34 


3.87±0.22 


0.491±0.250'''= 


-2 


53 


Ha 13^ 


CrA-74 


J190201. 7-370746 


4 


02 


19 


02 


01 


96 


-37 


07 


43 


5 


129. ±7. 


76.6±4.3 


56.2±2.7 


31.6±1.6 


3.88±0.38 


-2 


79 


V702 CrA; CrAPMS 2^ 


CrA-75 


J190222. 0-365541 


1 


41 


19 


02 


22 


12 


-36 


55 


40 


9 


65.4±3.4 


43.2±2.2 


30.5±1.5 


19.0±0.9 


2.22±0.25 


-2 


66 


VSSt 24'^; HBC 679" 


CrA-76 


J190249. 7-364618 


5 


79 


19 


02 


49 


80 


-36 


46 


23 


7 


37.0±2.0 




15.7±0.8 




1.07±0.17 


-2 


69 


ISO-CrA 174^ 


CrA-77 


.1190408.2-370249 


3 


20 


19 


04 


07 


95 


-37 


02 


47 


8 


36.7±2.0 




16.7±0.8 




1.04±0.15 


-2 


76 


VSS VIII-21* 



^From Patten ct al. (private communication). The shorter coordinate names are sources from the PSPC (Position Sensitive Proportional Counter) catalog (they are less precise than the 
longer coordinate names). All the ROSAT IDs have a prefix of "RX" 

'^Offset is calculated as the distance between our Spitzer position and the ROSAT position from Patten et al. (priv. comm). Note that the ROSAT positions are computed from centroiding 
using the PROS package in IRAF. 

'^Coordinates from our Spitzer catalog. 



'^References arc as follows: ^Patten) |l998ll: ^Walter et ahl Il997f): "iMarraco &: Rydgrenl Il98lll; ' Vrba et all jl976bl h lOlofsson et al.l <199gl h * lForbrich &z PreibischI <200'1 . CXO is the listing 
number from Table 2 their paper); IVrba et al.| jl976a| "): iGlass &: PenstonI jl975l ); iHerbig 6*6111 ]l98g| ) 



^In Table 13 of lFernandez et al.| j 20081 ) . CrA-59 (RX J1857. 5-3732) is referred to as a binary, where the eastern component has a spectral type of M5 and the western component has a spectral 
type of M6. Only one source is detected in our Spitzer observations. Note that this source does not have any IRAC fluxes because it is off the field in the IR,AC images (the MIPS field extends 
a little further). 

iPattenI j 19981 ") says this is not a member because it is an eclipsing binary and if you work out the numbers, it is likely foreground. This could be true because these *are* very bright fluxes.] 

^This flux is flagged "K" in the catalog meaning that it is complex, i.e. there are two or more detections within 2 in the same band in the same epoch. The quality is likely still quite good, 
however. 

^■^There is also a PSPC coordinate name for this source, RX J1858. 9-3640, and although it is about 42" away, it is the same source due to imprecise PSPC coordinates. 
^The flux for this source at this wavelength has been "band-filled" (described in J3). 

■'This flux has a quality of detection of "D" meaning it has a 2 < S/N < 3. Fluxes that do not have this note all have a quality of detection of "A" (S/N >7), "B" (S/N >5), or "C" (S/N > 3). 
^This flux is flagged "U" in the catalog indicating that it is an upper limit because it has been band-filled and has a S/N < 5. 



Table 7. Spitzer Matches to Known YSOs and Candidate YSOs in CrA 



Name 




ra'' 




Dec 


b 




3.6 


4.5 


5.8 fj.m 


8.0 fj.m 


24 fj.m 


a 




cza Oatalog 


Other Names'' 




(J2000) 


{J2000) 




[nrJy] 


[mJy] 


[irijy] 


[injy] 


[mJy] 


(Class) 


Designation 




R CrA 


19 


01 


53. 


67 


-36 


57 


08 







,d,e 


10500. ±1980 


,d,e 


76.1±644.0'*' = 


-2.55 


(III)* 


star 


TS 2.10^; VSSt 1^ 


S CrA 


19 


01 


08. 


,62 


-36 


57 


20, 


3 


2580. ±187. 


3040. ±218. 


3510. ±224. 


2640. ±234. 


<195'*'*^ 


-0.80 


(11) 


star+dust(IRl) 


VSSt 3^ ; SMM 7^ 


T CrA 


19 


01 


58. 


78 


-36 


57 


49 


9 


1930. ±131. 


2110. ±127. 


2360. ±130. 


3520. ±202. 




-0.35 


(11) 




TS 4.3^, VSSt 2^ 


VV CrA 


19 


03 


06. 


,80 


-37 


12 


49, 


1 


1570. ±203. 


2900. ±479. 


12200. ±915. 


6770. ±678. 


<5800'*-'^ 


0.36 


(I) 


rising 


VSSt 4^ ; SMM 9^ 


TY CrA 


19 


01 


40. 


81 


-36 


52 


33, 


8 


844. ±82. 


631. ±59. 


1060. ±114. 


2900. ±372. 




-0.76 


(II) 


PAH-em 


VSSt 5^ 


HD 176270 


19 


01 


04. 


30 


-37 


03 


41, 


7 


684. ±42. 


468. ±25. 


359. ±18. 


190. ±9. 


21.7±2.1 


-2.81 


(III) 




HR 7170"^ 


IRS 1^ 


19 


01 


50. 


68 


-36 


58 


09, 


7 


2960. ±352. 


4520. ±526. 


14000. ±1150. 


<150'*'^ 


22300. ±4600.'* 


0.92 


(I) 




VSSt 15^; TS 2.6^ 


IRS 2 


19 


01 


41. 


56 


-36 


58 


31, 


2 


1770. ±235. 


2330. ±239. 


6580. ±375. 


7940. ±506. 


<1450'*''' 


0.75 


(I) 


risin^ 


TS 13.1-'- ; SMM 5"^ 


IRS 7A ( — 7W) 


19 


01 


55. 


32 


-36 


57 


21, 


9 


50.9±4.5 


239. ±24. 


546. ±43. 


1050. ±78. 


<2790'*''^ 


2.64 


(I) 


risin 




IRS 7B (=7E) 


19 


01 


56. 


40 


—36 


57 


28, 


1 


36.7±2.6 


257. ±18. 


697. ±38. 


957. ±51. 




2.78 


(I) 


rising 




IRS 10 


19 


02 


04. 


09 


— 36 


57 


01, 


2 


5.28±0.31 


4.49±0.24 


3.70±0.37 


2.66±0.16 


1.27±0.83'*''' 


— 1.48 


(II) 


star 


TS 4.2^; [HHD2008]A 569^ 


IRS 11 


19 


01 


41. 


89 


— 36 


56 


58, 


4 


3.77±0.20 


3.02±0.15 


2.32±0.17 


1.47±0.13 


<1.9'*''' 


— 1.77 


(III) 


star 


TS 2.5^; [HHD2008]A 367^ 


IRS 15 


19 


01 


47. 


89 


— 36 


59 


30, 


1 


7.76±0.39 


5.29±0.27 


4.66±0.24 


2.65±0.14 


0.0966±0.5590'*''' 


— 2.01 


(III) 


star 


TS 13.4^; GP n^ 


VSSt 8 


19 


01 


37. 


27 


— 36 


49 


21, 





206. ±11. 


131. ±7. 


100. ±5. 


62.3±3.2 


6.45±0.61 


— 2.67 


(III) 


star 


TS 1.8^ GP 


Ha 10 


18 


58 


51. 


86 


— 37 


19 


22, 


7 


7.80±0.39 


4.70±0.23 


3.18±0.17 


2.04±0.11 


0.132±0.136'*' = 


— 2.68 


(III) 


star 


HBC 673*^ 


Ha 15 


19 


04 


17. 


25 


— 36 


59 


03, 


2 


20.5±1.1 




13.4±0.7 




17.3±1.6 


— 1 .24 


(II) 


star+dust(IR,3) 




ISO-CrA 46 


18 


59 


31. 


79 


— 36 


50 


26, 


,3 


39.1±2.0 


24.6±1.2 


17.7±0.9 


10.5±0.5 


1.11±0.17 




(III) 


star 




ISO-CrA 53 
















4 


18.9±0.9 


11.5±0.6 


7.85±0.39 


4.79±0.24 


0.399±0.154'*''' 


~2 73 


(iii) 


star 


[P98e] R13a^ | 


ISO-CrA 140 


19 


02 


10 


03 


~36 


56 


49 


,0 


1.32±0.08 


1.40±0.09 


1.35±0.31 


0.662±0.097=' 


0.375±0.290'*' = 


~1 66 


(iii) 


star 


185847.7-370111° ^ 


ISO-CrA 177 


19 


02 


54 


64 


3Q 


46 


19 


2 


14.5±0.7 




11.9±0.6 




11.2±1.0 


1 




star+ dust (IRS) 


CrA-4107^" O 


CXO 23 


19 


01 


43 


07 


3Q 


50 


18 


,9 


0.549±0.032 


0.370±0.023 


0.279±0.054 


0.121±0.104'*' = 


0.178±0.230'*''= 


2 18 


f TTT1 


star 


lN2 


CXO 34 


19 


01 


55. 


77 


— 36 


57 


27, 


,5 


11.1±1.1 


41.3±2.7 


71.5±4.0 


66.5±5.0 




0.87 




cup-down 


1 


CXO 37 


19 


01 


57. 


40 


— 37 


03 


11, 


,8 


3.00±0.16 


2.31±0.12 


1.67±0.10 


1.03±0.07 


0.150±0.206'*-'= 


— 2.10 




star 


185834.7-370735° 


CXO 38 


19 


01 


58. 


32 


— 37 


00 


26, 


7 


3.25±0.16 


2.68±0.15 


2.19±0.14 


1.30±0.09 


0.045±0.258'*''= 


— 1 . 71 




star 


G-32^^ 185835.7-370450° 


CXO 42 


19 


02 


01. 


96 


— 36 


54 


00, 


,0 


0.292±0.026 


0.258±0.024 


0.217±0.140''' = 


1.04±0.10 


3.31±0.35 


0.38 


(I) 


red 


185839.6-365823° 


CrA-205 


19 


01 


11. 


72 


— 37 


22 


21, 


,5 




2.89±0.16 




1.06±0.09 


0.163±0.153'*' = 


— 2.49 




star 




CrA-452 


19 


00 


44. 


55 


— 37 


02 


11, 


,0 


24.4±1.2 


16.5±0.8 


11.7±0.6 


7.01±0.35 


0.951±0.236'* 


— 2.52 




star 


ISO-CrA 98^^ 


CrA-453 


19 


01 


04. 


62 


-37 


01 


29, 


4 


5.74±0.29 


4.31±0.22 


3.03±0.15 


1.82±0.10 


0.00174±0. 13000'*' = 


-2.36 


(III) 


star 


CXO 1^^ 


CrA-468 


19 


01 


49. 


,35 


-37 


00 


28, 


,6 


8.30±0.42 


5.94±0.29 


3.91±0.20 


2.32±0.12 


0.453±0.256'*' = 


-2.52 


(III) 


star 


CXO 26^^; LS-RCrA 2^^ 


CrA-4108 


19 


02 


09. 


,67 


-36 


46 


34, 


,5 


7.85±0.42 




3.88±0.20 




0.351±0.173'*' = 


-2.40 


(III) 


star 


ISO-CrA 141^2 


CrA-4111 


19 


01 


20. 


,85 


-37 


03 


02, 


,9 


4.16±0.22 


3.08±0.15 


2.16±0.12 


1.35±0.09 


0.568±0.186 


-2.29 


(III) 


star_no_MPl 


CXO 7^^ 


Haas 4 


19 


01 


40. 


,67 


-36 


56 


04, 


,9 


0.978±0.068 


1.41±0.09 


1.38±0.12 


1.77±0.29 


0.458±0.338'*' = 


-0.28 


(F) 


cup- down 




Haas 7 


19 


01 


43. 


,67 


-36 


57 


17, 


2 


0.477±0.032 


0.426±0.030 


0.284±0.071 


0.127±0.095'*' = 


<35'*'» 


-1.70 


(III) 


star 




Haas 17 


19 


01 


51. 


74 


-36 


55 


14, 


2 


0.724±0.066 


0.579±0.038 


0.484±0.081 


0.236±0.125'*' = 


<8.3'*- = 


-1.48 


(II) 


star 




Haas 33 


19 


02 


00. 


,85 


-36 


58 


50, 


,6 


2.02±0.11 


1.36±0.07 


1.11±0.09 


0.619±0.069 


<5'*' = 


-2.17 


(III) 


star 


185838.4-370312° 


Haas 34 


19 


02 


01. 


41 


-36 


58 


51, 


7 


0.685±0.042 


0.445±0.031 


0.407±0.065 


0.220±0.059 


<4'*- = 


-2.02 


(III) 


star 


185838.9-370314° 


LS-RCrA 1' 


19 


01 


33. 


,56 


-37 


00 


30, 


,3 


1.25±0.07 


1.20±0.06 


1.33±0.09 


2.01±0.12 


12.2±1.2 


-0.08 


(F) 


Galc_star-hdust(IR2) 


G09 CrA-18^^ 


G09 CrA-9 


19 


01 


58. 


,32 


-37 


01 


06. 


,0 


0.302±0.019 


0.453±0.028 


0.469±0.055 


0.336±0.051 


0.0422±0.24S0'*' = 


-0.55 


(II) 


cup-down 




G09 CrA-11 


19 


01 


01. 


11 


-36 


58 


34. 


5 


1.08±0.06 


0.777±0.042 


0.544±0.044 


0.162±0.06l'> 


0.0151±0.1610'*' = 


-2.38 


(III) 


star 




CP p> 


19 


01 


38. 


,91 


-36 


53 


26. 


,4 


583. ±33. 


366. ±23. 


272. ±50.'* 


296. ±74.'* 




-2.63 


(III) 


star 


HD 176386 


CP wa"* 


19 


02 


22. 


,43 


-36 


55 


38. 


,6 


14.2±l.l' 


10.3±0.5 


6.78±0.35 


4.14±0.23 


0.316±0.208'*- = 


-2.44 


(III) 


star 


CrAPMS 36^"^ 


LMlO-2 


19 


01 


39. 


11 


-37 


00 


17. 


,3 


6.24±0.32 


4.70±0.24 


3.33±0.19 


2.07±0.12 


0.317±0.224'*' = 


-2.19 


(III) 


star 


185816.6-370439° 


LMlO-3 


19 


01 


51. 


,80 


-37 


10 


47. 


,9 


14.8±1.0 


10.1±0.5 


6.77±0.34 


3.99±0.21 


0.088±0.20S'*' = 


-2.52 


(III) 


star 




LMlO-4 


19 


02 


07. 


,68 


-37 


11 


55. 


,9 


4.35±0.23 


3.14±0.16 


2.40±0.13 


1.33±0.08 


0.302±0.169'*-'= 


-2.30 


(III) 


star 




LMlO-5 


19 


02 


39. 


,40 


-36 


53 


11. 


2 


6.11±0.31 


4.23±0.21 


3.02±0.19 


1.74±0.10 


0.229±0.153'*''= 


-2.44 


(III) 


star 


ISO-CrA 163^^ 



^T he name listed in this column can be used t o search for the source in SIMB AD, except in the following cases: IRS jTaylor &; Storevlll984| ). Ha^^£a£raco_fc_^^^|^^|ren| |l98ll ) ■ CXO jForbrich &c PreibischI 
l2007l . listing number from their Table 2), CrA- jLopez Marti' et al1l2005l h Haas iHaas et al1l200g| . listing number from their Table 3); G09 CrA- j Gut ermiitl^^ni 20091 . listing number from their Table 



4); GP iGlass fc PenstQij[T973) ; LMIO- jLopez Marti et al.||2010l). There are 3 known YSOs that do not appear in this list because they do not have Spitzer counterparts. These sources are (SIMBAD 
coordinates in J200 0): IRS 9 {RA: 19 01 52.63, Dec: -36 57 00.2: lTavlor StQre^ll984l ') . SMA 2 (RS 9) (RA: 19 01 55.28, Dec: -36 57 16.6: lBrow'^ll9S7l^ . and Anon 2 (RA: 19 01 06.92, Dec: -36 58 07.3; 
iKnacke et al.||1973l ). 

'^Coordinates from our Spitzer catalog. 

''Re fer ences are as fo l lows: ^Taylor &: StorevI il984l^: HVrba et al.l il976al^: ■iNutter et al.l ^2005^: iKnacke et"aD il973l'): ^Haas et al.l f200i^: *] Glass &: Penstonl ^1975^: iHerbig &z Belli il988l'): fatten 
1991} ; ^Wilking et al.Nl997ll : ^ HLopez Marti et al.l ll2005l '>: ^ ^Sicilia-Aguilar et al.l i20o'l ): HOlofsson et al.l il999l '): ^ iForbrich &c Preibis Ji2007l l : ^ iFernandez Comeronl i200ll) : ^ iGutermuth et al. 



20091) : ^ Walter et al.Ml99:^) . 



'^The flux for this source at this wavelength has been "band-filled" (described in ^3l . 

^This flux is flagged "U" in the catalog indicating that it is an upper limit because it has been band-filled and has a S/N < 5. 

^Due to band-filling and saturation, the spectral index for R Cr A has b een derived from th e Kg and 5.8 fim bands alone. Its value and resulting clas sification arc likely incorrect. It should be noted that 
R CrA has long been known to be a PMS star with an accretion disk f e.g. [Knacke et al.|fl973l ). and an SED slope of a Class II jHillenbrand et al.iri992l ). 

^IRS 1 is the only source in this table that has a 70 fira flux in the catalog; it is an excellent "A" quality flux with a value of: 33700. it4700. 

^This flux has a quality of detection of "D" meaning it has a 2 < S/N < 3. Fluxes that do not have this note all have a quality of detection of "A" (S/N >7), "B" (S/N >5), or "C" (S/N > 3). 
^Discussed in lFernandez &c Comero 

JThis source is a double with CrA-54 (HD 176386B). 

This source name comes from IWilking et al it is located to the northeast of source GP w (called GP wb in the same paper, where GP is from [Glass ^^^ensto and was also seen by 

IWalter et al ] i 199711 who call the two sources CrAPMS 3 and 3B. CrAPMS 3 is also a YSO candidate, and can be found in Table [6] as CrA- 75. ft seems that in ^^eyer fc WilkinJ <2009l 'l there was some 
confusion when referring to GP wa - in that paper, it actually corresponds to our CrA-75. 

^This flux is flagged "K" in the catalog meaning that it is complex, i.e. there are two or more detections within 2" in the same band in the same epoch. The quality is likely still quite good, however. 
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Table 8. Clustering Analysis Comparison 



G09^ All YSO candidates Class I/II YSOs 
(35 sources) (116 sources) (62 sources) 



Urit 0.116 pc 0.164 pc 0.152 pc 

Ncores (# members) 1 (24) 1 (68) 1 (34) 

Rhuiiicore) 0.200 pc 0.385 pc 0.237 pc 

Rcirc{core) 0.317 pc 0.590 pc 0.326 pc 

Aspect Ratio^ 2.51 2.36 1.89 

Density^ 190 pc'^ 150 pc^^ ^gg pp-s 



'^Values from lGutermuth et al.l (120091 ). Note that they have been modified for 
a distance of 130 pc (a distance of 170 pc was used in their paper). 

Aspect Ratio = (Rcirc/R^n/z)^ 

■^Density = Number of members in each core / (vrx R^^^J. 
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Table 9. Submillimeter Array Results 



Source 


Other Name 


RA (J2000)^ 


Dec (J2000)^ 


-f^ 1.3mm 








h m s 


o / // 


[mJy] 


[Mo] 


CrA-20 


IRS 5N 


19 01 48.47 


-36 57 14.9 


95 (7) 


0.023 




R CrA 


19 01 53.65 


-36 57 08.6 


36 (25) 


0.009 




SMA2 


19 01 55.29 


-36 57 17.0 


269 (20) 


0.063 




IRS 7B 


19 01 56.40 


-36 57 28.3 


320 (19) 


0.076 


CrA-43 


WMB 55 


19 01 58.52 


-36 57 09.3 


135 (27) 


0.032 


CrA-44 


IRAS 32 


19 02 58.68 


-37 07 36.7 


103 (8) 


0.024 



^RA, Dec values from point source fits to the continuum image data. 

''Point source fit to continuum images using baselines longer than 20 kA. 

'^The uncertainties are reported in braces. These are the uncertainties 
returned by the fitting program IMFIT in MIRIAD. 

•^Disk mass assuming a dust temperature of 30 K. 



- 106 - 



Table 10. Known HH objects in CrA 



Knot R.A. Dec. Knot R.A. Dec. 

Name (J2000) (J2000) Name (J2000) (J2000) 



HH82A 


19 


01 


12. 


,480 


-36 


57 


17, 


,50 


HH729C 


19 01 15.679 


-36 57 56.41 


HH82B 


19 


01 


16. 


,090 


-36 


57 


29, 


,30 


HH730A 


19 01 25.26 


-36 59 30.0 


HH96 


19 


01 


41. 


,309 


-37 


00 


51, 


,08 


HH730B 


19 


01 


24.079 


-36 


59 


59.78 


HH97 


19 


01 


45, 


,780 


-37 


00 


01, 


,51 


HH730C 


19 


01 


15.720 


-37 


00 


29.99 


HH97E 


19 


01 


46, 


,200 


-37 


00 


03, 


,89 


HH731A 


19 


01 


40.030 


-36 


58 


57.29 


HH97W 


19 


01 


45, 


,749 


-37 


00 


02, 


,92 


HH731B 


19 


01 


38.270 


-36 


58 


45.59 


HH98 


19 


01 


52, 


,829 


-36 


57 


36, 


,22 


HH732A 


19 


01 


57.199 


-36 


52 


07.21 


HH98B 


19 


01 


52, 


,150 


-36 


57 


28, 


,40 


HH732B 


19 


01 


56.410 


-36 


51 


38.41 


HH99 


19 


02 


05, 


,304 


-36 


54 


39, 


,96 


HH732C 


19 


01 


57.730 


-36 


52 


17.00 


HH99A 


19 


02 


05, 


,410 


-36 


54 


39, 


,60 


HH733 


19 


02 


09.470 


-36 


54 


46.19 


HH99B 


19 


02 


05, 


,609 


-36 


54 


34, 


,99 


HH734A 


19 


02 


11.011 


-36 


58 


28.60 


HH99C 


19 


02 


00, 


,910 


-36 


55 


17, 


,51 


HH734B 


19 


02 


11.729 


-36 


58 


22.69 


HH99E 


19 


01 


59, 


,040 


-36 


57 


18, 


,00 


HH735A 


19 


02 


12.559 


-36 


56 


37.79 


HHIOO 


19 


01 


49, 


,104 


-36 


58 


15, 


,96 


HH735B 


19 


02 


12.761 


-36 


56 


42.11 


HHIOON 


19 


01 


49, 


,104 


-36 


58 


17, 


,04 


HH736 


19 


02 


14.410 


-36 


55 


22.19 


HHIOOS 


19 


01 


50, 


,208 


-36 


58 


.39, 


,00 


HH860 


19 


01 


39.401 


-36 


58 


03.79 


HHlOl 


19 


01 


42, 


,000 


-37 


00 


43, 


,20 


HH861A 


19 


02 


21.569 


-37 


16 


54.98 


HHIOIN 


19 


01 


35, 


,090 


-37 


02 


55, 


,39 


HH861B 


19 


02 


24.300 


-37 


15 


38.02 


HHIOIS 


19 


01 


34, 


,930 


-37 


03 


01, 


,19 


HH862 


19 


02 


36.581 


-37 


11 


57.19 


HH104 


19 


01 


58, 


,992 


-36 


57 


15, 


,84 


HH863A 


19 


03 


09.010 


-37 


23 


55.72 


HH104D 


19 


01 


45, 


,890 


-36 


57 


03, 


,71 


HH863B 


19 


03 


10.819 


-37 


19 


12.29 


HH104C 


19 


01 


45, 


,439 


-36 


56 


57, 


,59 


HH863C 


19 


03 


18.600 


-37 


23 


31.20 


HH729A 


19 


01 


10, 


,560 


-36 


57 


32, 


,40 


HH863D 


19 


03 


26.909 


-37 


21 


54.40 


HH729B 


19 


01 


12, 


,941 


-36 


57 


48, 


,20 


HH863E 


19 


03 


31.699 


-37 


20 


32.60 
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Table 11. Identified knots in the CrA star- forming region 



Knot 
Name 



R.A. 
(J2000) 



Dec. 
(J2000) 



Comments 



Detected in IRAC maps - EGOs 



Spitzer outflow 1 A 


19 


00 


24. 


,139 


-36 


58 


07, 


.20 


Part of 


a larger flow with 


outflow 


2 and 


MHO 


2005? 


Spitzer outflow 1 B 


19 


00 


24, 


,744 


-36 


58 


05, 


,19 


Part of 


a larger flow with 


outflow 


2 and 


MHO 


2005? 


Spitzer outflow 1 C 


19 


00 


27, 


,379 


-36 


59 


00, 


,99 


Part of 


a larger flow with 


outflow 


2 and 


MHO 


2005? 


Spitzer outflow 1 D 


19 


00 


28, 


,033 


-36 


59 


18, 


,27 


Part of 


a larger flow with 


outflow 


2 and 


MHO 


2005? 


Spitzer outflow 1 E 


19 


00 


28, 


,304 


-36 


58 


53, 


,48 


Part of 


a larger flow with 


outflow 


2 and 


MHO 


2005? 


Spitzer outflow 1 F 


19 


00 


28, 


,602 


-36 


59 


13, 


,01 


Part of 


a larger flow with 


outflow 


2 and 


MHO 


2005? 


Spitzer outflow 1 G 


19 


00 


29, 


,218 


-36 


59 


00, 


.20 


Part of 


a larger flow with 


outflow 


2 and 


MHO 


2005? 


Spitzer outflow 1 H 


19 


00 


30, 


,132 


-36 


59 


16, 


.10 


Part of 


a larger flow with 


outflow 


2 and 


MHO 


2005? 


Spitzer outflow 1 I 


19 


00 


30, 


,648 


-36 


59 


01, 


.46 


Part of 


a larger flow with 


outflow 


2 and 


MHO 


2005? 


Spitzer outflow 1 J 


19 


00 


31, 


,198 


-36 


59 


25, 


,28 


Part of 


a larger flow with 


outflow 


2 and 


MHO 


2005? 


Spitzer outflow 1 K 


19 


00 


31, 


,308 


-36 


59 


08, 


,80 


Part of 


a larger flow with 


outflow 


2 and 


MHO 


2005? 


Spitzer outflow 1 L 


19 


00 


31, 


,850 


-36 


59 


47, 


.26 


Part of 


a larger flow with 


outflow 


2 and 


MHO 


2005? 


Spitzer outflow 1 M 


19 


00 


31, 


,962 


-36 


59 


30, 


.18 


Part of 


a larger flow with 


outflow 


2 and 


MHO 


2005? 


Spitzer outflow 1 N 


19 


00 


32, 


,892 


-36 


59 


05, 


.16 


Part of 


a larger flow with 


outflow 


2 and 


MHO 


2005? 


Spitzer outflow 1 O 


19 


00 


33, 


,592 


-36 


59 


32, 


,03 


Part of 


a larger flow with 


outflow 


2 and 


MHO 


2005? 


Spitzer outflow 1 P 


19 


00 


33, 


,913 


-36 


59 


00, 


,30 


Part of 


a larger flow with 


outflow 


2 and 


MHO 


2005? 


Spitzer outflow 1 Q 


19 


00 


43, 


,098 


-36 


58 


59, 


,21 


Part of 


a larger flow with 


outflow 


2 and 


MHO 


2005? 


Spitzer outflow 2 A 


19 


01 


31, 


,574 


-36 


57 


03, 


,69 


Jet-like 


morphology, part 


of MHO 2005 






Spitzer outflow 2 B 


19 


01 


32, 


,635 


-36 


56 


55, 


.58 


Jet-like 


morphology, part 


of MHO 2005 






Spitzer outflow 2 C 


19 


01 


34, 


,272 


-36 


56 


56, 


.71 


Jet-like 


morphology, part 


of MHO 2005 






IRAS 32 A 


19 


03 


25, 


,207 


-37 


01 


07, 


,39 


Part of 


IRAS 32 outflow - 


NE, blue lobe 






IRAS 32 B 


19 


03 


22, 


159 


-37 


03 


53, 


,27 


Part of 


IRAS 32 outflow - 


NE, blue lobe 






IRAS 32 C 


19 


03 


06, 


,643 


-37 


05 


33, 


,85 


Part of 


IRAS 32 outflow - 


NE, blue lobe 






IRAS 32 H 


19 


02 


34, 


,231 


-37 


11 


19, 


,14 


Part of 


IRAS 32 outflow - 


SW, red lobe 






IRAS 32 I 


19 


02 


34, 


,128 


-37 


11 


39, 


,90 


Part of 


IRAS 32 outflow - 


SW, red lobe 






IRAS 32 J 


19 


02 


30. 


,290 


-37 


12 


46, 


,61 


Part of 


IRAS 32 outflow - 


SW, red lobe 






IRAS 32 K 


19 


02 


25. 


,219 


-37 


14 


14, 


,13 


Part of 


IRAS 32 outflow - 


SW, red lobe 






IRAS 32 L 


19 


02 


24, 


,948 


-37 


14 


20, 


,23 


Part of 


IRAS 32 outflow - 


SW, red lobe 






IRAS 32 M 


19 


02 


23, 


,923 


-37 


14 


42, 


,25 


Part of 


IRAS 32 outflow - 


SW, red lobe 






IRAS 32 N 


19 


02 


24, 


,840 


-37 


16 


05, 


,23 


Part of 


IRAS 32 outflow - 


SW, red lobe 







Detected in H2 and IRAC maps 



IRAS 32 D 
IRAS 32 E 
IRAS 32 F 
IRAS 32 G 
MHO 2000 A 
MHO 2000 B 
MHO 2001 A 
MHO 2001 B 
MHO 2001 G 



19 03 04.142 -37 06 21.01 MHO 2010 - Part of IRAS 32 outflow - NE, blue lobe 

19 03 04.092 -37 06 30.78 MHO 2010 - Part of IRAS 32 outflow - NE, blue lobe 

19 02 53.472 -37 08 47.56 MHO 2010 - Part of IRAS 32 outflow - SW, red lobe 

19 02 51.478 -37 09 01.72 MHO 2010 - Part of IRAS 32 outflow - SW, red lobe 

19 02 05.561 -36 54 37.60 HH 99 B - Driven by IRS 9 or GrA 24. 

19 02 04.668 -36 54 56.40 HH 99 A - Same flow of A? Possibly driven by IRS 6. 

19 02 01.085 -36 56 37.47 Driven by GXO 34? 

19 02 00.523 -36 56 40.25 Driven by GXO 34? 

19 01 59.873 -36 56 32.08 Driven by YSO in the SMA 2 region 
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Table 11 — Continued 



Knot 
Name 



R.A. 
(J2000) 



Dec. 
(J2000) 



Comments 



MHO 2001 


D 


19 


01 


59 


568 


-36 


56 


50 


98 


Driven by CXO 34? 


MHO 2001 


E 


19 


01 


58 


512 


-36 


56 


55 


32 


Driven by CXO 34? 


MHO 2001 


F 


19 


01 


57 


900 


-36 


56 


53 


74 


Driven by YSO in the SMA 2 region 


MHO 2001 


G 


19 


01 


57 


828 


-36 


57 


05 


84 


Driven by CrA-43 (SMM 2) 


MHO 2001 


H 


19 


01 


57 


660 


-36 


56 


59 


79 


Driven by CrA-43 (SMM 2)? 


MHO 2001 


I 


19 


01 


56 


626 


-36 


56 


57 


79 


Driven by YSO in the SMA 2 region 


MHO 2001 


J 


19 


01 


56 


412 


-36 


57 


17 


09 


Driven by SMA 2? 


MHO 2001 


K 


19 


01 


56 


110 


-36 


57 


02 


84 


Driven by YSO in the SMA 2 region 


MHO 2001 


L 


19 


01 


55 


792 


-36 


57 


11 


62 


Driven by YSO in the SMA 2 region 


MHO 2001 


M 


19 


01 


55 


475 


-36 


57 


18 


80 


Driven by IRS 7A? 


MHO 2001 


N 


19 


01 


55 


355 


-36 


57 


16 


21 


Driven by SMA 2? 


MHO 2001 


O 


19 


01 


55 


330 


-36 


57 


12 


20 


Driven by YSO in the SMA 2 region 


MHO 2001 


P 


19 


01 


55 


474 


-36 


57 


08 


74 


Driven by YSO in the SMA 2 region 


MHO 2001 


Q 


19 


01 


52 


520 


-36 


57 


16 


12 


Driven by SMA 2? 


MHO 2002 


A 


19 


01 


59 


462 


-36 


57 


17 


63 


Driven by SMA 2? 


MHO 2002 


B 


19 


01 


59 


222 


-36 


57 


14 


46 


Driven by CrA-43 (SMM 2) 


MHO 2002 


C 


19 


01 


58 


622 


-36 


57 


19 


08 


Driven by SMA 2 or IRS7B? 


MHO 2003 


A 


19 


01 


50 


081 


-36 


58 


30 


23 


Driven by IRS 1 


MHO 2003 


B 


19 


01 


50 


394 


-36 


58 


33 


69 


Driven by IRS 1 


MHO 2003 


C 


19 


01 


49 


620 


-36 


59 


10 


20 


Driven by IRS 1 


MHO 2004 


A 


19 


01 


45 


432 


-36 


56 


57 


74 


HH 104C - bow-shock ~east-west oriented - Part of a larger flow 


MHO 2004 


B 


19 


01 


45 


864 


-36 


57 


03 


79 


HH 104D - part of a curved jet ne-sw oriented. Driven by IRS 6 


MHO 2004 


C 


19 


01 


46 


849 


-36 


56 


53 


99 


Part of a curved jet ne-sw oriented. Driven by IRS 6 


MHO 2004 


D 


19 


01 


47 


929 


-36 


56 


47 


64 


Part of a curved jet ne-sw oriented. Driven by IRS 6 


MHO 2005 


A 


19 


01 


28 


879 


-36 


57 


10 


66 


Part of a larger east-west flow? 


MHO 2005 


B 


19 


01 


25 


702 


-36 


57 


21 


11 


Part of a larger east- west flow? 


MHO 2005 


C 


19 


01 


16 


202 


-36 


57 


28 


73 


HH 82B - Part of a larger east-west flow? 


MHO 2005 


D 


19 


01 


15 


790 


-36 


57 


26 


35 


Part of a larger east- west flow? 


MHO 2005 


E 


19 


01 


12 


038 


-36 


57 


08 


52 


Part of a larger east-west flow? 


MHO 2005 


F 


19 


01 


09 


670 


-36 


59 


23 


75 


Part of a larger east-west flow? 


MHO 2005 


G 


19 


01 


04 


706 


-36 


58 


09 


74 


Part of a larger east-west flow? 


MHO 2005 


H 


19 


01 


12 


426 


-36 


57 


17 


09 


HH 82A - Part of a larger east-west flow? 


MHO 2005 


I 


19 


00 


59 


969 


-36 


58 


15 


00 


Part of a larger east-west flow? 


MHO 2006 


A 


19 


01 


40 


923 


-36 


58 


47 


18 


Driven by IRS 2 


MHO 2006 


B 


19 


01 


40 


827 


-36 


58 


52 


08 


Driven by IRS 2 


MHO 2006 


C 


19 


01 


40 


033 


-36 


58 


57 


85 


Driven by IRS 2 


MHO 2006 


D 


19 


01 


33 


924 


-37 


01 


03 


54 


Driven by IRS 2 


MHO 2007 


A 


19 


01 


34 


937 


-37 


02 


55 


89 


HH lOlN - Driven by IRS 5? 


MHO 2007 


B 


19 


01 


34 


282 


-37 


03 


01 


41 


HH lOlS - Driven by IRS 5? 


MHO 2008 


A 


19 


01 


47 


393 


-36 


59 


52 


30 


Part of a larger outflow incl. MHO 2003? 


MHO 2008 


B 


19 


01 


47 


198 


-37 


00 


02 


71 


Part of a larger outflow incl. MHO 2003? 


MHO 2008 


C 


19 


01 


46 


214 


-37 


00 


05 


50 


HH 97E - Part of a larger outflow incl. MHO 2003? 


MHO 2008 


D 


19 


01 


40 


080 


-37 


00 


02 


97 


Bow-shock ENE-WSW oriented, hkely driven by CrA-41. 


MHO 2009 


A 


19 


02 


12 


343 


-36 


57 


10 


80 


Driven by SMA 2? 
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Table 11 — Continued 



Knot 




R.A. 




Dec. 


Comments 


Name 




(.12000) 


(.12000) 




MHO 2009 


B 


19 02 


11.357 


-36 


57 09.12 


Driven by SMA 2? 


MHO 2009 


C 


19 02 


11.031 


-36 


57 09.67 


Driven by SMA 2? 


MHO 2009 


D 


19 02 


10.786 


-36 


57 12.80 


Driven by SMA 2? 


MHO 2009 


E 


19 02 


09.240 


-36 


57 13.97 


Driven by SMA 2? 


MHO 2009 


F 


19 02 


06.463 


-36 


57 23.72 


Driven by SMA 2? 


MHO 2009 


G 


19 02 


04.601 


-36 


57 10.32 


Driven by SMA 2? 


MHO 2011 


A 


19 02 


21.247 


-36 


55 53.97 


Part of a large west-east precessing outflow 


MHO 2011 


B 


19 02 


19.037 


-36 


56 11.80 


Part of a large west-east precessing outflow 


MHO 2011 


C 


19 02 


17.796 


-36 


56 24.24 


Part of a large west-east precessing outflow 


MHO 2011 


D 


19 02 


12.622 


-36 


56 43.94 


HH 735B - Part of a large west-east precessing outflow 


MHO 2011 


E 


19 02 


11.002 


-36 


56 49.38 


Part of a large west-east precessing outflow 


MHO 2011 


F 


19 02 


10.284 


-36 


56 52.56 


Part of a large west-east precessing outflow 


MHO 2012 


A 


19 02 


13.980 


-36 


56 05.99 


Driven by IRS 7B? 


MHO 2012 


B 


19 02 


07.356 


-36 


56 36.37 


Driven by IRS 7B? 


MHO 2012 


C 


19 02 


06.360 


-36 


56 41.61 


Driven by IRS 7B? 


MHO 2013 


A 


19 02 


09.447 


-36 


54 41.50 


Driven by CXO 34? 


MHO 2013 


B 


19 02 


09.454 


-36 


54 49.02 


HH 733 - Driven by CXO 34? 


MHO 2013 


C 


19 02 


09.221 


-36 


54 52.45 


Driven by CXO 34? 


MHO 2014 


A 


19 02 


00.641 


-36 


52 27.15 


Driven by IRS 5? 


MHO 2014 


B 


19 02 


00.281 


-36 


52 20.94 


Driven by IRS 5? 


MHO 2014 


C 


19 01 


59.801 


-36 


52 22.11 


HH 732C - Driven by IRS 5? 


MHO 2014 


D 


19 01 


57.662 


-36 


52 13.50 


Driven by IRS 5? 


MHO 2014 


E 


19 01 


57.158 


-36 


52 04.57 


HH 732A - Driven by IRS 5? 


MHO 2014 


F 


19 01 


56.916 


-36 


52 17.29 


Driven by IRS 5? 
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Table 12. Measured H2 proper motions in the CrA knots 



Knot 


P.M. 


^tan 


P.A. 


± 'I CXi.L±l.\j 


(-// -n 


I Ja.±±± J 


(°) 
K ) 


MHO 2000 A 


0.079±0.009 


49±6 


43±6 


MHO 2000 B 


0.024±0.007 


15±4 


91±16 


MHO 2001 A 


0.065±0.007 


40±5 


81±5 


MHO 2001 B 


0.065±0.007 


40±5 


81±5 


MHO 2001 C 


0.040±0.010 


24±6 


43±15 


MHO 2001 D 


0.074±0.010 


46±7 


54±8 


MHO 2001 E 


O.llOiO.OlO 


70±7 


41±5 


MHO 2001 F 


0.090±0.015 


55±10 


45±10 


MHO 2001 G 


0.039±0.007 


24±4 


328±10 


MHO 2001 J 


0.175±0.006 


108±4 


79±2 


MHO 2001 K 


0.067±0.009 


41±6 


45±8 


MHO 2001 M 


0.084±0.018 


52±11 


97±8 


MHO 2001 N 


0.023±0.009^ 


14±6^ 


62±23 


MHO 2001 


0.044±0.012 


27±7 


59±14 


MHO 2001 P 


O.llOiO.OlO 


70±7 


40±5 


MHO 2001 Q 


0.058±0.008 


36±5 


289±8 


MHO 2002 A 


0.083±0.014 


52±9 


81±9 


MHO 2002 B 


0.056±0.008 


35±5 


133±9 


MHO 2002 C 


0.040±0.008 


25±5 


65±12 


MHO 2003 A 


0.110±0.013 


68±8 


199±5 


MHO 2003 B 


0.056±0.013 


34±8 


220±13 


MHO 2003 C 


0.186±0.028 


115±17 


225±9 


MHO 2004 A 


0.038±0.011 


24±7 


286±11 


MHO 2006 C 


0.075±0.012 


46±8 


243±9 


MHO 2007 A 


0.073±0.012 


45±8 


215±10 


MHO 2007 B 


0.097±0.012 


60±8 


201±7 


MHO 2008 C 


0.085±0.012 


53±7 


215±8 


MHO 2008 D 


0.071±0.012 


44±8 


261±10 


MHO 2009 D 


0.147±0.045 


91±28 


63±15 


MHO 2009 E 


0.137±0.018 


85±11 


66±7 
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Table 12— Continued 



Knot 


P.M. 


Vtan 


P.A. 


Name 




(kms-1) 


(°) 


MHO 2009 F 


0.127±0.019 


79±12 


80±8 


MHO 2009 G 


0.147±0.036 


91±22 


63±12 


MHO 2011 D 


0.068±0.018 


42±11 


68±15 


MHO 2011 E 


0.256±0.019 


158±12 


84±4 


MHO 2011 F 


0.080±0.019 


49±12 


76±13 


MHO 2012 B 


0.134±0.018 


83±11 


68±8 


MHO 2012 C 


0.140±0.018 


87±11 


62±7 


MHO 2013 A 


0.086±0.017 


53±11 


45±11 


MHO 2013 B 


0.041±0.017^ 


26±1P 


101±24 


MHO 2013 C 


0.091±0.017 


56±11 


63±11 


MHO 2014 A 


0.065±0.017 


40±11 


32±15 


MHO 2014 B 


0.133±0.017 


82±11 


34±7 


MHO 2014 C 


0.071±0.017 


44±11 


31±14 


MHO 2014 F 


0.037±0.017^ 


23±11^ 


347±27 



^Below 3(7 



- 112 - 



Table 13. Candidate driving sources and possible flows 



YSO 


Flow* 


RA.(°) - blue lobe 


IRSl 


MHO 2003, MHO 2008 A-C (b) 


203 


IRS 2 


MHO 2006 (r?) 


26 


IRS 5N (?) 


MHO 2007 (b); MHO 2014 (r) 


25 


IRS 6 


MHO 2000 B (?), HH99C (?) (r); MHO 2004 B-D (b?) 


240 


IRS 9 or CrA-24 


MHO 2000 A (r) 


229 or 223 


IRS 7A (?) 


MHO 2001 F (?) 


55 


IRS 7B 


MHO 2012, HH736 (?), MHO 2002 C (r?); HH731B (?) (r?) 


248 


CXO 34 (?) 


MHO 2013 (?), MHO 2001 A, B, E (?), F (?) (b?) 


45 


SMA 2 region 


MHO 2013 (?), MHO 2001 O, P, K, I (?) (b?); MHO 2001 Q (r?) 


~48 


SMA2 


MHO 2009, MHO2001J, A (b?); MHO 2001 Q (r?) 


90 


IRS7A or SMA 2 


MHO 2011 (b); MHO 2004 D; Sptz. outflow 1, 2, MHO 2005 (r) 


~90 


CrA-43 (SMM 2) 


MHO 2002 B, HH734 (?) (b?); HH MHO 2001 G, H(?) (r?) 


125 


IRAS 32 


IRAS 32 outflow A-E (b); F-N, HH861 (r) 


39 


CrA-41 


MHO 2008 D (b?) 


242 


CrA-40 (?) 


HH 730 A-C (?) 


201 (?) 


S CrA (?) 


HH 762 (b?) 


115 


W CrA (?) 


HH862 (b?) 


274 



''b = blue lobe, r = red lobe 



